Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,009 T→A 6.8% intergenic (+169/+47) ybfK → / ← kdpE uncharacterized protein/response regulator in two‑component regulatory system with KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,0090TA 6.8% 306.6 / 7.7 114intergenic (+169/+47)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (46/60);  new base A (3/5);  total (49/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

ATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACG  >  NC_000913/720860‑721157
                                                                                                                                                     |                                                                                                                                                    
atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                                                                      >  1:848483/1‑150 (MQ=255)
 ttCTATCCATGAAATACGGTTGAACTGTTAATGTTATCAATAATAAATGATCAATATTATTTTATTTAATGATGATCTTTTTACTCAGTAATATAATATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCAGTATttatta                                                                                                                                                     <  2:176873/150‑1 (MQ=255)
 ttCTATCCAGGAAATACGGTTTAACTGTGAAGGTGATCAATAAAAAAAGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGATATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                                                                                                                                     <  2:1054308/150‑1 (MQ=255)
   cTAACCAGGAAATACGGTTGAACTGTGAAGGTGATCAAAAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatt                                                                                                                                                   <  2:1160843/150‑1 (MQ=255)
    tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                                                                                                                                                  <  1:1467976/150‑1 (MQ=255)
       ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGATATTTTTGTGTGTTGTTTAAGATAAAAAGACGTATttattattaaag                                                                                                                                               >  1:80719/1‑147 (MQ=255)
         atgAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                                                                             <  1:18269/148‑1 (MQ=255)
         aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                                                                             >  2:1481055/1‑150 (MQ=255)
         aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                                                                             <  2:222077/150‑1 (MQ=255)
          ggAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCACTCTTAATTTATTTAATGATGAGCTGTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTATTTATTATTACGGCtt                                                                                                                                            >  2:267692/1‑150 (MQ=255)
          ggAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                                                                            >  2:763879/1‑150 (MQ=255)
           gAAATTCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                                                                           <  2:170496/150‑1 (MQ=255)
           gAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                                                                           <  2:461154/150‑1 (MQ=255)
           gAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                                                                           <  2:223279/150‑1 (MQ=255)
            ataTACTGTTGAACTGTGAAGGTGTTCAATTAATAATGATCAATCTTAATTTATTTAATGATGATGTTTTTAATCAGTAATTTAAAATTTTGTATTTTTATTTTTGTGTGTTTTTTAAGTTAAAAAGCCGTTATTATAAATACGGCtttt                                                                                                                                          <  2:568862/148‑2 (MQ=255)
            aaaTACGGTTGAACTGTGAAGGTTATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATCCGGCtttt                                                                                                                                          <  1:1530660/150‑2 (MQ=255)
            aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                                                                          <  2:1451331/150‑2 (MQ=255)
            aaaTACGGTTGAACTGAGAAGGAGATCAATAAAAAATGATCAATCTTAATTTATTTTATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAACGCCGTAATAATAACTACGGCtttt                                                                                                                                          <  1:95654/150‑2 (MQ=255)
               tACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                                                                       <  2:1601294/150‑1 (MQ=255)
                 cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                                                                                     >  2:1029768/1‑150 (MQ=255)
                 cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATAAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                                                                     <  1:611203/150‑7 (MQ=255)
                  ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                                                                    <  1:1377019/150‑1 (MQ=255)
                   gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                                                                   >  2:1709742/1‑150 (MQ=255)
                   gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGACGTATTTATTATTACGGCTTTAATtaata                                                                                                                                   >  1:1038018/1‑150 (MQ=255)
                      taaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                                                                                <  2:1706702/149‑1 (MQ=255)
                      gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                                                                                >  2:242984/1‑150 (MQ=255)
                      gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                                                                                >  2:631580/1‑150 (MQ=255)
                      gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATAGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                                                                                >  1:181595/1‑150 (MQ=255)
                          tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                                                                            <  2:773176/150‑1 (MQ=255)
                           gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                                                                           >  1:1032796/1‑150 (MQ=255)
                           gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGACGTATTTATTATTACGGATTTAATTAATAAAaggcag                                                                                                                           >  1:978335/1‑150 (MQ=255)
                            tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg                                                                                                                          <  1:1581042/150‑1 (MQ=255)
                               aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGATTTGTTATTTTTGTGAGTTGTGTTAGATAAAACGCCGTATTTATTATTACGGCTTTAAGTATTAAAAGGCAGTCtt                                                                                                                       >  2:300862/1‑149 (MQ=255)
                                ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGACGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                                                                      >  1:911991/1‑150 (MQ=255)
                                  tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                                                                                    >  2:1386030/1‑150 (MQ=255)
                                   taTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                                                                   <  1:635105/149‑1 (MQ=255)
                                      cAATAAAAAATGGTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaa                                                                                                                >  2:1587952/1‑150 (MQ=255)
                                       aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                                                                  <  1:61682/147‑1 (MQ=255)
                                       aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                                                                  >  2:61682/1‑147 (MQ=255)
                                          aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                                                                            <  2:1148925/150‑1 (MQ=255)
                                          aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                                                                            <  1:1531061/150‑1 (MQ=255)
                                           aaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCATTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                                                                           <  1:108263/150‑1 (MQ=255)
                                                 gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATTAAAGGAAGGATTTATTAAAATTAATATTCa                                                                                                     >  1:285909/1‑150 (MQ=255)
                                                  aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                                                                                    <  2:1305327/150‑1 (MQ=255)
                                                    cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                                                                  <  1:1139471/150‑1 (MQ=255)
                                                        cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                                                                              >  1:839632/1‑150 (MQ=255)
                                                           aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                           <  2:1612688/150‑1 (MQ=255)
                                                            atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                                                                          <  1:1723766/150‑1 (MQ=255)
                                                            atttatttAATGATGAGCTTTTTACTAAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGGTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                                                                          <  1:1301137/150‑1 (MQ=255)
                                                                atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                      >  2:841692/1‑150 (MQ=255)
                                                                atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                      <  2:1418641/150‑1 (MQ=255)
                                                                atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                      >  2:891587/1‑150 (MQ=255)
                                                                   taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                                                                   >  1:386533/1‑150 (MQ=255)
                                                                   taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACAGATAGcc                                                                                   >  1:1261040/1‑150 (MQ=255)
                                                                    aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                                                                                  >  2:113975/1‑150 (MQ=255)
                                                                        atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                                                                              <  2:539079/150‑1 (MQ=255)
                                                                        atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                                                                              <  2:545953/150‑1 (MQ=255)
                                                                         tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAcc                                                                             <  1:137895/150‑1 (MQ=255)
                                                                                 ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                                                                     <  1:358724/150‑1 (MQ=255)
                                                                                 ttACTCAGAAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                                                                     <  1:309869/150‑1 (MQ=255)
                                                                                  tACTCATTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                                                                    <  2:1063863/150‑1 (MQ=255)
                                                                                  tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                                                                    <  1:481696/150‑1 (MQ=255)
                                                                                  tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                                                                    >  1:576813/1‑150 (MQ=255)
                                                                                    cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                                                                  >  2:1296760/1‑150 (MQ=255)
                                                                                     tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                                                                 <  1:789106/150‑1 (MQ=255)
                                                                                      cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                                                                <  2:1442143/150‑1 (MQ=255)
                                                                                      cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                                                                <  2:839632/150‑1 (MQ=255)
                                                                                      cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                                                                >  2:1498418/1‑150 (MQ=255)
                                                                                       aGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAATTAAGAATTTATTAAAATTAATATTCAAAGAATAAAACTATAGACAATACCGGTTTCAGTAATGa                                                               >  1:1228866/1‑150 (MQ=255)
                                                                                        ttAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAATGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa                                                              <  1:1710193/149‑1 (MQ=255)
                                                                                          aaTATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt                                                            >  2:381876/1‑150 (MQ=255)
                                                                                                   tatTGAATTGTTATTTTTGTGTGTTGTTTAAGATTAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggc                                                   <  2:513726/150‑1 (MQ=255)
                                                                                                      tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggcggg                                                >  2:6410/1‑150 (MQ=255)
                                                                                                       gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggcgggc                                               >  2:1319998/1‑150 (MQ=255)
                                                                                                       gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggcgggc                                               <  2:848483/150‑1 (MQ=255)
                                                                                                        aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGggcgggcg                                              <  1:1587952/150‑1 (MQ=255)
                                                                                                          ttGTTTTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGcgggcggg                                            <  2:80719/150‑1 (MQ=255)
                                                                                                          ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGcgggcggg                                            <  2:591492/150‑1 (MQ=255)
                                                                                                          ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGcgggcggg                                            <  2:1032796/150‑1 (MQ=255)
                                                                                                           tGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                                                                          >  1:1255541/1‑87 (MQ=255)
                                                                                                           tGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                                                                          <  2:1255541/87‑1 (MQ=255)
                                                                                                            gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGAt                                          <  1:423479/150‑1 (MQ=255)
                                                                                                              tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCTGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATcc                                        <  1:936777/150‑1 (MQ=255)
                                                                                                              tATTTTTGGGGGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAGTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTGCAGTAATAAAATGGCGTGGGCGGGCGGGATcc                                        <  1:626898/150‑1 (MQ=255)
                                                                                                                    tgtgtgTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTcc                                  >  2:1398689/1‑150 (MQ=255)
                                                                                                                            gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTg                          <  1:1697202/150‑1 (MQ=255)
                                                                                                                              ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTc                        <  2:630753/150‑1 (MQ=255)
                                                                                                                              ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTc                        <  2:1354191/150‑1 (MQ=255)
                                                                                                                                 aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCa                     >  1:1439175/1‑150 (MQ=255)
                                                                                                                                 aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCAAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATACTTTTCAAGTTTTTGTCGCa                     >  1:847543/1‑150 (MQ=255)
                                                                                                                                     aaaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATg                 >  2:776897/1‑150 (MQ=255)
                                                                                                                                       aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTc               >  2:1414331/1‑150 (MQ=255)
                                                                                                                                           ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCTGGCTGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCat           <  1:775591/150‑1 (MQ=255)
                                                                                                                                           ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCat           <  1:1709742/150‑1 (MQ=255)
                                                                                                                                            cGTATTTATTATTACTGATTTAATTAATAAATGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAACGATAGCCAATACGGGTTTAAGTAATGAACTGGGGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCata          <  2:1228866/150‑1 (MQ=25)
                                                                                                                                             gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCatat         <  1:891587/150‑1 (MQ=255)
                                                                                                                                              ttttttttATTACGGGTTTAATTAAAAAAAGGCATGCTTTATTAAAATTAAAATTCAAACCAAAAACCGAAACCCAAAACCGGTTTCATTAATGAAATGGGGGGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCAtata        <  2:386533/148‑1 (MQ=18)
                                                                                                                                                  tattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTATTGTCGCAGATGTCCCATATAAATa    >  2:162147/1‑150 (MQ=255)
                                                                                                                                                    ttattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACg  <  1:1273625/150‑1 (MQ=255)
                                                                                                                                                    ttattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGAGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACg  >  1:954458/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACG  >  NC_000913/720860‑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: