Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,011 T→A 11.4% intergenic (+171/+45) ybfK → / ← kdpE uncharacterized protein YbfK/DNA‑binding transcriptional activator KdpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,0110TA11.4% 594.1 / 54.3 220intergenic (+171/+45)ybfK/kdpEuncharacterized protein YbfK/DNA‑binding transcriptional activator KdpE
Reads supporting (aligned to +/- strand):  ref base T (57/138);  new base A (5/20);  total (62/158)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 228 total aligned reads displayed.

TCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACG  >  NC_000913/720862‑721157
                                                                                                                                                      |                                                                                                                                                  
        ggAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                                                                            <  2:970153/150‑1 (MQ=255)
        ggAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtt                                                                                                                                            <  1:781057/150‑1 (MQ=255)
         taaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                                                                           <  2:409130/149‑1 (MQ=255)
         gAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                                                                           <  2:892213/150‑1 (MQ=255)
         gAAATACGGTTGAACTGTGAACGGGATCAATAAAAACTGATCAATCTTAATTTATTTAATGATGAGCTTTTAACTCAGTAATCAAAAATATTGAATTGTTATTTTTGTGTGTTGTGTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                                                                           <  2:1399122/150‑1 (MQ=255)
          aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                                                                                                                                          <  2:434190/150‑1 (MQ=255)
          aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                                                                          <  2:47226/150‑2 (MQ=255)
          aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                                                                          <  2:1473654/150‑2 (MQ=255)
          aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTATTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                                                                          <  1:1525951/150‑2 (MQ=255)
           aaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                                                                                                                                         <  2:1092836/150‑1 (MQ=255)
           aaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttt                                                                                                                                         <  1:296062/150‑3 (MQ=255)
            aTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                                                                        <  2:75046/150‑1 (MQ=255)
            aTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                                                                        <  2:825715/150‑1 (MQ=255)
             tACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                                                                       <  2:1008884/150‑1 (MQ=255)
             tACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                                                                       <  1:430078/150‑1 (MQ=255)
             tACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                                                                       <  2:1288645/150‑1 (MQ=255)
               cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                                                                     <  2:1363571/150‑7 (MQ=255)
                ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                                                                    <  2:546748/150‑1 (MQ=255)
                ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                                                                    <  2:10557/150‑8 (MQ=255)
                ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                                                                    <  2:1244340/150‑8 (MQ=255)
                ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                                                                    <  1:349796/150‑8 (MQ=255)
                 gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                                                                   <  2:1188742/150‑1 (MQ=255)
                    gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                                                                                <  2:802130/150‑1 (MQ=255)
                      aCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                                                                              <  2:463854/150‑1 (MQ=255)
                       ctctgAAGGTGATCCATAAAAAATGAACAATCAAAATTTATTTAATGATGAGCTTTTTACTCAGTAAAATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                                                                             <  2:715372/147‑1 (MQ=255)
                       cTGTGAAGGTGATCAATAAAAAATGATCAATCTAAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                                                                             <  2:1229905/150‑1 (MQ=255)
                        tgtgAAGGTTATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                                                                            <  2:808774/150‑1 (MQ=255)
                        tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                                                                            <  2:1218063/150‑1 (MQ=255)
                        tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                                                                            <  2:1163568/150‑1 (MQ=255)
                         gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcat                                                                                                                           >  1:1009882/1‑149 (MQ=255)
                         gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                                                                           >  1:1498404/1‑150 (MQ=255)
                          tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg                                                                                                                          <  2:303693/150‑1 (MQ=255)
                             aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTg                                                                                                                       >  1:874223/1‑150 (MQ=255)
                             aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTg                                                                                                                       <  1:1126062/150‑1 (MQ=255)
                              tgttaTCAATTAAAAATTCTTAATTTTAATTTATTTAATTTATTTGTTTTTAATCATTAAAAAAAAATATTTTATTTTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                                                                      <  2:113724/146‑1 (MQ=255)
                              ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                                                                      >  2:1114008/1‑150 (MQ=255)
                                 gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                                                                   <  2:1180660/150‑1 (MQ=255)
                                  aTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                                                                  >  1:1099453/1‑150 (MQ=255)
                                    cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                                                                       <  2:610011/127‑1 (MQ=255)
                                    cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaa                                                                                                                >  2:1150597/1‑150 (MQ=255)
                                    cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGACGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGATGTATTaaa                                                                                                                >  1:1362130/1‑150 (MQ=255)
                                     aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                                                                               <  2:621723/150‑1 (MQ=255)
                                     aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                                                                               <  2:8978/150‑1 (MQ=255)
                                      aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                                                                              >  2:833208/1‑150 (MQ=255)
                                      aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                                                                              >  1:705968/1‑150 (MQ=255)
                                       tAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                                                                             <  2:196386/150‑1 (MQ=255)
                                       tAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                                                                             <  2:527043/150‑1 (MQ=255)
                                        aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                                                                            >  1:1286582/1‑150 (MQ=255)
                                         aaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                                                                           >  2:403289/1‑150 (MQ=255)
                                           aaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTATTTATTATTACGTCTTTAATTAATAAAATGAAGGCTTTATTAAAATTAata                                                                                                         <  2:100286/150‑1 (MQ=255)
                                           aaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                                                                         <  2:960921/150‑1 (MQ=255)
                                               taTCAATATTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTGATAAAAGGCAGGCTGTATTAAAATTAATATTCa                                                                                                     <  2:91020/149‑1 (MQ=255)
                                                  cAATCTTAATTTATTTAATGATGACCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                                                                  <  1:330122/150‑1 (MQ=255)
                                                   aaTCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                                                                 <  2:1211708/150‑1 (MQ=255)
                                                    aTCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                                                                                <  1:1447673/150‑1 (MQ=255)
                                                          atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCATGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                                                                          <  2:222754/150‑1 (MQ=255)
                                                          atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                                                                          <  1:61811/150‑1 (MQ=255)
                                                            ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                                                                        <  2:976965/150‑1 (MQ=255)
                                                            ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                                                                        <  2:1099453/150‑1 (MQ=255)
                                                              atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                      >  1:1077633/1‑150 (MQ=255)
                                                              atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                      <  1:163223/150‑1 (MQ=255)
                                                                ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                                                                          <  2:1123336/144‑1 (MQ=255)
                                                                ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGc                                                                                    >  1:537806/1‑150 (MQ=255)
                                                                 taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                                                                   >  2:218275/1‑150 (MQ=255)
                                                                   atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGATTTAATTAATAAAAGGCATGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGACaa                                                                                 >  1:1021420/1‑150 (MQ=255)
                                                                      gatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                                                                               <  1:403289/150‑1 (MQ=255)
                                                                          aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                                                                           >  2:1411536/1‑150 (MQ=255)
                                                                              ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                                                                       <  2:656107/150‑1 (MQ=255)
                                                                                ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                                                                     <  1:84185/150‑1 (MQ=255)
                                                                                   cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACAGATAGCCAATACCGGTTTCAGTaa                                                                  >  1:1434549/1‑150 (MQ=255)
                                                                                     cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                                                                  >  2:397524/1‑148 (MQ=255)
                                                                                     cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                                                                >  2:251791/1‑150 (MQ=255)
                                                                                     cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                                                                >  2:646743/1‑150 (MQ=255)
                                                                                             ttaaaTATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg                                                        <  1:111520/148‑1 (MQ=255)
                                                                                             taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg                                                        <  1:829198/150‑1 (MQ=255)
                                                                                               aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTg                                                      <  2:537806/150‑1 (MQ=255)
                                                                                                     tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggcggg                                                <  2:1498404/150‑1 (MQ=255)
                                                                                                        aTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGgcgggcgg                                             <  1:246246/150‑1 (MQ=255)
                                                                                                             tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATcc                                        >  2:101749/1‑150 (MQ=255)
                                                                                                               tttttGTGTGTTGTTTAACATAAAAAGCCGTATTTATTTTTACGGCTTTAAATAATAAAAGGCAGGCTGTTTTAAAATTAATTTTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTg                                      <  2:1487514/150‑1 (MQ=38)
                                                                                                               tttttGTGTGTTGTGTAACATAAAAAGCCGTATTTATTATTACGGCTTTAATAAATAAAACGCAGGCGGTATTAAAATTAAAATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGGGGGCGGGCGGGATCCTg                                      <  2:1513104/150‑1 (MQ=37)
                                                                                                                 tttGTGTTTTGTTTAACATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGtt                                    <  2:1445494/150‑1 (MQ=255)
                                                                                                                         ttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGtttt                            <  1:72468/150‑1 (MQ=255)
                                                                                                                          tgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGttttt                           <  2:730391/150‑1 (MQ=255)
                                                                                                                           gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTg                          <  2:345047/150‑1 (MQ=255)
                                                                                                                           gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTg                          <  1:218275/150‑1 (MQ=255)
                                                                                                                           gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTg                          <  1:518276/150‑1 (MQ=255)
                                                                                                                             ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGttt                             <  1:1234094/145‑1 (MQ=255)
                                                                                                                               aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGc                      >  2:504826/1‑150 (MQ=255)
                                                                                                                               aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGc                      <  1:940365/150‑1 (MQ=255)
                                                                                                                               aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCGGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGc                      <  1:1026103/150‑1 (MQ=255)
                                                                                                                                aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCa                     <  1:59969/150‑1 (MQ=255)
                                                                                                                                 gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAg                    <  1:1411536/150‑1 (MQ=255)
                                                                                                                                   attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGAt                  >  2:875058/4‑150 (MQ=25)
                                                                                                                                     taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAAGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGt                >  1:181469/2‑150 (MQ=25)
                                                                                                                                         gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCa            <  1:1104095/150‑1 (MQ=255)
                                                                                                                                          ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCat           <  1:1471579/150‑1 (MQ=255)
                                                                                                                                             tatttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCAtata        >  2:1064188/1‑150 (MQ=255)
                                                                                                                                                 tattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCTATAGCCAATACCGGTTTCAGTAATGAAATGGCTTGGGAGGGCGGGATCCTGTTCAAGTTTTTGTCGCAGATGTCCCATATAAATa    >  1:385608/1‑150 (MQ=255)
                                                                                                                                                  attattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAACGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATAc   >  2:259704/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACG  >  NC_000913/720862‑721157

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