Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,539,917 A→G 4.0% R217R (CGT→CGC narY ← nitrate reductase 2 (NRZ) subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,539,9170A→G 4.0% 12.8 ‑1.7 ‑0.0223R217R (CGT→CGCnarYnitrate reductase 2 (NRZ) subunit beta
Reads supporting (aligned to +/- strand):  new base (2/7):  ref base (105/109):  total (107/116)
Fisher's exact test strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

Only 100 of 242 total aligned reads displayed.

CAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTT  >  NC_007779/1539817‑1540017
                                                                                                    |                                                                                                    
cAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgccagcca                                                                                                      >  1:9092514/1‑101 (MQ=255)
 aaTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCAc                                                                                                     >  1:10220984/1‑101 (MQ=255)
  aTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACg                                                                                                    >  1:4179311/1‑101 (MQ=255)
  aTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACg                                                                                                    <  1:11764157/101‑1 (MQ=255)
   ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                                                                   <  1:3389598/101‑1 (MQ=255)
   ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                                                                   <  1:1382296/101‑1 (MQ=255)
     cGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCa                                                                                                 <  1:9763837/101‑1 (MQ=255)
     cGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCa                                                                                                 <  1:5535264/101‑1 (MQ=255)
       tGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTAGGACAACCGCCTTTTCACCAACCGCAACCACCGCCttt                                                                                              >  1:515364/1‑101 (MQ=21)
          gTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTAt                                                                                            >  1:10553962/1‑101 (MQ=255)
          gTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTAt                                                                                            >  1:11470523/1‑101 (MQ=255)
          gTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTAt                                                                                            >  1:125750/1‑101 (MQ=255)
           tAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATc                                                                                           <  1:6700310/101‑1 (MQ=255)
            aaCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATcc                                                                                          <  1:11102179/101‑1 (MQ=255)
            aaCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATcc                                                                                          >  1:13801777/1‑101 (MQ=255)
             aCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCt                                                                                         >  1:2443656/1‑101 (MQ=255)
              cAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTg                                                                                        >  1:12291687/1‑101 (MQ=255)
               agaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGa                                                                                       <  1:13980475/101‑1 (MQ=255)
                gaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGat                                                                                      >  1:7344888/1‑101 (MQ=255)
                gaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGat                                                                                      <  1:12634630/101‑1 (MQ=255)
                gaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGat                                                                                      <  1:12918031/101‑1 (MQ=255)
                   gaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGatcaa                                                                                   <  1:13163858/101‑1 (MQ=255)
                      gCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGatcaatca                                                                                <  1:14151333/101‑1 (MQ=255)
                      gCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGatcaatca                                                                                <  1:2194598/101‑1 (MQ=255)
                       cATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTAGCCTGATCAATCAg                                                                               <  1:7112681/101‑1 (MQ=255)
                        aTTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGc                                                                              <  1:11820388/101‑1 (MQ=255)
                         tttttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGca                                                                             <  1:8200898/101‑1 (MQ=255)
                          ttttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcac                                                                            >  1:7665387/1‑101 (MQ=255)
                           tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcaca                                                                           >  1:13297577/1‑101 (MQ=255)
                           tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcaca                                                                           >  1:8927182/1‑101 (MQ=255)
                           tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcaca                                                                           >  1:6381007/1‑101 (MQ=255)
                             tCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAAt                                                                         >  1:11928571/1‑101 (MQ=255)
                              cTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATg                                                                        >  1:7810532/1‑101 (MQ=255)
                               tGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGc                                                                       >  1:2595791/1‑101 (MQ=255)
                               tGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGc                                                                       <  1:5068123/101‑1 (MQ=255)
                               tGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGc                                                                       >  1:10454366/1‑101 (MQ=255)
                                  cTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGt                                                                    >  1:13552269/1‑101 (MQ=255)
                                  cTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGt                                                                    <  1:1136220/101‑1 (MQ=255)
                                  cTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGt                                                                    >  1:12501475/1‑101 (MQ=255)
                                  cTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGt                                                                    >  1:1081269/1‑101 (MQ=255)
                                   ttGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtc                                                                   >  1:14269469/1‑101 (MQ=255)
                                   ttGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtc                                                                   >  1:5216053/1‑101 (MQ=255)
                                    tGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtct                                                                  >  1:3125971/1‑101 (MQ=255)
                                     gCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGTCACGGCATTTATCCTGATCAATCAGCACAATGCCGtctt                                                                 >  1:9698414/1‑101 (MQ=255)
                                       cGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtcttct                                                               >  1:4342558/1‑101 (MQ=255)
                                       cGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtcttct                                                               >  1:4379240/1‑101 (MQ=255)
                                        ggtcttgtggTTGAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGAGGCAGCGGTCCTGGTCGACCAGCACAATGCCGtcttctg                                                              >  1:13728520/10‑100 (MQ=255)
                                        gCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTcttctt                                                              <  1:1155546/101‑1 (MQ=255)
                                         cTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                                                             <  1:4999503/101‑1 (MQ=255)
                                         cTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                                                             <  1:9416577/101‑1 (MQ=255)
                                          ttttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcg                                                            >  1:10679353/1‑101 (MQ=255)
                                          ttttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcg                                                            >  1:1979732/1‑101 (MQ=255)
                                          ttttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcg                                                            <  1:6570632/101‑1 (MQ=255)
                                          ttttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcg                                                            <  1:9540493/101‑1 (MQ=255)
                                            ttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcgcg                                                          <  1:7519384/101‑1 (MQ=255)
                                             tCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGt                                                         >  1:5936822/1‑101 (MQ=255)
                                                aGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTg                                                      <  1:6826055/101‑1 (MQ=255)
                                                    gAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGa                                                  <  1:1754699/101‑1 (MQ=255)
                                                    gAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGa                                                  >  1:6496252/1‑101 (MQ=255)
                                                    gAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTaca                                                  <  1:3350705/101‑3 (MQ=255)
                                                    gAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTaca                                                  <  1:5889810/101‑3 (MQ=255)
                                                    gAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTaca                                                  <  1:5359370/101‑3 (MQ=255)
                                                     aaGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGAt                                                 <  1:10085256/101‑1 (MQ=255)
                                                      agtttttttttttGTAGGGACACCCTTTTATGCACAACCGCCAGCCACGTCTTTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATa                                                <  1:9198226/95‑1 (MQ=255)
                                                      aGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATa                                                <  1:1506551/101‑1 (MQ=255)
                                                       gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg                                               <  1:1563382/101‑1 (MQ=255)
                                                       gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg                                               >  1:13804491/1‑101 (MQ=255)
                                                        tAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgc                                              >  1:31701/1‑101 (MQ=255)
                                                         aGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgcg                                             <  1:7471357/101‑1 (MQ=255)
                                                          gATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgcgc                                            <  1:13637418/101‑1 (MQ=255)
                                                           aTTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGcgcc                                           <  1:5479907/101‑1 (MQ=255)
                                                             ttttttGTACGGACAACCGCTTATGCACAAACGCCGGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGcgccgc                                         >  1:9333625/1‑101 (MQ=255)
                                                               ttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCtt                                       <  1:4650319/101‑1 (MQ=255)
                                                                cttGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGcgcccgatg                                      <  1:8572484/100‑6 (MQ=255)
                                                                 ttGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGcgcccgatgg                                     <  1:8981800/101‑7 (MQ=255)
                                                                      cGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAgg                                >  1:3861324/1‑101 (MQ=255)
                                                                        gACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTc                              <  1:8612613/101‑1 (MQ=255)
                                                                        gACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTc                              <  1:3467626/101‑1 (MQ=255)
                                                                           aaCCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGcc                           >  1:2844596/1‑101 (MQ=255)
                                                                            aCCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCa                          <  1:3531098/101‑1 (MQ=255)
                                                                             ccGCTTATGCACAAACGCCAGCCACTGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCAc                         >  1:6517619/1‑101 (MQ=255)
                                                                              cGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACg                        >  1:957351/1‑101 (MQ=255)
                                                                              cGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACg                        >  1:11718987/1‑101 (MQ=255)
                                                                                cTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCa                      <  1:8108694/101‑1 (MQ=255)
                                                                                cTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATACCGCCGCTTGGGCAGGTCGCCACGCa                      <  1:12734826/101‑1 (MQ=255)
                                                                                  tATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGc                    >  1:11915170/1‑101 (MQ=255)
                                                                                     gCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTgg                 <  1:7970805/101‑1 (MQ=255)
                                                                                      cacaAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTggg                >  1:6655268/1‑101 (MQ=255)
                                                                                          aaCGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTg            >  1:6112883/1‑101 (MQ=255)
                                                                                           aCGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGa           <  1:1675496/101‑1 (MQ=255)
                                                                                            cGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAg          <  1:12875141/101‑1 (MQ=255)
                                                                                            cGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAg          <  1:6969919/101‑1 (MQ=255)
                                                                                             gccagccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAgg         <  1:3639028/101‑1 (MQ=255)
                                                                                             gccagccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAgg         >  1:13503584/1‑101 (MQ=255)
                                                                                             gccagccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAgg         >  1:10737130/1‑101 (MQ=255)
                                                                                               cagccaCGGCATTTGTCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCa       <  1:10067514/101‑1 (MQ=255)
                                                                                               cagccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCa       >  1:4739777/1‑101 (MQ=255)
                                                                                                agccaCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAg      >  1:6418460/1‑101 (MQ=255)
                                                                                                agccaCGGCATTTATCCTGATCAACCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAg      >  1:13477661/1‑101 (MQ=255)
                                                                                                    aCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGtt  >  1:4277330/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTT  >  NC_007779/1539817‑1540017

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.