Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,288 A→G 8.4% intergenic (+64/+78) PP_1501 → / ← PP_1502 hypothetical protein/OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,707,2880AG 8.4% 205.4 / 7.0 71intergenic (+64/+78)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base A (30/35);  new base G (3/3);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

AGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCG  >  NC_002947/1707139‑1707437
                                                                                                                                                     |                                                                                                                                                      
agAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                                        <  1:935749/150‑1 (MQ=255)
  aaaGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGAGCGAGCGCCGCGCGGGCGGCGCTCGCTCTcaca                                                                                                                                                      >  2:1874354/1‑150 (MQ=255)
          gACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGcac                                                                                                                                              <  1:133066/150‑1 (MQ=255)
           aCCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGcaca                                                                                                                                             >  2:881906/1‑150 (MQ=255)
             ccgCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCacaac                                                                                                                                           >  1:1673802/3‑150 (MQ=255)
                   gCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCCCACGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATc                                                                                                                                     <  1:645236/150‑1 (MQ=255)
                       gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCg                                                                                                                                 >  2:2148755/1‑150 (MQ=255)
                       gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCg                                                                                                                                 >  1:1761207/1‑150 (MQ=255)
                          gACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAAc                                                                                                                              >  2:1517187/1‑150 (MQ=255)
                            cAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACaa                                                                                                                            >  1:1851855/1‑150 (MQ=255)
                             ccccGAAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCTGCACAACATTATCGGCGAACaaa                                                                                                                           <  1:2238499/148‑1 (MQ=255)
                               ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                                      >  1:1128828/1‑121 (MQ=255)
                               ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                                      <  2:1128828/121‑1 (MQ=255)
                                    aGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGc                                                                                                                    <  2:691363/150‑1 (MQ=255)
                                     gcgcgcGGTGCGCTGAATGGGAGAGCAAACAGCGAGATACATACATTGAGGAGTGGGCGCATCATGGTGGGCGCCCGCGAGACCGAGCGCCGCGCAGGCCGCACTCGATCCAACAAGCTTCACAACATTATCGGCGAACAAATAAAATcc                                                                                                                   <  2:1673802/150‑1 (MQ=255)
                                                 cTGCATGGGCCAGGCAAAACACCGCAAACTGCAGTCAGCTATCGGACTAGAATTTTTGCCTCTCGTGATATCGTGTGTCTCTCGGGCTTCGCTCGTTCTCACAGGTGTCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATcg                                                                                                      <  1:2089267/150‑1 (MQ=255)
                                                  tGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATcgc                                                                                                     <  1:112643/150‑1 (MQ=255)
                                                       gggCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGc                                                                                                <  1:1259415/150‑1 (MQ=255)
                                                              gAAACACCGCGGCAGGTGCAGTTAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                                         >  2:1668981/1‑150 (MQ=255)
                                                              gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGACGCGGGGCttttctt                                                                                         <  2:2059438/150‑6 (MQ=38)
                                                              gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                                         <  1:665698/150‑1 (MQ=255)
                                                              gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGACTGTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                                         <  1:1874354/150‑1 (MQ=255)
                                                               aaaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCtt                                                                                        <  2:1851855/150‑1 (MQ=255)
                                                                aaCACCGCGGCAGGTGCAGTTAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCttt                                                                                       <  1:1668981/150‑1 (MQ=255)
                                                                 acacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTc                                                                                      >  2:1103933/1‑150 (MQ=255)
                                                                   acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                                         <  1:389044/145‑1 (MQ=255)
                                                                   acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                                         >  2:389044/1‑145 (MQ=255)
                                                                   acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAg                                                                                    >  2:645742/1‑150 (MQ=255)
                                                                      gcgGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTa                                                                                 <  1:2148755/150‑1 (MQ=255)
                                                                             gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCggg                                                                                                  <  2:2197713/126‑1 (MQ=255)
                                                                             gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCggg                                                                                                  >  1:2197713/1‑126 (MQ=255)
                                                                                 aGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAg                                                                      >  2:823955/1‑150 (MQ=255)
                                                                                  gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgc                                                                     >  1:1899315/1‑150 (MQ=255)
                                                                                   tGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTCTCAGCTAGCCAGCCTTAgcg                                                                    >  1:1304931/1‑150 (MQ=255)
                                                                                     aGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCt                                                                  >  1:1597687/1‑150 (MQ=255)
                                                                                     aGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCt                                                                  <  2:436907/150‑1 (MQ=255)
                                                                                      ctcgcgggcgccaaaaatgctgcgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                                                  <  2:1698703‑M2/60‑1 (MQ=255)
                                                                                      ctcgcgggcgccaaaaatgctgcgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                                                  >  1:1698703‑M2/27‑86 (MQ=255)
                                                                                       tcgcgggcgccaaaaatgctgcgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                                      >  1:2319903‑M2/26‑65 (MQ=255)
                                                                                       tcgcgggcgccaaaaatgctgcgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                                      <  2:2319903‑M2/40‑1 (MQ=255)
                                                                                       gTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                                   >  1:284580/1‑115 (MQ=255)
                                                                                       gTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                                   <  2:284580/115‑1 (MQ=255)
                                                                                         ttCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAAc                                                              <  1:823955/150‑1 (MQ=255)
                                                                                          tCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCGTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAAcc                                                             >  1:1410728/1‑150 (MQ=255)
                                                                                             GCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                                          >  1:881906‑M1/1‑150 (MQ=255)
                                                                                                 agcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGCCGCCAAAGCGATCCAGTCTTGGCGCTTAGCCTTGCTCa                                                      >  2:610843/1‑150 (MQ=25)
                                                                                                 agcaTTGTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCa                                                      <  2:2281086/150‑1 (MQ=255)
                                                                                                    aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTc                                                                                           <  2:2102110/110‑1 (MQ=255)
                                                                                                    aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTc                                                                                           >  1:2102110/1‑110 (MQ=255)
                                                                                                    aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCg                                                   >  2:1507341/1‑150 (MQ=255)
                                                                                                            cgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                           <  2:1761207‑M2/146‑1 (MQ=255)
                                                                                                            cgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                           <  1:802076‑M2/146‑1 (MQ=255)
                                                                                                            cgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                           <  1:1024164‑M2/146‑1 (MQ=255)
                                                                                                                   gagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCg                                    <  2:1595949/150‑1 (MQ=255)
                                                                                                                    agaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCgg                                   <  2:1604251/150‑1 (MQ=255)
                                                                                                                     gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGc                                  <  1:1507341/150‑1 (MQ=255)
                                                                                                                     gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGc                                  <  1:1846131/150‑1 (MQ=255)
                                                                                                                       tCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGc                                                                                                >  1:1865755/1‑86 (MQ=37)
                                                                                                                       tCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGc                                                                                                <  2:1865755/86‑1 (MQ=37)
                                                                                                                         gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAg                                                                                                                          >  2:2286385/1‑59 (MQ=255)
                                                                                                                         gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAg                                                                                                                          <  1:2286385/59‑1 (MQ=255)
                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                                   <  1:1375518/81‑1 (MQ=25)
                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                                   >  2:1375518/1‑81 (MQ=25)
                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt                              >  1:1509367/1‑150 (MQ=255)
                                                                                                                          aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTCTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt                             <  2:1304931/150‑1 (MQ=255)
                                                                                                                            cgccgcCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                           <  1:1621589/134‑1 (MQ=255)
                                                                                                                            cgccgcCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                           >  2:1621589/1‑134 (MQ=255)
                                                                                                                             gccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                           >  1:1013685/1‑133 (MQ=255)
                                                                                                                             gccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                           <  2:1013685/133‑1 (MQ=255)
                                                                                                                                  gcgggcggCGCTCGATCTCACAGGCGGCACACCATAGTTAGAGGCCGATGAAAAGCCCCGCGTCGATCGCGCGACTGGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTgg                     <  1:1103933/150‑1 (MQ=255)
                                                                                                                                         gcgcTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGt              <  2:1509367/150‑1 (MQ=255)
                                                                                                                                           gcacGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCa            >  1:2329358/4‑150 (MQ=255)
                                                                                                                                            cTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCAc           >  2:1878249/1‑150 (MQ=255)
                                                                                                                                             tCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCAcc          >  1:1945740/1‑150 (MQ=255)
                                                                                                                                               gATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCt         >  1:1489539/1‑149 (MQ=255)
                                                                                                                                               gATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCt         <  2:1489539/149‑1 (MQ=255)
                                                                                                                                                 tctcACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCa      >  1:1410267/1‑150 (MQ=255)
                                                                                                                                                     acaGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCg  >  1:2192859/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCG  >  NC_002947/1707139‑1707437

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: