Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,141,153 A→G 5.7% intergenic (+167/+29) flgL → / ← rne flagellar hook‑filament junction protein 2/ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,141,1530AG 5.7% 271.7 / 4.9 88intergenic (+167/+29)flgL/rneflagellar hook‑filament junction protein 2/ribonuclease E
Reads supporting (aligned to +/- strand):  ref base A (39/44);  new base G (3/2);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGACG  >  NC_000913/1141005‑1141298
                                                                                                                                                    |                                                                                                                                                 
cATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTAc                                                                                                                                                  >  1:895061/1‑150 (MQ=255)
 aTGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTAcc                                                                                                                                                 >  1:585680/1‑150 (MQ=255)
   gAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAg                                                                                                                                               >  2:787501/1‑150 (MQ=255)
   gAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAg                                                                                                                                               <  1:1283725/150‑1 (MQ=255)
    aaaCTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAgg                                                                                                                                              <  2:734481/150‑1 (MQ=255)
    aaaCTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAgg                                                                                                                                              >  1:205231/1‑150 (MQ=255)
     aaCTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAggg                                                                                                                                             <  1:1809700/150‑1 (MQ=255)
       cTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCt                                                                                                                                           <  1:743559/150‑1 (MQ=255)
        tGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCtt                                                                                                                                          <  1:1996213/150‑1 (MQ=255)
        tGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCtt                                                                                                                                          <  2:895061/150‑1 (MQ=255)
         gggTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTg                                                                                                                                         >  1:981815/1‑150 (MQ=255)
         gggTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTg                                                                                                                                         <  1:386925/150‑1 (MQ=255)
           gTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGAt                                                                                                                                       <  2:1773825/150‑1 (MQ=255)
           gTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGAt                                                                                                                                       <  2:408182/150‑1 (MQ=255)
            tATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCtttttt                                                                                                                                      <  2:354045/150‑5 (MQ=255)
              tGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTAc                                                                                                                                    >  1:1385329/1‑150 (MQ=255)
              tGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTAc                                                                                                                                    >  2:1958017/1‑150 (MQ=255)
               gTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACt                                                                                                                                   >  1:731786/1‑150 (MQ=255)
                  ttGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTg                                                                                                                                <  1:733006/150‑1 (MQ=255)
                   tGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGa                                                                                                                               <  2:1668289/150‑1 (MQ=255)
                   tGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGa                                                                                                                               <  2:1642715/150‑1 (MQ=255)
                      ctgcctgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTGGATTACTTTGAGCt                                                                                                                            >  2:969193/1‑150 (MQ=255)
                      ctgcctgcTCTGGGATCGCTGGGGCGGACATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCt                                                                                                                            <  1:427494/150‑1 (MQ=255)
                       tgcctgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTa                                                                                                                           >  2:476748/1‑150 (MQ=255)
                        gcctgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTaa                                                                                                                          <  2:1385329/150‑1 (MQ=255)
                             ctctGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAattatt                                                                                                                     >  1:7473/1‑150 (MQ=255)
                               ctGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTAc                                                                                                                   >  2:908770/1‑150 (MQ=255)
                                 gggATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTc                                                                                                                 <  2:753783/150‑1 (MQ=255)
                                    aTCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAAc                                                                                                              <  2:625683/150‑1 (MQ=255)
                                       gCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAgg                                                                                                           >  1:1056389/1‑150 (MQ=255)
                                       gCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAgg                                                                                                           >  2:656781/1‑150 (MQ=255)
                                        cgggggCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGt                                                                                                          <  2:1286448/148‑1 (MQ=255)
                                                ggcATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAcg                                                                                                  >  1:1676201/1‑150 (MQ=255)
                                                         ttGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGc                                                                                         <  2:1056389/150‑1 (MQ=255)
                                                          tGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCg                                                                                        <  2:731786/150‑1 (MQ=255)
                                                           gCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCgg                                                                                       >  1:1394196/1‑150 (MQ=255)
                                                             cTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAg                                                                                             <  1:1785661/142‑1 (MQ=255)
                                                             cTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAg                                                                                             >  2:1785661/1‑142 (MQ=255)
                                                              tATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCag                                                                                    <  2:1867155/150‑1 (MQ=255)
                                                                 tttGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAgg                                                                                 <  1:1694324/150‑1 (MQ=255)
                                                                     cATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatg                                                                             <  1:888634/150‑1 (MQ=255)
                                                                       ttgttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatgat                                                                           <  2:205231/150‑1 (MQ=255)
                                                                       ttgttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatgat                                                                           <  1:1958017/150‑1 (MQ=255)
                                                                       ttgttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatgat                                                                           >  2:221629/1‑150 (MQ=255)
                                                                            ggTTAGCAAGGATGCCATTCGACGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGtt                                                                      <  1:392904/150‑1 (MQ=255)
                                                                             gTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATgtgtag                                                                     <  1:476748/150‑3 (MQ=255)
                                                                             gTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTg                                                                     <  1:1020520/150‑1 (MQ=255)
                                                                               tAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGcc                                                                   >  1:878392/1‑150 (MQ=255)
                                                                               tAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTGGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGcc                                                                   <  1:969193/150‑1 (MQ=255)
                                                                                aGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCg                                                                  <  2:981815/150‑1 (MQ=255)
                                                                                  cAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTa                                                                >  1:1004999/1‑150 (MQ=255)
                                                                                  cAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTa                                                                >  2:715347/1‑150 (MQ=255)
                                                                                  cAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTa                                                                <  1:1872534/150‑1 (MQ=255)
                                                                                   aaGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCtt                                                                                                                                          <  2:1494348/75‑1 (MQ=255)
                                                                                   aaGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCtt                                                                                                                                          >  1:1494348/1‑75 (MQ=255)
                                                                                         gCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGaccac                                                         >  2:1210196/1‑150 (MQ=255)
                                                                                          ccATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGaccacc                                                        <  1:117586/150‑1 (MQ=255)
                                                                                          ccATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGaccacc                                                        >  2:476669/1‑150 (MQ=255)
                                                                                              tCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCg                                                    <  1:113391/150‑1 (MQ=255)
                                                                                               cGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCgg                                                   >  2:1492277/1‑150 (MQ=255)
                                                                                                   gAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCAcc                                               <  2:7473/150‑1 (MQ=255)
                                                                                                      ttttAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCttt                                            >  2:1026286/1‑150 (MQ=255)
                                                                                                               gTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCa                                                 >  2:1564903/1‑136 (MQ=255)
                                                                                                               gTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCa                                                 <  1:1564903/136‑1 (MQ=255)
                                                                                                                ttGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAg                                  <  1:1210196/150‑1 (MQ=255)
                                                                                                                  gATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGc                                >  1:1899723/1‑150 (MQ=255)
                                                                                                                   aTTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCa                               <  2:1004999/150‑1 (MQ=255)
                                                                                                                    ttCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCaa                              >  1:782878/1‑150 (MQ=255)
                                                                                                                     tCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCaaa                             <  1:506761/150‑1 (MQ=255)
                                                                                                                      cAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATTACCACCTGCGGCACCTTTACCTTCGAAGGCaaaa                            >  2:97373/1‑150 (MQ=255)
                                                                                                                      cAAAGATGAAATAAAAAAGCCCTGCCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCaaaa                            <  2:878392/150‑1 (MQ=255)
                                                                                                                               aaTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGgcgc                   >  2:1665879/1‑150 (MQ=255)
                                                                                                                                aTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCt                  >  2:224383/1‑150 (MQ=255)
                                                                                                                                aTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCt                  >  2:1604172/1‑150 (MQ=255)
                                                                                                                                  aaaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGc                >  1:1425383/1‑150 (MQ=255)
                                                                                                                                  aaaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGc                <  2:1536858/150‑1 (MQ=255)
                                                                                                                                   atcaaGCCCTGGTAACTGCCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATCACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc               >  2:993201/4‑150 (MQ=255)
                                                                                                                                   aaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc               >  2:418430/1‑150 (MQ=255)
                                                                                                                                   aaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc               <  2:1899723/150‑1 (MQ=255)
                                                                                                                                   aaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc               >  1:1104152/1‑150 (MQ=255)
                                                                                                                                     caaGCCCTGGTAACTGCCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAg             >  2:905526/2‑150 (MQ=255)
                                                                                                                                     aaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAg             <  1:224383/150‑1 (MQ=255)
                                                                                                                                       aGCCCTGGTAACTGCCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTc           >  1:290088/1‑150 (MQ=255)
                                                                                                                                         cccTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCAc         <  1:1746524/150‑1 (MQ=255)
                                                                                                                                         cccTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCAc         <  1:115089/150‑1 (MQ=255)
                                                                                                                                             ggCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGTGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACtgtg     >  2:17750/1‑150 (MQ=255)
                                                                                                                                               cAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGATGCCAGTCACTGTGAc   >  2:409756/1‑150 (MQ=255)
                                                                                                                                                aGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGAcg  >  1:1005189/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGACG  >  NC_000913/1141005‑1141298

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