Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,983 A→T 7.4% intergenic (+46/+30) matP → / ← ompA Ter macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,9830AT 7.4% 594.8 / 27.8 174intergenic (+46/+30)matP/ompATer macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (97/64);  new base T (6/7);  total (103/71)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 177 total aligned reads displayed.

GATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGC  >  NC_000913/1018838‑1019131
                                                                                                                                                  |                                                                                                                                                     
gATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCgcagca                                                                                                                                                    >  3:773513/1‑149 (MQ=255)
    gTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCgg                                                                                                                                                <  1:1590793/149‑1 (MQ=255)
     ttaaGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCggg                                                                                                                                               <  4:1088109/146‑1 (MQ=255)
     ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCggg                                                                                                                                               >  4:75879/1‑149 (MQ=255)
     ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCggg                                                                                                                                               <  3:1099345/149‑1 (MQ=255)
     ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCggg                                                                                                                                               <  3:125285/149‑1 (MQ=255)
        aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                            <  2:1339637/149‑1 (MQ=255)
         gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttt                                                                                                                                           >  1:673229/1‑149 (MQ=255)
           tGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                                         <  3:948878/149‑1 (MQ=255)
              ttGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGtttttttt                                                                                                                                      <  3:1216860/149‑4 (MQ=255)
              ttGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTa                                                                                                                                      >  1:1459704/1‑149 (MQ=255)
              ttGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACACCGCTGCTGCGGGGtttttttt                                                                                                                                      <  1:969185/149‑4 (MQ=255)
                gAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                                    >  2:35512/1‑149 (MQ=255)
                gAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                                    >  3:156314/1‑149 (MQ=255)
                 aaGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCa                                                                                                                                   <  3:13392/149‑1 (MQ=255)
                       cGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACgag                                                                                                                             >  2:1535607/1‑149 (MQ=255)
                        ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATATCTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACgaga                                                                                                                            >  4:1028137/1‑149 (MQ=255)
                        ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACgaga                                                                                                                            >  2:407048/1‑149 (MQ=255)
                          aaaaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAAc                                                                                                                          >  1:898196/1‑149 (MQ=255)
                           aaaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACt                                                                                                                         >  2:927455/1‑149 (MQ=255)
                              cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                      <  4:295888/149‑1 (MQ=255)
                              cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                      <  3:1471398/149‑1 (MQ=255)
                               aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAg                                                                                                                     >  2:1120296/1‑149 (MQ=255)
                                aaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                    >  4:1226705/1‑149 (MQ=255)
                                   agaAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTg                                                                                                                 >  3:1518485/1‑149 (MQ=255)
                                    gaAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGc                                                                                                                <  1:349221/149‑1 (MQ=255)
                                          cgcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGa                                                                                                         <  2:804313/149‑1 (MQ=255)
                                             gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGtt                                                                                                      >  3:265254/1‑149 (MQ=255)
                                               aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTAc                                                                                                    >  2:1292209/1‑149 (MQ=255)
                                                tAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACa                                                                                                   >  3:1301981/1‑149 (MQ=255)
                                                 aaaaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACaa                                                                                                  <  2:747184/149‑1 (MQ=255)
                                                     tGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCggg                                                                                                                                               >  3:247180/1‑101 (MQ=255)
                                                       tcttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGAGGGGTTTTTCTACCAGACGAGAACTTAAACCTGCGG‑CTGAGTTACAACGTCt                                                                                             >  4:636204/1‑149 (MQ=255)
                                                        cttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCttt                                                                                           <  4:1522348/149‑1 (MQ=255)
                                                         ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                    >  2:832787/1‑124 (MQ=255)
                                                         ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                    <  1:832787/124‑1 (MQ=255)
                                                         ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTg                                                                                          >  2:1267203/1‑149 (MQ=255)
                                                          tctcTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGa                                                                                         >  2:285842/1‑149 (MQ=255)
                                                          tctcTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGa                                                                                         <  2:1459704/149‑1 (MQ=255)
                                                            tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                      >  4:230165/1‑119 (MQ=255)
                                                            tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATa                                                                                       >  4:1554662/1‑149 (MQ=255)
                                                              tGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGt                                                                                                       >  1:1334611/1‑131 (MQ=255)
                                                              tGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATAcc                                                                                     >  3:704827/1‑149 (MQ=255)
                                                               gAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCt                                                                                    <  1:1535607/149‑1 (MQ=255)
                                                                aaGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCtt                                                                                   >  4:1525255/1‑149 (MQ=255)
                                                                aaGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCtt                                                                                   <  3:289870/149‑1 (MQ=255)
                                                                 aGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCttt                                                                                  <  1:927455/149‑1 (MQ=255)
                                                                  gCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTa                                                                                 >  2:632773/1‑149 (MQ=255)
                                                                  gCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTa                                                                                 >  3:1284857/1‑149 (MQ=255)
                                                                  gCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTa                                                                                 >  2:944393/1‑149 (MQ=255)
                                                                    aGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAAc                                                                               >  1:1347491/1‑149 (MQ=255)
                                                                     ggATCTGCAGGCATCGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTACGG‑CTGAGTTACAACGTCTTTGATACCTTTAACt                                                                              >  4:860901/1‑149 (MQ=255)
                                                                           gCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCCCATCAGTGGGGTTTTTCAACTAACCAAGACCTTCAGCCGGGGC‑TTGTATTACAACGTCTTGGTAACTTTTAACTCCGact                                                                        >  3:551057/1‑147 (MQ=255)
                                                                            cAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtct                                                                       <  2:984454/149‑1 (MQ=255)
                                                                             aGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtctc                                                                      >  1:287196/1‑149 (MQ=255)
                                                                               gCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTa                                                                    >  2:893564/1‑149 (MQ=255)
                                                                               gCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTa                                                                    >  1:651819/1‑149 (MQ=255)
                                                                                 aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTc                                                                                                                                        <  1:183635/80‑1 (MQ=255)
                                                                                 aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTc                                                                                                                                        >  2:183635/1‑80 (MQ=255)
                                                                                 aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcg                                                                  >  1:690980/1‑149 (MQ=255)
                                                                                 aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcg                                                                  >  1:113760/1‑149 (MQ=255)
                                                                                    gCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcga                                                               >  3:219432/1‑149 (MQ=255)
                                                                                          tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                    <  3:448148/91‑1 (MQ=255)
                                                                                          tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑‑‑AAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCgg                                                       >  1:499604/1‑149 (MQ=255)
                                                                                           aaGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                                                         <  4:25102/132‑1 (MQ=255)
                                                                                             ggAATAACTGACGCAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCgaa                                                      >  4:1435082/1‑147 (MQ=255)
                                                                                              gAATAACTGACGAAAGTCAGTTCAGTTTACTAAAGGC‑AAAAAAAACGCTGCACGAGCGGGGCTTTTTGACCGGACGAG‑CTTTAAGCTCGCGGTGTGGGTTACTACTTCTTTGTTTCGTTTAACTCGGATCTCGACGCGGCGTTccctgg                                                     >  4:1152655/1‑145 (MQ=255)
                                                                                                 tAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCa                                                  >  4:670735/1‑149 (MQ=255)
                                                                                                  aaCTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTACGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                                 <  3:860901/149‑1 (MQ=255)
                                                                                                   aCTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAgg                                                >  1:261863/1‑149 (MQ=255)
                                                                                                     tGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                              <  2:261863/149‑1 (MQ=255)
                                                                                                     tGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                              <  4:1518485/149‑1 (MQ=255)
                                                                                                       aCGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGt                                            >  2:522890/1‑149 (MQ=255)
                                                                                                        cGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                           <  4:1221266/149‑1 (MQ=255)
                                                                                                         gAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                           <  3:636204/149‑1 (MQ=255)
                                                                                                         gAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCg                                          >  4:404965/1‑149 (MQ=255)
                                                                                                           aaGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAt                                        <  3:1532386/149‑1 (MQ=255)
                                                                                                             gTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                      >  2:759482/1‑149 (MQ=255)
                                                                                                              tCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAg                                     >  4:680092/1‑149 (MQ=255)
                                                                                                                           cTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACTAGAACTTAAGCCTGCGA‑CTGGGTCACAACGTCTTTGATACCGTTAACTTCGATCTAGACGTGAAGATCGGGAGCGAGGCAGTCGATCAGTGCAGCACGCTGt                        >  2:515870/1‑149 (MQ=255)
                                                                                                                            tAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGtt                       >  2:1044795/1‑149 (MQ=255)
                                                                                                                            tAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCCTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGtt                       >  4:814174/1‑149 (MQ=255)
                                                                                                                             aaaGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt                      <  4:423245/149‑1 (MQ=255)
                                                                                                                              aaGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTc                     <  4:462433/149‑1 (MQ=255)
                                                                                                                                ggC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCAc                   <  4:1296645/149‑1 (MQ=255)
                                                                                                                                ggC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCAc                   <  1:1044795/149‑1 (MQ=255)
                                                                                                                                    aaaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt                <  3:370009/149‑1 (MQ=255)
                                                                                                                                     aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTTAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTg               >  1:1037541/1‑149 (MQ=255)
                                                                                                                                       aaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGAAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTc             >  4:208970/1‑149 (MQ=255)
                                                                                                                                        aaaaCCCCGCTGCTGCGGGGTTTCTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                                                         >  1:586193/1‑88 (MQ=39)
                                                                                                                                        aaaaCCCCGCTGCTGCGCGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTca            >  1:456456/1‑149 (MQ=255)
                                                                                                                                         aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTaa                                                                                >  4:1176016/1‑80 (MQ=255)
                                                                                                                                         aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTaa                                                                                <  3:1176016/80‑1 (MQ=39)
                                                                                                                                          acccccGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGAGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcacc          >  2:1176573/3‑148 (MQ=255)
                                                                                                                                          aaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGCGCCGGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca          >  4:216052/1‑149 (MQ=255)
                                                                                                                                          aaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca          >  3:1129816/1‑149 (MQ=255)
                                                                                                                                            ccccGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTTCAGCACGCTGTTTCACGTTGTCACAgg        >  2:1567374/1‑149 (MQ=255)
                                                                                                                                            ccccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAgg        <  1:1267203/149‑1 (MQ=255)
                                                                                                                                              ccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgtg      >  1:845826/1‑149 (MQ=255)
                                                                                                                                                  agcagcGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGc  <  2:502183/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGC  >  NC_000913/1018838‑1019131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑