Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,599 T→A 6.6% intergenic (+32/+387) hdeD → / ← arrS acid‑resistance membrane protein/antisense sRNA ArrS, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,657,5990TA 6.6% 533.9 / 21.2 152intergenic (+32/+387)hdeD/arrSacid‑resistance membrane protein/antisense sRNA ArrS, function unknown
Reads supporting (aligned to +/- strand):  ref base T (65/77);  new base A (8/2);  total (73/79)
Fisher's exact test for biased strand distribution p-value = 4.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

Only 100 of 158 total aligned reads displayed.

TTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTA  >  NC_000913/3657333‑3657855
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                
       tgtgTTGCAGGGGTTATTCGCTTGATTAGTTGGTATCGTCAGCGTTCAAGGAAAGGCAGCTGGCTACAGCTCGTTATTGGGGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                   <  1:23450/275‑1 (MQ=255)
       ggtgtTGCAGGGGTAATTCGCCTGATGAGGTGGTATCGTAAGCGTTAAATAAAAGGCAGCTGGCTACAGCTCGTTATTGGTGGGCTGGATATTGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACGGGTGGGAATTGAACTGAAATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                   <  2:539498/274‑1 (MQ=255)
         tgtTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTACAACGCCGATGGGGTATGTAACGCTGGTTTCCACACTGGTGGGCATTGCACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAACACGCAATAAGTGCCCCGGTTTTCACCCGGACCATAGTTACTTGATGTTACTAttact                                                                                                                                                                                                                                                  >  4:546949/1‑270 (MQ=255)
                                tgaGTGGGTCTCGTCCGCGTTAAAGGAAAGCCAGCTGGGTACAGTCAGTTGTTGGTGTGTTGGATTTCGCCATGTCTCGGATCTTCCTTGGTGACACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGCTATTTAGCGCCGCCAGCCTGTTCGGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGTATCATAGTCGCTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACg                                                                                                                                                                                                                           <  4:510083/274‑1 (MQ=255)
                                         aTCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGACTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGACCCTAGTCCCTTGCTGTGACTATGATCTGATTCAAACTCTTTCCCCTa                                                                                                                                                                                                                             >  2:42072/1‑263 (MQ=255)
                                            ggtaGCGTTTAATGAAAGGATGCTGTGTACAGTTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTTTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                                                                                                                                                                              <  2:81864/272‑1 (MQ=255)
                                              cAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATCCTCCTCCCGCTACGCAGTGTTGTAgtaca                                                                                                                                                                                                            >  4:318599/1‑271 (MQ=255)
                                                       aTGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAATCTCTCCGCTACGCAGTGCTGTAGATCAATTGCGCAc                                                                                                                                                                                                   >  2:324813/1‑275 (MQ=255)
                                                                aGCTGGTTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                                                                                                                                                                           <  2:320430/274‑1 (MQ=255)
                                                                     gCTCCAGCTGGTTATTGGTGTGCTGGTTATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                                                                                                                                                                                     <  2:105810/275‑1 (MQ=255)
                                                                                 tATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTa                                                                                                                                                                                                                             >  2:496496/1‑223 (MQ=255)
                                                                                         tgCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTAttt                                                                                                                                                                 >  2:576160/1‑275 (MQ=255)
                                                                                          gCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATCCTCGCCCCGCTACCAAGTGTTGTAGCTCAATTGCGCACTATCATTGAAATAATTTCCCGCCAGTGATTAttt                                                                                                                                                                 >  4:138281/1‑274 (MQ=255)
                                                                                                            cTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCCCCCTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                                                                                                                                                                           <  2:8105/230‑1 (MQ=255)
                                                                                                            cTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCa                                                                                                                                                               <  2:518659/258‑1 (MQ=255)
                                                                                                             tGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATc                                                                                                                                              <  1:102144/274‑1 (MQ=255)
                                                                                                                           tGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTTGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCAct                                                                                                                                >  4:494447/1‑274 (MQ=255)
                                                                                                                             cAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                                                                                                                                                                                            >  6:800657/1‑148 (MQ=255)
                                                                                                                               aCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                                                                                                                                                                          >  5:84924/1‑148 (MQ=255)
                                                                                                                               aCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                                                                                                                                                                          <  6:767954/148‑1 (MQ=255)
                                                                                                                                cGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTa                                                                                                                                                                                                                             >  4:628976/1‑176 (MQ=255)
                                                                                                                                cGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCAGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                                                                                                                                                                                     >  4:450703/1‑216 (MQ=255)
                                                                                                                                 gCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACt                                                                                                                                                                                                                                                        >  6:388460/1‑148 (MQ=255)
                                                                                                                                 gCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACg                                                                                                                                                                                                                           <  2:742836/177‑1 (MQ=255)
                                                                                                                                  ccGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTa                                                                                                                                                                                                                                                       <  6:285286/148‑1 (MQ=255)
                                                                                                                                  ccGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                                                                                                                                                                                      <  7:512506/149‑1 (MQ=255)
                                                                                                                                    gATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATg                                                                                                                                                                                                                                                     <  5:535256/148‑1 (MQ=255)
                                                                                                                                    gATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATg                                                                                                                                                                                                                                                     <  5:879964/148‑1 (MQ=255)
                                                                                                                                           tCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                                                                                                                                                                              <  5:200555/148‑1 (MQ=255)
                                                                                                                                            cTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                                                                                                                                                                             >  5:739228/1‑148 (MQ=255)
                                                                                                                                             tGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCGTCTAATTTTTTTCACTCTATGGCAAATTAGCCAtt                                                                                                             >  4:609387/1‑275 (MQ=255)
                                                                                                                                                 aCGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAc                                                                                                                                                                                                                                        <  5:388460/148‑1 (MQ=255)
                                                                                                                                                 aCGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCGTCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCa                                                                                                          <  2:609389/274‑1 (MQ=255)
                                                                                                                                                  cGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACtt                                                                                                                                                                                                                                                                 <  6:104913/122‑1 (MQ=255)
                                                                                                                                                   gCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCACGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                     <  7:655503/149‑1 (MQ=255)
                                                                                                                                                      ggTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctct                                                                                                                                                                                                                                   <  5:140847/148‑1 (MQ=255)
                                                                                                                                                      ggTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctct                                                                                                                                                                                                                                   >  6:645324/1‑148 (MQ=255)
                                                                                                                                                             acacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATgg                                                                                             <  4:576158/275‑1 (MQ=255)
                                                                                                                                                                  ggtggGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtg                                                                                                                                                                                                                      <  7:743472/149‑1 (MQ=255)
                                                                                                                                                                    tggGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                                                                                                                                                                                     <  5:800657/148‑1 (MQ=255)
                                                                                                                                                                       gAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                                                                                                                                                 <  8:671905/149‑1 (MQ=255)
                                                                                                                                                                        aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAg                                                                                                                                                                                                                >  8:525354/1‑149 (MQ=255)
                                                                                                                                                                         aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAg                                                                                                                                                                                                                <  6:84924/148‑1 (MQ=255)
                                                                                                                                                                           tGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                                                                                                                                                                         >  4:520345/1‑137 (MQ=255)
                                                                                                                                                                             aaCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCAt                                                                                                              >  2:13822/1‑242 (MQ=255)
                                                                                                                                                                              aCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACtt                                                                                                                                                                                                                                                                 >  2:465266/1‑94 (MQ=255)
                                                                                                                                                                              aCTGATATTTAGCGCCGCCAGCCTATTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATAATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCaa                                                                                                                                                                                                           >  5:285929/1‑148 (MQ=255)
                                                                                                                                                                               cTGATATTTAGCGCCGCCCGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                     <  2:64465/121‑1 (MQ=255)
                                                                                                                                                                               cTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                                                                                                                                                                                          <  5:645324/148‑1 (MQ=255)
                                                                                                                                                                               cTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                                                                                                                                                                                          <  6:541944/148‑1 (MQ=255)
                                                                                                                                                                                tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                                                                                                                                                                                         >  5:340594/1‑148 (MQ=255)
                                                                                                                                                                                 gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCg                                                                                                                                                                                                                                <  4:378560/124‑1 (MQ=255)
                                                                                                                                                                                 gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                                                                                                                                                                                         <  2:294239/163‑1 (MQ=255)
                                                                                                                                                                                 gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCa                                                                                                          <  4:13822/242‑1 (MQ=255)
                                                                                                                                                                                  atatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                                                                                                                                                                             <  4:465264/94‑1 (MQ=255)
                                                                                                                                                                                  atatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                                                                                                                                                                         >  4:817255/1‑130 (MQ=255)
                                                                                                                                                                                  atatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcg                                                                                                                                                                                                      >  4:517214/1‑149 (MQ=255)
                                                                                                                                                                                    atTTAGCGCCGCCAGCCTATTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                                                                                                                                                                                     <  6:285929/148‑1 (MQ=255)
                                                                                                                                                                                      ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACt                                                                                                                                                                                                  <  2:517216/149‑1 (MQ=255)
                                                                                                                                                                                            ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATa                                                                                                                                                                                                                                         >  4:760906/1‑104 (MQ=255)
                                                                                                                                                                                               ccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGc                                                                                                                                                                                                                               >  4:726629/1‑111 (MQ=255)
                                                                                                                                                                                                cAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                     <  2:760908/104‑1 (MQ=255)
                                                                                                                                                                                                    cTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataat                                                                                                                                                                                    <  8:305052/149‑1 (MQ=255)
                                                                                                                                                                                                       ttCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                                                                                                                                                                                 >  7:111313/1‑149 (MQ=255)
                                                                                                                                                                                                        tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                                                                                                                                                                                 >  5:398698/1‑148 (MQ=255)
                                                                                                                                                                                                           gCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGc                                                                                                                                                                             <  8:283652/149‑1 (MQ=255)
                                                                                                                                                                                                            cTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgt                                                                                                                                                                                                                       >  6:354782/1‑106 (MQ=255)
                                                                                                                                                                                                             ttCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt                                                                                                                                                                            <  6:739228/148‑1 (MQ=255)
                                                                                                                                                                                                             ttCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                                                                                                                                                                           <  8:822004/149‑1 (MQ=255)
                                                                                                                                                                                                             ttCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGaa                                             <  2:453241/275‑1 (MQ=255)
                                                                                                                                                                                                               cGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTa                                                                                                                                                                                                                             >  2:81427/1‑97 (MQ=255)
                                                                                                                                                                                                                gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                                                                                                                                                 >  7:827725/1‑108 (MQ=255)
                                                                                                                                                                                                                 ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACATGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTg                                         >  2:434518/1‑275 (MQ=255)
                                                                                                                                                                                                                   aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                                                                                                                                                                         <  4:81426/97‑1 (MQ=255)
                                                                                                                                                                                                                        gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTc                                                                                                                                                                >  7:235056/1‑149 (MQ=255)
                                                                                                                                                                                                                         ttCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCa                                                                                                                                                               >  7:451184/1‑149 (MQ=255)
                                                                                                                                                                                                                            gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAcc                                                                                                                                                            >  8:709580/1‑149 (MQ=255)
                                                                                                                                                                                                                             ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAcc                                                                                                                                                            >  5:241300/1‑148 (MQ=255)
                                                                                                                                                                                                                                           ttACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAg                                                                                                                                                                                                                >  6:323438/1‑82 (MQ=255)
                                                                                                                                                                                                                                              ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctct                                                                                                                                                                                                                                   >  5:377095/1‑60 (MQ=255)
                                                                                                                                                                                                                                              ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAt                                                                                                                                          >  7:145739/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   gTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                                                                                                                                                 >  7:545912/1‑73 (MQ=255)
                                                                                                                                                                                                                                                   gTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                                                                                                                                                 <  8:545912/73‑1 (MQ=255)
                                                                                                                                                                                                                                                   gTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTTAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGc       >  2:234597/1‑275 (MQ=255)
                                                                                                                                                                                                                                                    ttGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                                                                                                                                                                                  >  2:200210/1‑71 (MQ=255)
                                                                                                                                                                                                                                                    ttGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAttttttt                                                                                                                                    <  8:111313/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        cACCCGGTTCATAGTACCTTGATGTGACTATGTACCGATTAAAACTTTCTCCGCTTCGCAGTGTTGTAGTTAAATTGCGCACTATCATTGAAATAATTACCTGTTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGCTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTAGGGGCGGTTa  <  4:234596/275‑1 (MQ=255)
                                                                                                                                                                                                                                                        cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                                                                                                                                                                              <  4:200209/71‑1 (MQ=255)
                                                                                                                                                                                                                                                             ggATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTAt                                                                                                                            <  6:398698/148‑1 (MQ=255)
                                                                                                                                                                                                                                                              gATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTa                                                                                                                                                                                                                             >  2:718255/1‑50 (MQ=255)
                                                                                                                                                                                                                                                                tCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTAt                                                                                                                            <  7:425514/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                tCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTAt                                                                                                                            >  8:425514/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                tCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGc                                                                                                                         >  5:701878/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 cATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCa                                                                                                                        >  6:528667/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                  aTAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                                                                                                                                                                         <  4:718253/50‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAAt                                                                                                                     >  6:369446/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     gTCACATCAAGTGACTATGATCCGTTTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAAtt                                                                                                                    >  6:123607/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                      tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTa                                                                                                                   >  5:147741/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                        aCTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGcc                                                                                                                >  7:44557/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGcc                                                                                                                >  5:336109/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                
TTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTA  >  NC_000913/3657333‑3657855

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