Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,603 T→A 18.3% intergenic (+36/+383) 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,6030TA18.3% 220.2 / 42.1 85intergenic (+36/+383)hdeD/arrSacid‑resistance membrane protein/antisense sRNA ArrS, function unknown
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base A (8/7);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-02

CGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTAAATAAGTAATCCGGGTT  >  NC_000913/3657330‑3657872
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                             
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)
                                                                                                                        ccTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAtttttt                                                                                                                                                      >  4:329481/1‑275 (MQ=255)
                                                                                                                        ccTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAtttttt                                                                                                                                                      >  2:452852/1‑275 (MQ=255)
                                                                                                                                 aaCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                                                                                                                                                                                            <  8:886642/148‑1 (MQ=255)
                                                                                                                                 aaCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTAt                                                                                                                                             >  4:395948/1‑275 (MQ=255)
                                                                                                                                               cTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATg                                                                                                                                                                                                                                                                      <  5:881836/140‑1 (MQ=255)
                                                                                                                                               cTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATg                                                                                                                                                                                                                                                                      >  6:881836/1‑140 (MQ=255)
                                                                                                                                                 gTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACt                                                                                                                                                                                                                                                                         >  6:827868/1‑135 (MQ=255)
                                                                                                                                                 gTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATAGTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACt                                                                                                                                                                                                                                                                         <  5:827868/135‑1 (MQ=255)
                                                                                                                                                     cGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                                                                                                                                                        >  5:497974/1‑132 (MQ=255)
                                                                                                                                                     cGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                                                                                                                                                        <  6:497974/132‑1 (MQ=255)
                                                                                                                                                         ggTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                                    <  6:139842/132‑1 (MQ=255)
                                                                                                                                                         ggTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                                                                                                                                                    >  5:139842/1‑132 (MQ=255)
                                                                                                                                                                acacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                                                                                                                                                                                            >  8:498287/1‑117 (MQ=255)
                                                                                                                                                                acacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                                                                                                                                                                                            <  7:498287/117‑1 (MQ=255)
                                                                                                                                                                acacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATTATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTANTGATTATTTAAACCTACTGATTTTCATCTAATTTTTTTCACTTCATGGCAATTTAGCCATTTCAAACATTATCATgt                                                                                                              >  3:133151/1‑274 (MQ=255)
                                                                                                                                                                acacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATgg                                                                                                              <  2:395949/275‑1 (MQ=255)
                                                                                                                                                                 cacTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATgg                                                                                                              <  4:446399/274‑1 (MQ=255)
                                                                                                                                                                  acTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                                                                                                                                                                                            >  5:258318/1‑115 (MQ=255)
                                                                                                                                                                  acTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                                                                                                                                                                                            <  6:258318/115‑1 (MQ=255)
                                                                                                                                                                            aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                                       <  5:586774/110‑1 (MQ=255)
                                                                                                                                                                            aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                                       <  5:1120126/110‑1 (MQ=255)
                                                                                                                                                                            aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                                       >  6:1120126/1‑110 (MQ=255)
                                                                                                                                                                            aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                                                                                                                                                       >  6:586774/1‑110 (MQ=255)
                                                                                                                                                                             ttGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCg                                                                                                 <  2:479664/275‑1 (MQ=255)
                                                                                                                                                                                    gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACt                                                                                                                                                                                                                                                                         >  8:612495/1‑100 (MQ=255)
                                                                                                                                                                                    gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACt                                                                                                                                                                                                                                                                         <  7:612495/100‑1 (MQ=255)
                                                                                                                                                                                     atatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAgg                                                                                         >  4:97068/1‑275 (MQ=255)
                                                                                                                                                                                      tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACGCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACt                                                                                                                                                                                                                                                                         >  5:1084507/1‑98 (MQ=39)
                                                                                                                                                                                      tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACt                                                                                                                                                                                                                                                                         <  6:1084507/98‑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)
                                                                                                                                                                                                                 ttgCCAGTTTGTTCGTTAAGCAGCACTAATTACCCCGGTTGTAACCCGGAGCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTGTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGaa                                                              <  4:339436/272‑1 (MQ=255)
                                                                                                                                                                                                                     cattttGTTCGTGAAGCAGCAATACTTACCCCGTTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGa                                                         <  1:133151/272‑1 (MQ=255)
                                                                                                                                                                                                                     cAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                                                                                                                                                                                                <  8:540970/76‑1 (MQ=255)
                                                                                                                                                                                                                     cAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                                                                                                                                                                                                >  7:540970/1‑76 (MQ=255)
                                                                                                                                                                                                                         ttGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAg                                                                                                                                                                                                                                         >  8:436385/1‑95 (MQ=255)
                                                                                                                                                                                                                         ttGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAg                                                                                                                                                                                                                                         <  7:436385/95‑1 (MQ=255)
                                                                                                                                                                                                                                ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                                                                                                                                                                                                 >  8:1109921/1‑64 (MQ=255)
                                                                                                                                                                                                                                ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                                                                                                                                                                                                 <  7:1109921/64‑1 (MQ=255)
                                                                                                                                                                                                                                     cagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                                      >  5:1023332/1‑70 (MQ=255)
                                                                                                                                                                                                                                     cagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                                                                                                                                                                                                      <  6:1023332/70‑1 (MQ=255)
                                                                                                                                                                                                                                                  cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTAt                            >  2:264901/1‑275 (MQ=255)
                                                                                                                                                                                                                                                  cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTAt                            >  2:70455/1‑275 (MQ=255)
                                                                                                                                                                                                                                                   ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                                                                                                                                                                                                           >  5:223188/1‑83 (MQ=255)
                                                                                                                                                                                                                                                   ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                                                                                                                                                                                                           <  6:223188/83‑1 (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)
                                                                                                                                                                                                                                                                              tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                                                                                                                                                                                                          >  5:921179/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                              tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                                                                                                                                                                                                          <  6:921179/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tgtgACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTa                                                                                                                                              >  8:534892/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                                 tgtgACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTa                                                                                                                                              <  7:534892/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                             
CGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGATATTCCTTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCATTTCAAACATTATCATGGCTGATATTTTCCGTAGTCAGGTTTAATGTTTTAAAAGTGCTGTGGGAAAGTGAACAAAGAGTTCCGTAAGCGTTGATGCTATGGGCGGTTAAATAAGTAATCCGGGTT  >  NC_000913/3657330‑3657872

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