Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,853,676 A→T 5.4% intergenic (+34/‑246) tisB → / → emrD toxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,853,6760AT 5.4% 1020.7 / 30.8 291intergenic (+34/‑246)tisB/emrDtoxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (158/117);  new base T (8/8);  total (166/125)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 302 total aligned reads displayed.

GTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATA  >  NC_000913/3853527‑3853824
                                                                                                                                                     |                                                                                                                                                    
gTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGGCGGGaa                                                                                                                                                      >  4:2324283/1‑149 (MQ=255)
    tCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAAcc                                                                                                                                                   >  4:1124684/1‑149 (MQ=255)
       tGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCt                                                                                                                                                >  2:1983748/1‑149 (MQ=255)
         aaCCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTcc                                                                                                                                              <  4:1785856/149‑1 (MQ=255)
            cacaGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCt                                                                                                                                           >  4:765268/1‑149 (MQ=255)
              caGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCttt                                                                                                                                         <  3:2277360/149‑1 (MQ=255)
                 gagaCGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTtat                                                                                                                                      >  4:105685/1‑149 (MQ=255)
                   gaCGCGTATGAACCTGGTGGATATCTCCATTCTTATCCGCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGGCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTatat                                                                                                                                    <  1:112907/149‑1 (MQ=255)
                    aCGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATAtt                                                                                                                                   >  2:1138327/1‑149 (MQ=255)
                      gcgTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTGTATTTa                                                                                                                                 <  2:1785675/149‑1 (MQ=255)
                      gcgTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTa                                                                                                                                 <  4:1862919/149‑1 (MQ=255)
                         tATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCg                                                                                                                              >  1:204036/1‑149 (MQ=255)
                          aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGt                                                                                                                             <  3:258054/149‑1 (MQ=255)
                          aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGt                                                                                                                             <  3:1836278/149‑1 (MQ=255)
                          aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGt                                                                                                                             >  4:2145908/1‑149 (MQ=255)
                          aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGt                                                                                                                             >  2:899905/1‑149 (MQ=255)
                            gAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGc                                                                                                                           >  3:881678/1‑149 (MQ=255)
                            gAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGc                                                                                                                           >  1:1302739/1‑149 (MQ=255)
                              aCCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCgg                                                                                                                         <  3:655369/149‑1 (MQ=255)
                              aCCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCgg                                                                                                                         >  1:1355804/1‑149 (MQ=255)
                                   gtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTgg                                                                                                                    >  1:413147/1‑149 (MQ=255)
                                       atatCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTc                                                                                                                <  4:825428/149‑1 (MQ=255)
                                       atatCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTc                                                                                                                <  1:1522529/149‑1 (MQ=255)
                                       atatCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTc                                                                                                                <  2:1194341/149‑1 (MQ=255)
                                        tatCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCg                                                                                                               <  4:517103/149‑1 (MQ=255)
                                           cGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgat                                                                                                            <  1:250181/149‑1 (MQ=255)
                                           cGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgat                                                                                                            >  3:2444989/1‑149 (MQ=255)
                                           cGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgat                                                                                                            >  4:2374568/1‑149 (MQ=255)
                                            gCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAAc                                                                                                                                                    >  2:2325132/1‑108 (MQ=255)
                                            gCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCt                                                                                                                                           <  3:1329166/117‑1 (MQ=255)
                                            gCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgatg                                                                                                           >  3:729800/1‑149 (MQ=255)
                                            gCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTCTTATATTTAGCGTGCGGGTTGGTGTCGgatg                                                                                                           >  3:2434998/1‑149 (MQ=255)
                                             ccATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgatgc                                                                                                          >  3:1999135/1‑149 (MQ=255)
                                              cATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaa                                                                                                                                                      >  4:396419/1‑104 (MQ=255)
                                                     aTCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGa                                                                                                  <  3:1279276/149‑1 (MQ=255)
                                                       ccTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACg                                                                                                >  1:595071/1‑149 (MQ=255)
                                                        cTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGc                                                                                               >  1:645067/1‑149 (MQ=255)
                                                        cTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGc                                                                                               <  3:781086/149‑1 (MQ=255)
                                                         tCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGTTTCCCGCTCCCCttt                                                                                                                                         <  4:1319298/106‑1 (MQ=38)
                                                         tCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCa                                                                                              >  4:2395817/1‑149 (MQ=255)
                                                          cAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCAt                                                                                             <  4:2444989/149‑1 (MQ=255)
                                                          cAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCAt                                                                                             <  1:1029888/149‑1 (MQ=255)
                                                              cTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTa                                                                                         >  2:1230080/1‑149 (MQ=255)
                                                               tCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                     >  1:2441683/1‑88 (MQ=255)
                                                               tCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                     <  2:2441683/88‑1 (MQ=255)
                                                               tCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTAt                                                                                        >  3:568098/1‑149 (MQ=255)
                                                               tCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTAt                                                                                        >  4:330315/1‑149 (MQ=255)
                                                                 aTTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATcc                                                                                      >  1:797205/1‑149 (MQ=255)
                                                                 aTTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTCTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATcc                                                                                      >  2:710610/1‑149 (MQ=255)
                                                                 aTTGCTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATcc                                                                                      >  2:1355438/1‑149 (MQ=255)
                                                                  ttgttgCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCg                                                                                     >  2:2392146/1‑149 (MQ=255)
                                                                    gttgCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGcc                                                                                   >  4:1978500/1‑149 (MQ=255)
                                                                    gttgCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGcc                                                                                   >  1:1456164/1‑149 (MQ=255)
                                                                     ttgCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATACGATGCTGACGCATCTTATCCGccc                                                                                  >  2:968135/1‑149 (MQ=255)
                                                                       gcagcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTg                                                                                                                     >  2:1450275/1‑112 (MQ=255)
                                                                         agcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTAcc                                                                              <  2:1644400/149‑1 (MQ=255)
                                                                              tGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCAtctc                                                                         >  1:1817812/1‑149 (MQ=255)
                                                                              tGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCAtctc                                                                         >  4:1736013/1‑149 (MQ=255)
                                                                              tGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCAtctc                                                                         >  4:1240672/1‑149 (MQ=255)
                                                                                 aaCTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTcc                                                                      >  1:1687913/1‑149 (MQ=255)
                                                                                  aCTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTccc                                                                     <  1:2382492/149‑1 (MQ=255)
                                                                                      cTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGc                                                                                                                                                 >  2:1081768/1‑69 (MQ=255)
                                                                                       ttGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCaa                                                                >  2:185294/1‑149 (MQ=255)
                                                                                       ttGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCaa                                                                <  1:305309/149‑1 (MQ=255)
                                                                                        tGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAAc                                                               >  2:2158148/1‑149 (MQ=255)
                                                                                             cTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTa                                                          <  1:1550479/149‑1 (MQ=255)
                                                                                                ttCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTg                                                       >  1:525597/1‑149 (MQ=255)
                                                                                                ttCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTg                                                       >  3:1089863/1‑149 (MQ=255)
                                                                                                   tGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCg                                                    >  1:118082/1‑149 (MQ=255)
                                                                                                           cTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGgt                                                                                                                   >  2:2522408/1‑78 (MQ=255)
                                                                                                           cTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGtt                                            >  4:1835872/1‑149 (MQ=255)
                                                                                                           cTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGtt                                            <  1:1923089/149‑1 (MQ=255)
                                                                                                           cTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGtt                                            <  3:1776488/149‑1 (MQ=255)
                                                                                                            tGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGttt                                           <  1:899905/149‑1 (MQ=255)
                                                                                                             gAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTa                                          >  2:2015289/1‑149 (MQ=255)
                                                                                                              aaGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTAc                                         >  3:325935/1‑149 (MQ=255)
                                                                                                              aaGTAATTAAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTAc                                         >  2:831977/1‑149 (MQ=255)
                                                                                                                 tAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACata                                      >  2:621480/1‑149 (MQ=255)
                                                                                                                 tAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACata                                      >  4:1446450/1‑149 (MQ=255)
                                                                                                                 tAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACata                                      >  2:378632/1‑149 (MQ=255)
                                                                                                                        gATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGcc                               <  3:105685/149‑1 (MQ=255)
                                                                                                                         aTTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCg                              >  3:956461/1‑149 (MQ=255)
                                                                                                                         aTTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCg                              >  2:1538995/1‑149 (MQ=255)
                                                                                                                              aGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATCCTGCCGCTaaa                         <  3:1446450/149‑1 (MQ=255)
                                                                                                                                   tataaaaGGGGAGCGGTTTCCCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAAtt                    >  1:2344673/5‑149 (MQ=255)
                                                                                                                                    aCCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTc                   >  1:1237330/1‑149 (MQ=255)
                                                                                                                                     taaaaGGGGAGCGGTTTCCCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTcc                  <  2:2344673/147‑1 (MQ=255)
                                                                                                                                     ccAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTcc                  <  3:2145908/149‑1 (MQ=255)
                                                                                                                                     ccAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTcc                  >  1:2300684/1‑149 (MQ=255)
                                                                                                                                      cAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTccc                 >  1:2059743/1‑149 (MQ=255)
                                                                                                                                      cAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTccc                 >  4:609140/1‑149 (MQ=255)
                                                                                                                                       aaaGGGGAGCGGTTTCCCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTcccc                >  3:1591168/1‑149 (MQ=255)
                                                                                                                                        aaGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCa               >  1:138533/1‑149 (MQ=255)
                                                                                                                                           gggAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTaaa                         >  1:76337/1‑136 (MQ=255)
                                                                                                                                              aGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTg         <  4:205789/149‑1 (MQ=255)
                                                                                                                                              aGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTg         <  3:1736013/149‑1 (MQ=255)
                                                                                                                                               gCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTgg        >  2:1699628/1‑149 (MQ=255)
                                                                                                                                               gCGGCAAACCGCTCCCCTTTTATATTTTTCGTGCGTGTTGGCGTCGGATGCGATGTTGACGCATCTTATTCGCCCTATCATCTCTCCCTGCAACACTTATTGCCGCTTGTGTTTACAGATTCTGCCGCTAAACAATTCCCCATCCCTgg        <  2:1096188/149‑1 (MQ=255)
                                                                                                                                                  ggAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGt     <  1:1091854/149‑1 (MQ=255)
                                                                                                                                                     aaCCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGtata  <  4:1563742/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATA  >  NC_000913/3853527‑3853824

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: