Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,601 G→C 6.6% intergenic (+34/+385) 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,6010GC 6.6% 897.3 / 43.2 259intergenic (+34/+385)hdeD/arrSacid‑resistance membrane protein/antisense sRNA ArrS, function unknown
Reads supporting (aligned to +/- strand):  ref base G (109/132);  new base C (9/8);  total (119/140)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

Only 100 of 261 total aligned reads displayed.

TTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTAA  >  NC_000913/3657333‑3657856
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                               
tttttttgtggtggCAGGGGTTATGGTCCTGNNGAGTTGGGATTGTATGCTTGTAATGGAATGNAACTGGGTACCGGTCGTTTTTGTTGGGCTGGATATGGTCATGGCCGGGATATTCCGTGGTGAAACGGAGATTGTNTCTGGAACGCTGGTTTCAACACTGGTGGGAATTGAACTGATATTGAGCGCCGCCAGCCTGTTCAGCTTGGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACt                                                                                                                                                                                                                                                         <  3:107779/264‑1 (MQ=255)
         tgtTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                                                                                                                                                                                   >  2:671903/1‑274 (MQ=255)
              cAGGGGTTATTCGCCTGATGAGTTGGTACCGTCAGCGTTCAATGAAAGGCAGCTGTCTACAGCTCGTTATTGGTGTTCTGCATATCGTCATCGCCTGACTATTTCTTGGTGCAATTCCGATTGTGTCTCTAACGCTGTTTTCCTCTCTGGTGGGAATTGAACTGATATTAAGCGCCGCCAGCCTTGTCAGCTTCGCCAGTTTGTTTCTTAAGAAGCAAACATTACCCCGGTTTGCAACCGGCGCATATTCACTTGAAGTGACTATGATACGagt                                                                                                                                                                                                                                              >  4:620531/1‑272 (MQ=255)
              acggggttttctcccTGGTGGGTTGGTTTGGTCATGGTGCAATGAAAGGCAGAGGGTTACAGCTTGTTAATGGTGTGTTGGGTATGGTCATTGCCTGTATATTCCTTGGTGCAACGCCGTTGGTGTCTGTAACGCTGGTTTCCACACTGCTGGTAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTa                                                                                                                                                                                                                                             <  4:404413/262‑1 (MQ=255)
                      aTTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTTTGCTGGATAACGTCATTGCCTGGATATTCCTTGGTTCATCGCCTATGGTGTCTGTAACGCCGGTTTTCACACTGGTGGGAATTGAACCGATATTTAGCGCCGCCAGCCTGTTCAGCGTCGCCAGCCTGTTCGCTAAGCATCAATCCTTCACCCGGGTCTCCCCCGGAACATCGTCCCTCGAGGTGACTATGTTCCCACTAAtccact                                                                                                                                                                                                                                      >  4:154288/1‑269 (MQ=255)
                              gatgaGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                                                                                                                                                                          <  1:157997/278‑1 (MQ=255)
                                tgaGTTGGTATCGTCGGCGGTCAATGAAAGGCAGCTGGTTACAGCTCGTGATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGAGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACg                                                                                                                                                                                                                            <  2:191663/274‑1 (MQ=255)
                                 gaGTTGGTATCGTAAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGc                                                                                                                                                                                                                           <  2:140237/274‑1 (MQ=255)
                                     gggTATCGTCAGCGTTCAATGAAAGGCAGCTGGTTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCCGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                                                                                                                                                                                      <  2:590024/274‑1 (MQ=255)
                                                                  cTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGa                                                                                                                                                                                                                >  2:508183/1‑252 (MQ=255)
                                                                                           cTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTTAAATAATTACCTTCTAGTGATTATTTCAACCt                                                                                                                                                            >  1:80117/1‑279 (MQ=255)
                                                                                                                         ggTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATATGACTATGATCCGATTAATACTCTCTCCGCTACg                                                                                                                                                                                                                            >  1:252309/1‑185 (MQ=255)
                                                                                                                            gCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                                                                                                                                                                                             <  6:707623/149‑1 (MQ=255)
                                                                                                                            gCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                                                                                                                                                                                             <  5:1008558/149‑1 (MQ=255)
                                                                                                                              aaCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                                                                                                                                                                           <  6:619939/149‑1 (MQ=255)
                                                                                                                               aCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGACCATAGTCACTTGATGGGAc                                                                                                                                                                                                                                                          >  6:1442528/1‑149 (MQ=255)
                                                                                                                                 gCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                                                                                                                                        <  5:625881/149‑1 (MQ=255)
                                                                                                                                  ccGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                       <  6:468202/149‑1 (MQ=255)
                                                                                                                                  ccGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                       <  6:1299059/149‑1 (MQ=255)
                                                                                                                                  ccGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACCGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                                                                                                                                                                                       <  5:6963/149‑1 (MQ=255)
                                                                                                                                   cGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                       <  5:904996/148‑1 (MQ=255)
                                                                                                                                    gATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGa                                                                                                                                                                                                                                                     <  5:1346806/149‑1 (MQ=255)
                                                                                                                                     aTGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                    <  7:1068477/149‑1 (MQ=255)
                                                                                                                                     aTGGTGTCCGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATACTTACCCCGGTTGTCACCCGGATCATAGTCACTTGATTTGACTATGATCCGATTAATACTCTCTCCGCTACGAAGTGTTGTAGCTCAATTgcgcc                                                                                                                                                                                                     >  1:235134/1‑195 (MQ=255)
                                                                                                                                       ggTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATcc                                                                                                                                                                                                                                                  <  5:1366275/149‑1 (MQ=255)
                                                                                                                                         tgtCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGTGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                                                                                                                                                                                >  6:169426/1‑149 (MQ=255)
                                                                                                                                             tGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTaa                                                                                                                                                                                                                                            >  8:1159234/1‑149 (MQ=255)
                                                                                                                                               tAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATa                                                                                                                                                                                                                                          >  5:1524413/1‑149 (MQ=255)
                                                                                                                                                aaCGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAc                                                                                                                                                                                                                                         >  6:899781/1‑149 (MQ=255)
                                                                                                                                                  cGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActc                                                                                                                                                                                                                                       <  8:1313231/149‑1 (MQ=255)
                                                                                                                                                    cTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctc                                                                                                                                                                                                                                     <  8:1490755/149‑1 (MQ=255)
                                                                                                                                                    cTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctc                                                                                                                                                                                                                                     <  6:410924/149‑1 (MQ=255)
                                                                                                                                                    cTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctc                                                                                                                                                                                                                                     >  5:410924/1‑149 (MQ=255)
                                                                                                                                                     tGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATcc                                                                                                                                                                                                                                                  >  5:1122290/1‑135 (MQ=255)
                                                                                                                                                      ggTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACtctctc                                                                                                                                                                                                                                   <  5:838370/149‑1 (MQ=255)
                                                                                                                                                       gTTTCCNCACTGGTGNNAATTGAACTGATATTTAGCGCCGNCATCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                      <  3:63695/145‑1 (MQ=255)
                                                                                                                                                            cacaCTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTAc                                                                                                                                                                                                                             <  5:296165/149‑1 (MQ=255)
                                                                                                                                                              cacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGc                                                                                                                                                                                                                           >  6:785728/1‑149 (MQ=255)
                                                                                                                                                               acTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                                                                                                                                                                          <  6:211637/149‑1 (MQ=255)
                                                                                                                                                               acTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                                                                                                                                                                          >  7:373406/1‑149 (MQ=255)
                                                                                                                                                                  ggtggGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtg                                                                                                                                                                                                                       <  6:330495/149‑1 (MQ=255)
                                                                                                                                                                    tggGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgtt                                                                                                                                                                                                                     >  6:1091923/1‑149 (MQ=255)
                                                                                                                                                                      ggAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                                                                                                                                                                                   <  6:1524413/149‑1 (MQ=255)
                                                                                                                                                                      ggAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                                                                                                                                                                                   <  8:246527/149‑1 (MQ=255)
                                                                                                                                                                       gAATTGTACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAGCCGATTAATACTCTCTCCGCTACGCAGTGTTTTa                                                                                                                                                                                                                  <  7:1122005/149‑1 (MQ=255)
                                                                                                                                                                          ttGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                                                                                                                                                  >  7:1430550/1‑146 (MQ=255)
                                                                                                                                                                          ttGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                                                                                                                                                                               <  7:542167/149‑1 (MQ=255)
                                                                                                                                                                             aaCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCaa                                                                                                                                                                                                            >  7:798000/1‑149 (MQ=255)
                                                                                                                                                                                tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtg                                                                                                                                                                                                                       <  5:1566278/135‑1 (MQ=255)
                                                                                                                                                                                    atTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATa                                                                                                                                                                                                                                          <  8:1580661/112‑1 (MQ=255)
                                                                                                                                                                                    atTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATa                                                                                                                                                                                                                                          >  7:1580661/1‑112 (MQ=255)
                                                                                                                                                                                     tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCAc                                                                                                                                                                                                    >  8:832552/1‑149 (MQ=255)
                                                                                                                                                                                            ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                                                                                                                                                                             <  8:1021025/149‑1 (MQ=255)
                                                                                                                                                                                            ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                                                                                                                                                                             <  7:248157/149‑1 (MQ=255)
                                                                                                                                                                                             cgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                    >  3:36261/1‑93 (MQ=255)
                                                                                                                                                                                               ccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGTATCATAGTCACTTGATGTGTCTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTTCGCANTTTCATTGAAATCATTACCTGCTAGTGATTATTTCAACCTACTGACTTTCATTTAATTTTTTTCACTCTATGGCAAATTTNCCATTTCAAACATTACTTTTGCTTTNATTCTCCGTTATCAGGTTTAATGTTTTTAAAGt                                                        >  3:163975/1‑279 (MQ=255)
                                                                                                                                                                                                cAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaaa                                                                                                                                                                                         <  6:1213143/149‑1 (MQ=255)
                                                                                                                                                                                                  gCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaata                                                                                                                                                                                       <  5:1651915/149‑1 (MQ=255)
                                                                                                                                                                                                   ccTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                    <  2:523400/87‑1 (MQ=255)
                                                                                                                                                                                                    cTGTTCAGCTTGGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACGTGAGGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGAACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCAGGGCTGATATGTTCCGTAGGCAGGTTTAATGTTTTAAAAGTg                                                       <  2:385465/275‑1 (MQ=255)
                                                                                                                                                                                                       ttCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                                                                                                                                                                                  >  8:1621970/1‑149 (MQ=255)
                                                                                                                                                                                                        tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                                                                                                                                                                                 >  8:1008628/1‑149 (MQ=255)
                                                                                                                                                                                                         cAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCt                                                                                                                                                                                <  7:682641/149‑1 (MQ=255)
                                                                                                                                                                                                                gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg                                                                                                                                                                         <  7:1159234/149‑1 (MQ=255)
                                                                                                                                                                                                                  cAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCt                                                                                                                                                                                >  5:968336/1‑140 (MQ=255)
                                                                                                                                                                                                                  cAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCt                                                                                                                                                                                <  6:968336/140‑1 (MQ=255)
                                                                                                                                                                                                                        gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTc                                                                                                                                                                 >  7:205347/1‑149 (MQ=255)
                                                                                                                                                                                                                            gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAAt                                                                                                                                                                                                                                           <  6:1545331/71‑1 (MQ=255)
                                                                                                                                                                                                                            gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCACTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAcc                                                                                                                                                             >  5:458105/1‑149 (MQ=255)
                                                                                                                                                                                                                                agcagcAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACt                                                                                                                                                         >  6:1331511/1‑149 (MQ=255)
                                                                                                                                                                                                                                 gcagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                                                                                                                                                                                   >  5:346597/1‑90 (MQ=255)
                                                                                                                                                                                                                                   agcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGaa                                                                                                                                                      <  8:205347/149‑1 (MQ=255)
                                                                                                                                                                                                                                   agcaATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGaa                                                                                                                                                      >  8:929230/1‑149 (MQ=255)
                                                                                                                                                                                                                                     caATAATTCCCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAttttttt                                                                                                                                     <  1:37743/164‑1 (MQ=255)
                                                                                                                                                                                                                                      aataatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAAttt                                                                                                                                                   >  6:931162/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCAAGTGATTAttt                                                                                                                                                                  <  5:1272502/131‑1 (MQ=255)
                                                                                                                                                                                                                                           ttACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                                                                                                                                                                             <  8:769934/102‑1 (MQ=255)
                                                                                                                                                                                                                                           ttACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATc                                                                                                                                               <  5:208555/148‑1 (MQ=255)
                                                                                                                                                                                                                                                   gTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAACCAt                                                                                                      >  2:715318/1‑181 (MQ=255)
                                                                                                                                                                                                                                                    ttGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAttttttt                                                                                                                                     <  8:494550/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     tGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTaa  <  1:125218/279‑1 (MQ=255)
                                                                                                                                                                                                                                                       tCACCCGGATCATAGTCACTTGANGTGANTATGATCCTATTAAGACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGAGTATTTAAACCCACGGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATc                                                                                                  <  4:715313/181‑1 (MQ=255)
                                                                                                                                                                                                                                                           tcGGATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTa                                                                                                                              >  8:36818/2‑149 (MQ=255)
                                                                                                                                                                                                                                                           ccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                                                                                                                                                                                   <  5:1477223/96‑1 (MQ=255)
                                                                                                                                                                                                                                                           ccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGt                                                                                                                                                                          <  8:1377847/105‑1 (MQ=255)
                                                                                                                                                                                                                                                           ccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGt                                                                                                                                                                          >  7:1377847/1‑105 (MQ=255)
                                                                                                                                                                                                                                                           ccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTa                                                                                                                              >  5:1206629/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             ggATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATg                                                                                                                            <  4:37559/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                tCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                                                                                                                                                                                   <  8:677345/59‑1 (MQ=255)
                                                                                                                                                                                                                                                                tCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                                                                                                                                                                            <  7:1021190/98‑1 (MQ=255)
                                                                                                                                                                                                                                                                tCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCa                                                                                                                         >  8:1225869/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aTAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCaaa                                                                                                                       >  7:933306/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAAt                                                                                                                      >  6:1608059/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   tAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAAt                                                                                                                      >  8:1425620/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     gTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTa                                                                                                                    <  7:1008628/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTa                                                                                                                    >  5:515061/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAg                                                                                                                   >  7:760221/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGc                                                                                             <  3:23761/168‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ttGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCaaa                                                                                                                       <  6:369143/141‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCAttt                                                                                                             <  5:736620/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                               
TTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTAA  >  NC_000913/3657333‑3657856

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