Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,539,887 C→G 3.8% P227P (CCG→CCC narY ← nitrate reductase 2 (NRZ) subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,539,8870C→G 3.8% 14.5 ‑1.8 ‑0.0238P227P (CCG→CCCnarYnitrate reductase 2 (NRZ) subunit beta
Reads supporting (aligned to +/- strand):  new base (2/7):  ref base (113/114):  total (117/121)
Fisher's exact test strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.95e-01

Only 100 of 247 total aligned reads displayed.

TTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGG  >  NC_007779/1539787‑1539987
                                                                                                    |                                                                                                     
tttCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTAc                                                                                                       >  1:10055236/1‑101 (MQ=255)
  tCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACgg                                                                                                     <  1:9343275/101‑1 (MQ=255)
  tCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACgg                                                                                                     <  1:3211525/101‑1 (MQ=255)
   cTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGa                                                                                                    <  1:3904680/101‑1 (MQ=255)
   cTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGa                                                                                                    <  1:11694019/101‑1 (MQ=255)
   cTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGa                                                                                                    >  1:1327955/1‑101 (MQ=255)
   cTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGa                                                                                                    <  1:13407168/101‑1 (MQ=255)
         acacGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCg                                                                                              <  1:11504254/101‑1 (MQ=255)
            cggtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCtt                                                                                          >  1:113572/1‑101 (MQ=255)
             ggtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTa                                                                                         >  1:4899961/1‑101 (MQ=255)
             ggtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTa                                                                                         >  1:6898939/1‑101 (MQ=255)
             ggtcggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTa                                                                                         <  1:749781/101‑1 (MQ=255)
                cggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGc                                                                                      >  1:2087585/1‑101 (MQ=255)
                 ggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGca                                                                                     <  1:5985391/101‑1 (MQ=255)
                 ggtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGca                                                                                     <  1:2972993/101‑1 (MQ=255)
                  gtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGcac                                                                                    <  1:10573470/101‑1 (MQ=255)
                  gtTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGcac                                                                                    >  1:2657761/1‑101 (MQ=255)
                   ttGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGcaca                                                                                   >  1:11915564/1‑101 (MQ=255)
                   ttGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGcaca                                                                                   >  1:11576920/1‑101 (MQ=255)
                    tGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTATGGACAACCGCTTATGCACaa                                                                                  >  1:7171134/1‑101 (MQ=255)
                    tGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaa                                                                                  <  1:4327354/101‑1 (MQ=255)
                     gACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaaa                                                                                 >  1:6063417/1‑101 (MQ=255)
                     gACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACaaa                                                                                 >  1:13783005/1‑101 (MQ=255)
                      aCCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAAc                                                                                >  1:6277722/1‑101 (MQ=255)
                       ccGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACg                                                                               >  1:12622390/1‑101 (MQ=255)
                       ccGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACg                                                                               <  1:3228974/101‑1 (MQ=255)
                        cGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgc                                                                              >  1:1689349/1‑101 (MQ=255)
                         ggACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgcc                                                                             >  1:145630/1‑101 (MQ=255)
                          gACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgcca                                                                            <  1:2460418/101‑1 (MQ=255)
                           aCTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgccag                                                                           >  1:12798481/1‑101 (MQ=255)
                            cTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTTGACCGACCACCCCTTTTTGCACAAAGgcaac                                                                           >  1:2800967/1‑98 (MQ=21)
                            cTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgccagc                                                                          >  1:512189/1‑101 (MQ=255)
                             tCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACgccagcc                                                                         >  1:9654894/1‑101 (MQ=255)
                                aTTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACg                                                                      >  1:4179311/1‑101 (MQ=255)
                                 ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                                     <  1:3389598/101‑1 (MQ=255)
                                 ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                                     <  1:1382296/101‑1 (MQ=255)
                                 ttCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACgg                                                                     >  1:7235817/1‑101 (MQ=255)
                                   cGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCa                                                                   <  1:10753340/101‑1 (MQ=255)
                                      gggTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTTCGGACAATCGCTTATGCACAAACGCCAGCCACGGCAttt                                                                >  1:12094714/1‑101 (MQ=255)
                                        gTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTAt                                                              >  1:11470523/1‑101 (MQ=255)
                                         tAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATc                                                             >  1:13641602/1‑101 (MQ=255)
                                         tAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATc                                                             <  1:6700310/101‑1 (MQ=255)
                                          aaCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCGAGCCACGGCATTTATcc                                                            >  1:10098287/1‑101 (MQ=255)
                                          aaCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATcc                                                            <  1:11102179/101‑1 (MQ=255)
                                           aCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCt                                                           >  1:4428646/1‑101 (MQ=255)
                                           aCAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCt                                                           >  1:3568924/1‑101 (MQ=255)
                                            cAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTg                                                          >  1:12291687/1‑101 (MQ=255)
                                             agaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGa                                                         <  1:13980475/101‑1 (MQ=255)
                                              gaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGat                                                        <  1:12918031/101‑1 (MQ=255)
                                              gaagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGat                                                        <  1:912884/101‑1 (MQ=255)
                                               aagaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGatc                                                       >  1:7568314/1‑101 (MQ=255)
                                                     cATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAg                                                 <  1:1599641/101‑1 (MQ=255)
                                                      aTTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGc                                                <  1:11820388/101‑1 (MQ=255)
                                                      aTTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGc                                                <  1:1360444/101‑1 (MQ=255)
                                                         tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcaca                                             >  1:8927182/1‑101 (MQ=255)
                                                         tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGcaca                                             >  1:6381007/1‑101 (MQ=255)
                                                          ttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACaa                                            >  1:13928117/1‑101 (MQ=255)
                                                           tCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAAt                                           <  1:12879291/101‑1 (MQ=255)
                                                            cTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATg                                          >  1:7810532/1‑101 (MQ=255)
                                                             tGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGc                                         >  1:10454366/1‑101 (MQ=255)
                                                                cTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGt                                      <  1:13802265/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)
                                                                 ttGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtc                                     <  1:9479365/101‑1 (MQ=255)
                                                                  tGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGtct                                    <  1:7456114/101‑1 (MQ=255)
                                                                       cTTTTCCAGTTGAAGTAGATTTTTTTGTCCGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                               <  1:5650597/101‑1 (MQ=255)
                                                                       cTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                               <  1:13485819/101‑1 (MQ=255)
                                                                       cTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCttcttc                               <  1:4999503/101‑1 (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:9224065/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:5359370/101‑3 (MQ=255)
                                                                                  gAAGTAGACCTTCTTGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTaca                    <  1:3350705/101‑3 (MQ=255)
                                                                                     gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg                 <  1:1314175/101‑1 (MQ=255)
                                                                                     gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg                 <  1:1563382/101‑1 (MQ=255)
                                                                                     gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg                 <  1:10243207/101‑1 (MQ=255)
                                                                                     gTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAg                 <  1:14152649/101‑1 (MQ=255)
                                                                                     gTAAATTTTTTTGTACGGCAACCCGCTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgc                <  1:2427278/101‑1 (MQ=25)
                                                                                       aGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAgcg               >  1:8759608/1‑101 (MQ=255)
                                                                                         aTTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGcgcc             <  1:5479907/101‑1 (MQ=255)
                                                                                         aTTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGcgcc             >  1:5103381/1‑101 (MQ=255)
                                                                                            tttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCt          <  1:9754391/101‑1 (MQ=255)
                                                                                             ttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCtt         <  1:4650319/101‑1 (MQ=255)
                                                                                             ttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGCCTTCTTCGCGTTTGTAGATAGCGCCGCtt         <  1:4285866/101‑1 (MQ=255)
                                                                                              tttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTg        >  1:13542449/1‑101 (MQ=255)
                                                                                              cttGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGcgcccgatg        <  1:8572484/100‑6 (MQ=255)
                                                                                               ttGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGcgcccgatgg       <  1:8981800/101‑7 (MQ=255)
                                                                                               ttGTAGGGGCAACCGGAGACGCACATGCGCCAACCGCGGCAGCGGTCCTGGTCGACCAGCACAATGCCGTCTTCTGTGCGTTTGTACATCGcgcccgatgg       <  1:6584987/101‑7 (MQ=255)
                                                                                               ttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTgg       <  1:9785547/101‑1 (MQ=255)
                                                                                                tGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTggg      >  1:794569/1‑101 (MQ=255)
                                                                                                 gTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGc     <  1:12603235/101‑1 (MQ=255)
                                                                                                 gTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGc     >  1:6529080/1‑101 (MQ=255)
                                                                                                    cGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAgg  <  1:3741152/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
TTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGG  >  NC_007779/1539787‑1539987

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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.