Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,539,877 T→C 2.7% I231V (ATC→GTC)  narY ← nitrate reductase 2 (NRZ) subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,539,8770T→C 2.7% 5.6 ‑0.4 ‑0.0223I231V (ATC→GTC) narYnitrate reductase 2 (NRZ) subunit beta
Reads supporting (aligned to +/- strand):  new base (2/4):  ref base (105/112):  total (107/116)
Fisher's exact test strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.84e-01

Only 100 of 242 total aligned reads displayed.

CCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCC  >  NC_007779/1539777‑1539976
                                                                                                    |                                                                                                    
cccACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAt                                                                                                      >  1:1692954/1‑101 (MQ=255)
cccACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAt                                                                                                      <  1:9110466/101‑1 (MQ=255)
   aCGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAtttt                                                                                                   >  1:4136130/1‑101 (MQ=255)
   aCGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAtttt                                                                                                   <  1:5489907/101‑1 (MQ=255)
   aCGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAtttt                                                                                                   <  1:7193288/101‑1 (MQ=255)
    cGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAttttt                                                                                                  <  1:4741335/101‑1 (MQ=255)
     gCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAtttttt                                                                                                 >  1:4521408/1‑101 (MQ=255)
     gCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAtttttt                                                                                                 >  1:5107347/1‑101 (MQ=255)
      cAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAttttttt                                                                                                <  1:8570078/101‑1 (MQ=255)
      cAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGAttttttt                                                                                                >  1:7474924/1‑101 (MQ=255)
       aGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTg                                                                                               <  1:8878928/101‑1 (MQ=255)
        ggTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGt                                                                                              <  1:2686365/101‑1 (MQ=255)
          tttCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTAc                                                                                            >  1:10055236/1‑101 (MQ=255)
          tttCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTAc                                                                                            >  1:10988995/1‑101 (MQ=255)
           ttCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACg                                                                                           <  1:14124267/101‑1 (MQ=255)
            tCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACgg                                                                                          <  1:3211525/101‑1 (MQ=255)
             cTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGa                                                                                         <  1:11694019/101‑1 (MQ=255)
              tGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGAc                                                                                        >  1:13251714/1‑101 (MQ=255)
                   acacGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCg                                                                                   <  1:11504254/101‑1 (MQ=255)
                      cggtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCtt                                                                               >  1:113572/1‑101 (MQ=255)
                       ggtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTa                                                                              >  1:6898939/1‑101 (MQ=255)
                        gtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTAt                                                                             <  1:7825044/101‑1 (MQ=255)
                          cggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGc                                                                           <  1:14327844/101‑1 (MQ=255)
                           ggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGca                                                                          <  1:2972993/101‑1 (MQ=255)
                           ggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGca                                                                          <  1:5985391/101‑1 (MQ=255)
                            gtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGcac                                                                         <  1:10573470/101‑1 (MQ=255)
                             ttGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACGACCGCTTATGcaca                                                                        >  1:12479706/1‑101 (MQ=255)
                             ttGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGcaca                                                                        >  1:11576920/1‑101 (MQ=255)
                             ttGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGcaca                                                                        >  1:11915564/1‑101 (MQ=255)
                               ggcctggcTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaaa                                                                      <  1:3839989/94‑1 (MQ=255)
                               gACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaaa                                                                      >  1:6063417/1‑101 (MQ=255)
                                aCCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAAc                                                                     >  1:6277722/1‑101 (MQ=255)
                                aCCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAAc                                                                     <  1:12905889/101‑1 (MQ=255)
                                 ccGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACg                                                                    <  1:3228974/101‑1 (MQ=255)
                                  cGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgc                                                                   <  1:4326277/101‑1 (MQ=255)
                                    gACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgcca                                                                 <  1:4660191/101‑1 (MQ=255)
                                    gACTCAATTCGTGGGCAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgcca                                                                 <  1:12180229/101‑1 (MQ=255)
                                     aCTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgccag                                                                >  1:12798481/1‑101 (MQ=255)
                                      cTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTTGACCGACCACCCCTTTTTGCACAAAGgcaac                                                                >  1:2800967/1‑98 (MQ=21)
                                         aaTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCAc                                                            >  1:4997089/1‑101 (MQ=255)
                                          aTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACg                                                           <  1:11764157/101‑1 (MQ=255)
                                           ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                          <  1:3389598/101‑1 (MQ=255)
                                           ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                          >  1:7235817/1‑101 (MQ=255)
                                             cGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCa                                                        <  1:9763837/101‑1 (MQ=255)
                                             cGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCa                                                        <  1:10753340/101‑1 (MQ=255)
                                               tGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTAGGACAACCGCCTTTTCACCAACCGCAACCACCGCCttt                                                     >  1:515364/1‑101 (MQ=21)
                                                gggTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCAttt                                                     >  1:89772/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)
                                                    aaCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCGAGCCACGGCATTTATcc                                                 >  1:10098287/1‑101 (MQ=255)
                                                    aaCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATcc                                                 >  1:13801777/1‑101 (MQ=255)
                                                     aCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCt                                                >  1:3568924/1‑101 (MQ=255)
                                                     aCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCt                                                >  1:2443656/1‑101 (MQ=255)
                                                      cAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTg                                               >  1:2550117/1‑101 (MQ=255)
                                                        gaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGat                                             <  1:12634630/101‑1 (MQ=255)
                                                         aagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGatc                                            >  1:7568314/1‑101 (MQ=255)
                                                              gCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGatcaatca                                       <  1:2194598/101‑1 (MQ=255)
                                                               cATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAg                                      <  1:1599641/101‑1 (MQ=255)
                                                               cATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTAGCCTGATCAATCAg                                      <  1:7112681/101‑1 (MQ=255)
                                                                  ttttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcac                                   >  1:7665387/1‑101 (MQ=255)
                                                                   tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcaca                                  >  1:13297577/1‑101 (MQ=255)
                                                                   tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcaca                                  >  1:6381007/1‑101 (MQ=255)
                                                                    ttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACaa                                 >  1:13928117/1‑101 (MQ=255)
                                                                    ttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACaa                                 <  1:1131575/101‑1 (MQ=255)
                                                                    ttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACaa                                 <  1:5346523/101‑1 (MQ=255)
                                                                     tCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAAt                                >  1:11928571/1‑101 (MQ=255)
                                                                     tCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAAt                                <  1:12879291/101‑1 (MQ=255)
                                                                       tGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGc                              >  1:10454366/1‑101 (MQ=255)
                                                                          cTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGt                           >  1:13552269/1‑101 (MQ=255)
                                                                          cTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGt                           >  1:1081269/1‑101 (MQ=255)
                                                                            tGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtct                         <  1:7456114/101‑1 (MQ=255)
                                                                             gCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGTCACGGCATTTATCCTGATCAATCAGCACAATGCCGtctt                        >  1:9698414/1‑101 (MQ=255)
                                                                               cGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtcttct                      >  1:4342558/1‑101 (MQ=255)
                                                                                ggtcttgtggTTGAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGAGGCAGCGGTCCTGGTCGACCAGCACAATGCCGtcttctg                     >  1:13728520/10‑100 (MQ=255)
                                                                                 cTTTTCCAGTTGAAGTAGATTTTTTTGTCCGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                    <  1:5650597/101‑1 (MQ=255)
                                                                                 cTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                    <  1:13485819/101‑1 (MQ=255)
                                                                                 cTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                    <  1:9416577/101‑1 (MQ=255)
                                                                                  ttttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcg                   <  1:9540493/101‑1 (MQ=255)
                                                                                  ttttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcg                   <  1:6570632/101‑1 (MQ=255)
                                                                                  ttttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcg                   >  1:1979732/1‑101 (MQ=255)
                                                                                   cttgtggTTGAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGtcttctgtgc                  >  1:5474064/7‑97 (MQ=255)
                                                                                    ttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTcgcg                 <  1:2840017/101‑1 (MQ=255)
                                                                                      ccAGTTGAAGTAGATTTTTTTGTACTGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGtt               >  1:313728/1‑101 (MQ=255)
                                                                                      ccAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGtt               <  1:8234762/101‑1 (MQ=255)
                                                                                        aGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTg             <  1:6826055/101‑1 (MQ=255)
                                                                                          ttGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTa           >  1:5047113/1‑101 (MQ=255)
                                                                                            gAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGa         <  1:1754699/101‑1 (MQ=255)
                                                                                            gAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGa         >  1:6496252/1‑101 (MQ=255)
                                                                                            gAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTaca         <  1:5889810/101‑3 (MQ=255)
                                                                                            gAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTaca         <  1:3350705/101‑3 (MQ=255)
                                                                                              agtttttttttttGTAGGGACACCCTTTTATGCACAACCGCCAGCCACGTCTTTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATa       <  1:9198226/95‑1 (MQ=255)
                                                                                              aGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATa       <  1:1506551/101‑1 (MQ=255)
                                                                                               gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg      <  1:1314175/101‑1 (MQ=255)
                                                                                               gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg      >  1:13804491/1‑101 (MQ=255)
                                                                                               gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg      <  1:10243207/101‑1 (MQ=255)
                                                                                               gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg      <  1:14152649/101‑1 (MQ=255)
                                                                                               gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg      <  1:1563382/101‑1 (MQ=255)
                                                                                                 aGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgcg    >  1:8759608/1‑101 (MQ=255)
                                                                                                  gATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgcgc   <  1:13637418/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCC  >  NC_007779/1539777‑1539976

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.