Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,053,512 G→A intergenic (+94/+135) hemN → / ← yshB coproporphyrinogen III dehydrogenase/uncharacterized protein YshB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,053,5120GA97.6% 143.6 / ‑6.1 42intergenic (+94/+135)hemN/yshBcoproporphyrinogen III dehydrogenase/uncharacterized protein YshB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/20);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CAAAAGATGGGCTGGTGGATGTGGATGAGAAGGGAATACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACTGCGCCACTCGATACCAGATTAATTATTGAT  >  NC_000913/4053276‑4053749
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             
cAAAAGATGGGCTGGTGGATGTGGATGAGAAGGGAATACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCt                                                                                                                                                                                                                                 >  1:271878/1‑251 (MQ=255)
cAAAAGATGGGCTGGTGGATGTGGATGAGAAGGGAATACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGAGTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCt                                                                                                                                                                                                                                 <  1:315415/251‑1 (MQ=255)
    aGATGGGCTGGTGGATGTGGATGAGAAGGGAATACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTg                                                                                                                                                                                                                             >  1:213938/1‑251 (MQ=255)
        gggCTGGTGGATGTGGATGAGAAGGGAATACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAgg                                                                                                                                                                                                                         <  2:280624/251‑1 (MQ=255)
          gCTGGTGGATGTGGATGAGAAGGGAATACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACAATTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGcc                                                                                                                                                                                                                       >  1:206804/1‑251 (MQ=255)
                  atgtggatgAGAAGGGAATACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAgg                                                                                                                                                                                                               >  2:26113/1‑251 (MQ=255)
                    gtggatgAGAAGGGAATACAGGCGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACATATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAACGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCg                                                                                                                                                                                                             >  2:12080/1‑251 (MQ=255)
                                 gAATACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                                                                <  1:155986/251‑1 (MQ=255)
                                   aTACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                                                                >  2:148625/1‑249 (MQ=255)
                                           gACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTg                                                                                                                                                                                       <  2:225217/250‑1 (MQ=255)
                                                         gCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAt                                                                                                                                                                        <  2:105185/251‑1 (MQ=255)
                                                                       aaCATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGATGGTcgc                                                                                                                                                          >  1:252149/1‑251 (MQ=255)
                                                                       aaCATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTGTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                                                                                                                          <  2:37243/251‑1 (MQ=255)
                                                                               cATGTGTTTTGATTCCTATCTGCGCTAGAAAGCGCGGATGCAGCAGTTCTCTCGTGTGATTTAAATAGTGTATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGTTGAGAACAAGATTTCCTGACGCGCTGGGCAACAAAACCTATTGAAATGATTGTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGGTGCCTTATGCGGCTCTGGTCGCGTCTTAt                                                                                                                                                   <  2:340447/250‑1 (MQ=255)
                                                                                  gggCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCAGGTGATTTAAATAGTGAAGTGCGCTTCTTTTACAAAGTAGAAAACGAAAGCTCCATCAATGTGCTGAGAACATGAGTGCCTGATGCGCTGTGCAACAAATCCTATTGAAATTATTGCTTGTGTAGGCCGGATAAGGCGTGTACGCCGCATCAGGCATCTGAGCATGATGCCTGATTCGGCGCTGGTCATGTCTTATCAg                                                                                                                                                <  2:239736/248‑1 (MQ=255)
                                                                                            aCCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGAGTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                                                                                                     <  1:215983/251‑1 (MQ=255)
                                                                                                 tCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                                                                                                >  1:331270/1‑251 (MQ=255)
                                                                                                      gcCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                                                                                                                           <  1:315618/251‑1 (MQ=255)
                                                                                                                gcgGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCAGCATCAGGCATCTGTGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                                                                                                                 >  1:47067/1‑251 (MQ=255)
                                                                                                                  ggATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGtt                                                                                                                                                                                                           <  1:224379/159‑1 (MQ=255)
                                                                                                                  ggATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGtt                                                                                                                                                                                                           >  2:224379/1‑159 (MQ=255)
                                                                                                                      gcagcaGTTCTCTCGGGCGAATTACATAGTGAATGGCGCTTAGCTCACTAAGTAGAAAACGAAAGCGCAATCAATGTGCTGAGAACAAGATTGACTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGCTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGACATAAACAACg                                                                                                             >  2:388728/1‑249 (MQ=255)
                                                                                                                          catttCTCTCGGTTGATTGAAATAGTGAATGGCGCTTCGTTTACAAGGTAGAAAACGAAAGCGGCATCAATTTGCTGAGAACAAGATTGCCTGATTCGGTGGGCAACAAATCCTATTGTAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCAGCAGGCATCTGTGCATGATGCCTGATGCTGCGCTAGTCGCGTCTTATCATGCCTACGACTTTCCTCGCGTCCGGCATCAACAACGAGCTTc                                                                                                       <  2:53056/248‑1 (MQ=255)
                                                                                                                              tctctcGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGTCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc                                                                                                    <  2:126055/250‑1 (MQ=255)
                                                                                                                                   cGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaaa                                                                                                <  1:41927/249‑1 (MQ=255)
                                                                                                                                      gTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGcc                                                                                           >  1:203446/1‑251 (MQ=255)
                                                                                                                                            tAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATc                                                                                      >  1:137067/1‑250 (MQ=255)
                                                                                                                                                     aaTGGCGCTTCGTTTACAAAGTAGAAAACTAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCATACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGca                                                                             <  2:284958/250‑1 (MQ=255)
                                                                                                                                                        ggCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcacac                                                                          >  1:225226/1‑250 (MQ=255)
                                                                                                                                                                           aGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAg                                                       <  1:326333/250‑1 (MQ=255)
                                                                                                                                                                             aaaaCGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGttt                                                    <  1:42035/251‑1 (MQ=255)
                                                                                                                                                                                     aGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgtg                                            <  1:201518/251‑1 (MQ=255)
                                                                                                                                                                                        gcCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCCCGGCAGCAACGGCAGTTTGCGGTGTGTgg                                         >  2:121226/1‑251 (MQ=255)
                                                                                                                                                                                            tCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGCCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGt                                     >  1:373451/1‑251 (MQ=255)
                                                                                                                                                                                              aaTGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTc                                    >  1:52764/1‑250 (MQ=255)
                                                                                                                                                                                                            aaGATTGCCTGCTGCGCTGGACAACAAAACCTATTGAAATTATTGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTATTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACACCGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGCGTGTGGCTGTCAACTGCGCCAcc                        >  2:220028/1‑247 (MQ=255)
                                                                                                                                                                                                                        tGCGCTGGGCAACAGAACCTATTGAAATTATTGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAAAGCCCAATCATTGCGGCACACAGCACGGAAGCAACGGCAGTTTGCGGTGTGTGGCTGTCCACTGCGCCACTCGCTACAAGatnaat         >  2:369088/1‑251 (MQ=255)
                                                                                                                                                                                                                           gcTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGCGGCTGTCAACTGCGCCACTCGATACCAGATTAattat      <  1:140129/251‑1 (MQ=255)
                                                                                                                                                                                                                            cTGGGCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCAGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACTGCGCCACTCGATACCAGATTAattatt     >  1:243329/1‑251 (MQ=255)
                                                                                                                                                                                                                              gggCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACTTCGCCACTCGATACCAGATTAATTATTg    <  2:243329/250‑1 (MQ=255)
                                                                                                                                                                                                                              gggCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACTGCGCCACTCGATACCAGATTAattatt     <  1:124034/249‑1 (MQ=255)
                                                                                                                                                                                                                               ggCAACAAAACCTATTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACTGCGCCACTCGATACCAGATTAATTATTGAt  >  1:202891/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             
CAAAAGATGGGCTGGTGGATGTGGATGAGAAGGGAATACAGGTGACGGCGAAAGGTCGCTTGCTGATCCGCAACATTTGCATGTGCTTTGATACCTATCTGCGCCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACTGCGCCACTCGATACCAGATTAATTATTGAT  >  NC_000913/4053276‑4053749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: