Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,986 A→T 7.5% intergenic (+49/+27) 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,9860AT 7.5% 586.3 / 28.7 174intergenic (+49/+27)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 = 8.55e-01

Only 100 of 178 total aligned reads displayed.

GATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCA  >  NC_000913/1018838‑1019133
                                                                                                                                                     |                                                                                                                                                    
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‑AAAAAAAACCCCGCTGCTGCggg                                                                                                                                                 <  1:937494/149‑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                                                                                                                                             <  3:1531989/149‑1 (MQ=255)
         gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttt                                                                                                                                             >  2:349221/1‑149 (MQ=255)
         gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttt                                                                                                                                             >  1:673229/1‑149 (MQ=255)
            gATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTc                                                                                                                                          >  2:912436/1‑149 (MQ=255)
            gATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTc                                                                                                                                          <  2:845143/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)
               tGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAc                                                                                                                                       <  3:5820/149‑1 (MQ=255)
                gAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                                      >  3:156314/1‑149 (MQ=255)
                 aaGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCa                                                                                                                                     <  3:13392/149‑1 (MQ=255)
                        ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATATCTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACgaga                                                                                                                              >  4:1028137/1‑149 (MQ=255)
                        ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACgaga                                                                                                                              >  2:407048/1‑149 (MQ=255)
                         gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                             <  4:1190329/149‑1 (MQ=255)
                          aaaaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAAc                                                                                                                            >  1:898196/1‑149 (MQ=255)
                              cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                        <  3:1471398/149‑1 (MQ=255)
                              cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                        <  4:295888/149‑1 (MQ=255)
                               aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAg                                                                                                                       >  3:913035/1‑149 (MQ=255)
                               aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAg                                                                                                                       >  2:1120296/1‑149 (MQ=255)
                               aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAg                                                                                                                       >  2:126958/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)
                                               aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTAc                                                                                                      >  2:1292209/1‑149 (MQ=255)
                                                  aaaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAAc                                                                                                   <  2:357225/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                                                                                                                      <  1:832787/124‑1 (MQ=255)
                                                         ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                      >  2:832787/1‑124 (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)
                                                                  gCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTa                                                                                   >  2:632773/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)
                                                                              ggCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATctct                                                                       >  4:370009/1‑149 (MQ=255)
                                                                               gCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTa                                                                      >  1:651819/1‑149 (MQ=255)
                                                                                 aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTc                                                                                                                                          >  2:183635/1‑80 (MQ=255)
                                                                                 aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcg                                                                    >  1:113760/1‑149 (MQ=255)
                                                                                   tGCTGGGTAAGGAATAACTGACGAAAGTCGGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcgcg                                                                  <  3:1110662/149‑1 (MQ=255)
                                                                                    gCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcga                                                                 >  3:219432/1‑149 (MQ=255)
                                                                                          tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                      <  3:448148/91‑1 (MQ=255)
                                                                                          tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTa                                                                      >  1:750934/1‑138 (MQ=255)
                                                                                          tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTa                                                                      <  2:750934/138‑1 (MQ=255)
                                                                                          tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑‑‑AAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCgg                                                         >  1:499604/1‑149 (MQ=255)
                                                                                           aaGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                                                           >  3:25102/1‑132 (MQ=255)
                                                                                           aaGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                                                           <  4:25102/132‑1 (MQ=255)
                                                                                             ggAATAACTGACGCAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCgaa                                                        >  4:1435082/1‑147 (MQ=255)
                                                                                             ggAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                                        <  3:1270051/149‑1 (MQ=255)
                                                                                              gAATAACTGACGAAAGTCAGTTCAGTTTACTAAAGGC‑AAAAAAAACGCTGCACGAGCGGGGCTTTTTGACCGGACGAG‑CTTTAAGCTCGCGGTGTGGGTTACTACTTCTTTGTTTCGTTTAACTCGGATCTCGACGCGGCGTTccctgg                                                       >  4:1152655/1‑145 (MQ=255)
                                                                                               aaTAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                                      >  3:1287734/1‑149 (MQ=255)
                                                                                                 tAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCa                                                    >  4:670735/1‑149 (MQ=255)
                                                                                                  aaCTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                                   <  3:1146530/149‑1 (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)
                                                                                                             gTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                        <  1:229320/149‑1 (MQ=255)
                                                                                                              tCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAg                                       >  4:680092/1‑149 (MQ=255)
                                                                                                                 gTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgc                                    >  2:1553118/1‑149 (MQ=255)
                                                                                                                         tACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCt                            <  4:1467367/149‑1 (MQ=255)
                                                                                                                           cTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACTAGAACTTAAGCCTGCGA‑CTGGGTCACAACGTCTTTGATACCGTTAACTTCGATCTAGACGTGAAGATCGGGAGCGAGGCAGTCGATCAGTGCAGCACGCTGt                          >  2:515870/1‑149 (MQ=255)
                                                                                                                            tAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGtt                         >  2:1044795/1‑149 (MQ=255)
                                                                                                                            tAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGtt                         >  3:26646/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)
                                                                                                                                        aaaaCCCCGCTGCTGCGGGGTTTCTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                                                           <  2:586193/88‑1 (MQ=39)
                                                                                                                                        aaaaCCCCGCTGCTGCGCGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTca              >  1:456456/1‑149 (MQ=255)
                                                                                                                                         aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTaa                                                                                  >  4:1176016/1‑80 (MQ=255)
                                                                                                                                          acccccGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGAGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcacc            >  2:1176573/3‑148 (MQ=255)
                                                                                                                                          aaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGCGCCGGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca            >  4:216052/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)
                                                                                                                                                    cagcGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCa  <  2:673229/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGG‑CTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCA  >  NC_000913/1018838‑1019133

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