Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,599 T→A 14.8% 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,5990TA14.8% 296.9 / 36.5 101intergenic (+32/+387)hdeD/arrSacid‑resistance membrane protein/antisense sRNA ArrS, function unknown
Reads supporting (aligned to +/- strand):  ref base T (41/45);  new base A (8/7);  total (49/52)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 106 total aligned reads displayed.

TGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTAAATAAGTAATCCGGGTT  >  NC_000913/3657327‑3657872
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
tggcGGTTTGGTCTGTGTTGAAGGGGTTATTGTCCGTATGAGTTGGTATCGTCCGCGTTCAATGAAAGGTAGATGGCTACAGCTCGTTATTGGTGTGCTGGGTATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTg                                                                                                                                                                                                                                                                                 <  2:283590/272‑1 (MQ=255)
   gtgTTTTTTCTGGTTTGTCGGGGTTGTTCGGCTTATTAGGTGGTATCGTCAGCGGTCAATGAAAGGCAGCTGGCTTCTGCTGGTTATTGGGGTGCTTGATATCGTCATTGCCTTGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACGGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                                                                                                                                                                                              <  4:160761/273‑1 (MQ=255)
     gTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTCGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                                                                                                                                                                                            >  2:243780/1‑275 (MQ=255)
             tgggtTGCGGGGGTTATTTGCCTGAGGGGTTGGTTTCGTCAGCGTTCAATGAAAGGCAGCGGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACGCTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                                    <  4:243779/272‑1 (MQ=255)
                           tATTTGCCTGATGGGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGGTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                                      <  2:480510/275‑1 (MQ=255)
                           tATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCCCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCAGGAACATAGTCCCTTGATGTGAATATGATCCGATTAATActct                                                                                                                                                                                                                                                      >  1:69724/1‑275 (MQ=255)
                           tATTCGCCGGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                                      <  1:93406/275‑1 (MQ=255)
                                      tgggTTGGTAGCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACg                                                                                                                                                                                                                                            <  2:209817/271‑1 (MQ=255)
                                      tgcgtgtgTATCGCCAGTGTTCAATGTAATGCCGCTGTGTACTGCTCGTTATTGGTTTGCTGGATATAGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTTGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTAACCCGGTTATCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACg                                                                                                                                                                                                                                            <  4:355407/267‑1 (MQ=255)
                                        aGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                                                                                                                                                                                          <  4:157897/274‑1 (MQ=255)
                                                     tgCGTTCAATGAAAGGCCGGTGTTTACATNNTGTTTGTGTTGTGCTGGATATCGTCATTGCNGGGATATTCCTTGGTGCAACGACGATGTTGTCTGCAACGCTGGTTTCCACCCTGGTGGGAATTGAACTGATATTNAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCaa                                                                                                                                                                                                                            <  3:69724/274‑1 (MQ=255)
                                                       cGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCATTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTAGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCAATTGATGGTAATATGATCCGATTAATCCTCTCTCCGCTACGCAGTGTTGTCAACCAAtt                                                                                                                                                                                                                          >  2:77540/1‑275 (MQ=255)
                                                          tCAATGAAAGGGAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcg                                                                                                                                                                                                                       <  2:393118/275‑1 (MQ=255)
                                                                        cTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaaa                                                                                                                                                                                                         >  2:446400/1‑275 (MQ=255)
                                                                           ggtACAGCTTGTTTATGGGGTGCTGGATATCGTCATTGTCTGGGTATTCCTTGGTGGAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCGGCCGGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                                                                                                                                                                                                      <  4:461451/273‑1 (MQ=255)
                                                                                            gtgtGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTa                                                                                                                                                                                                                                                                        >  2:328353/1‑192 (MQ=255)
                                                                                                gCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                                    <  4:328352/192‑1 (MQ=255)
                                                                                                   ggATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACt                                                                                                                                                                                                                                                                                   <  2:242426/174‑1 (MQ=255)
                                                                                                                         ttCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACt                                                                                                                                                                                                                                                                                   <  4:228940/152‑1 (MQ=255)
                                                                                                                           ccTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAtttttt                                                                                                                                                      >  4:329481/1‑275 (MQ=255)
                                                                                                                             ttGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACt                                                                                                                                                                                                                                                                                   <  7:218643/148‑1 (MQ=255)
                                                                                                                              tGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACtt                                                                                                                                                                                                                                                                                  >  6:261082/1‑148 (MQ=255)
                                                                                                                                gTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                <  5:261082/148‑1 (MQ=255)
                                                                                                                                gTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                <  6:1085876/148‑1 (MQ=255)
                                                                                                                                gTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                >  7:567678/1‑148 (MQ=255)
                                                                                                                                gTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                <  8:567678/148‑1 (MQ=255)
                                                                                                                                gTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACNTGa                                                                                                                                                                                                                                                                                <  7:1018016/148‑1 (MQ=255)
                                                                                                                                    aaCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTAt                                                                                                                                             >  4:395948/1‑275 (MQ=255)
                                                                                                                                                  cTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATg                                                                                                                                                                                                                                                                      >  6:881836/1‑140 (MQ=255)
                                                                                                                                                  cTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATg                                                                                                                                                                                                                                                                      <  5:881836/140‑1 (MQ=255)
                                                                                                                                                    gTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACt                                                                                                                                                                                                                                                                         >  6:827868/1‑135 (MQ=255)
                                                                                                                                                    gTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATAGTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACt                                                                                                                                                                                                                                                                         <  5:827868/135‑1 (MQ=255)
                                                                                                                                                        cGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                                                                                                                                                        >  5:497974/1‑132 (MQ=255)
                                                                                                                                                            ggTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                                    >  5:139842/1‑132 (MQ=255)
                                                                                                                                                            ggTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                                    <  6:139842/132‑1 (MQ=255)
                                                                                                                                                                   acacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                                                                                                                                                                                            <  7:498287/117‑1 (MQ=255)
                                                                                                                                                                   acacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                                                                                                                                                                                            >  8:498287/1‑117 (MQ=255)
                                                                                                                                                                   acacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATTATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTANTGATTATTTAAACCTACTGATTTTCATCTAATTTTTTTCACTTCATGGCAATTTAGCCATTTCAAACATTATCATgt                                                                                                              >  3:133151/1‑274 (MQ=255)
                                                                                                                                                                   acacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATgg                                                                                                              <  2:395949/275‑1 (MQ=255)
                                                                                                                                                                    cacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATgg                                                                                                              <  4:446399/274‑1 (MQ=255)
                                                                                                                                                                     acTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                                                                                                                                                                                            <  6:258318/115‑1 (MQ=255)
                                                                                                                                                                              aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                <  5:707076/102‑1 (MQ=255)
                                                                                                                                                                              aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                >  6:707076/1‑102 (MQ=255)
                                                                                                                                                                               aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                                       <  5:586774/110‑1 (MQ=255)
                                                                                                                                                                               aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                                       >  6:586774/1‑110 (MQ=255)
                                                                                                                                                                               aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                                       >  6:1120126/1‑110 (MQ=255)
                                                                                                                                                                               aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                                       <  5:1120126/110‑1 (MQ=255)
                                                                                                                                                                                ttGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCg                                                                                                 <  2:479664/275‑1 (MQ=255)
                                                                                                                                                                                       gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACt                                                                                                                                                                                                                                                                         >  8:612495/1‑100 (MQ=255)
                                                                                                                                                                                        atatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAgg                                                                                         >  4:97068/1‑275 (MQ=255)
                                                                                                                                                                                         tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACGCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACt                                                                                                                                                                                                                                                                         >  5:1084507/1‑98 (MQ=39)
                                                                                                                                                                                         tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACt                                                                                                                                                                                                                                                                         <  6:1084507/98‑1 (MQ=255)
                                                                                                                                                                                           tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                >  7:376128/1‑89 (MQ=255)
                                                                                                                                                                                           tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                <  8:376128/89‑1 (MQ=255)
                                                                                                                                                                                            ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                >  7:766812/1‑88 (MQ=255)
                                                                                                                                                                                            ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                                                                                                                                                <  8:766812/88‑1 (MQ=255)
                                                                                                                                                                                              aaacaggtcAGCCTGTTCAGCTTCGGTNTNTTGTTCGTTAAGCATCAATAATTACCCCGGTTGTCACCCGGATCANAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCANTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTaa                                                                                    <  3:52864/266‑1 (MQ=255)
                                                                                                                                                                                                             ttCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGa                                                                                                                                                                                                                                                                     >  6:821994/1‑82 (MQ=255)
                                                                                                                                                                                                             ttCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGa                                                                                                                                                                                                                                                                     <  5:821994/82‑1 (MQ=255)
                                                                                                                                                                                                              tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                                                                                                                                                        >  7:234000/1‑78 (MQ=38)
                                                                                                                                                                                                              tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                                                                                                                                                        <  8:234000/78‑1 (MQ=38)
                                                                                                                                                                                                               cAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCCAACATTATCATTGCTGATATTTTCCGTAGTCAGTTTTAATTTTTTaaaacctgttt                                                                   >  1:52864/1‑267 (MQ=255)
                                                                                                                                                                                                                        cattttGTTCGTGAAGCAGCAATACTTACCCCGTTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGa                                                         <  1:133151/272‑1 (MQ=255)
                                                                                                                                                                                                                        cAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                                                                                                                                                                                                <  8:540970/76‑1 (MQ=255)
                                                                                                                                                                                                                        cAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                                                                                                                                                                                                >  7:540970/1‑76 (MQ=255)
                                                                                                                                                                                                                            ttGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAg                                                                                                                                                                                                                                         <  7:436385/95‑1 (MQ=255)
                                                                                                                                                                                                                            ttGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAg                                                                                                                                                                                                                                         >  8:436385/1‑95 (MQ=255)
                                                                                                                                                                                                                                   ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                                                                                                                                                                                                 >  8:1109921/1‑64 (MQ=255)
                                                                                                                                                                                                                                   ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                                                                                                                                                                                                 <  7:1109921/64‑1 (MQ=255)
                                                                                                                                                                                                                                        cagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                                      <  6:1023332/70‑1 (MQ=255)
                                                                                                                                                                                                                                        cagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                                      >  5:1023332/1‑70 (MQ=255)
                                                                                                                                                                                                                                                     cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTAt                            >  2:264901/1‑275 (MQ=255)
                                                                                                                                                                                                                                                     cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTAt                            >  2:70455/1‑275 (MQ=255)
                                                                                                                                                                                                                                                      ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                                                                                                                                                                                                           <  6:223188/83‑1 (MQ=255)
                                                                                                                                                                                                                                                      ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                                                                                                                                                                                                           >  5:223188/1‑83 (MQ=255)
                                                                                                                                                                                                                                                         gTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACt                                                                                                                                                                         >  4:234551/1‑130 (MQ=255)
                                                                                                                                                                                                                                                         gTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTc                                              >  2:378399/1‑253 (MQ=255)
                                                                                                                                                                                                                                                         gTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGc                        >  2:39122/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAAt                                                                                                                                                                     <  2:234552/130‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTa                                          <  4:378398/253‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ggACCATAGTCACTTGACGTGACTATGATACGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACATGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTAAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTGATGTTTTAAAAGTGCTGTGGGAAAGTGAAAAAAGAGTACCGTAAACGTTGTAGGTATGGGCGGTAAAATaa              >  4:211671/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                    gATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATc                                                                                                                                                                                                                >  2:224913/1‑80 (MQ=255)
                                                                                                                                                                                                                                                                    gATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaat                                                                                                                                                                                                        >  4:490437/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                        aTAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                                                                                                                                                                                            <  4:224912/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aTAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAAtt                                                                                                                                                                                                    <  2:490438/88‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCaaa                                                                                                                                       >  8:537798/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                         tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTAAATAAGTAATc        <  4:70454/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcg                                                                                                                                                                                                                       >  7:1128534/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                          aGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcg                                                                                                                                                                                                                       <  8:1128534/67‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt                                                                                                                                                                                             <  7:366220/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt                                                                                                                                                                                             >  8:366220/1‑93 (MQ=39)
                                                                                                                                                                                                                                                                          aGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTa                                                                                                                                                                           >  8:4664/1‑111 (MQ=255)
                                                                                                                                                                                                                                                                           gTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTc                                                                                                                                                                  >  6:98526/1‑119 (MQ=255)
                                                                                                                                                                                                                                                                           gTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTc                                                                                                                                                                  <  5:98526/119‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGa                                                                                                                                                                                                                                <  5:239757/56‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGa                                                                                                                                                                                                                                >  6:239757/1‑56 (MQ=255)
                                                                                                                                                                                                                                                                              aCTTGATGTGACTATGATCCGAATAATACTCTCCCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCCTTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCAATTTCACACATTATCATGTCTGATATTTTCCGTAGTCAGGTGTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACACAGAGTTCCGTAAGCGTTGATCCTATGGGCGGTTAAATAAGTAATCCggtt   >  4:346263/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                               cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataat                                                                                                                                                                                                     >  7:830549/1‑80 (MQ=255)
                                                                                                                                                                                                                                                                               cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataat                                                                                                                                                                                                     <  8:830549/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTTTGGGAAAGTGAACAAAGAGTTCCGTACGCGTTGTTGCTATGGGCGTTTAAATAAGTAACCCGGGtt  >  4:377205/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
TGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTAAATAAGTAATCCGGGTT  >  NC_000913/3657327‑3657872

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