Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,877,011 G→A intergenic (‑59/+56) yidA ← / ← yidB sugar phosphate phosphatase/DUF937 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,877,0110GA98.5% 237.4 / NA 68intergenic (‑59/+56)yidA/yidBsugar phosphate phosphatase/DUF937 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/34);  minor base C (0/1);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGT  >  NC_000913/3876866‑3877151
                                                                                                                                                 |                                                                                                                                             
aaTCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATac                                                                                                                                            <  6:23534/149‑1 (MQ=255)
 aTCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATaca                                                                                                                                           <  7:8577/149‑1 (MQ=255)
    gCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATacaaca                                                                                                                                        <  7:85487/149‑1 (MQ=255)
     cATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAacac                                                                                                                                       >  2:37836/1‑149 (MQ=255)
      aTTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACAcc                                                                                                                                      <  8:44758/149‑1 (MQ=255)
      aTTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACAcc                                                                                                                                      >  7:44758/1‑149 (MQ=255)
            aaCGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAgg                                                                                                                                 <  8:2685/148‑1 (MQ=255)
                       aaaTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGAt                                                                                                                    <  1:37836/149‑1 (MQ=255)
                                 aTCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAAcgcgcg                                                                                                          <  5:92128/149‑1 (MQ=255)
                                     ggCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGc                                                                                                      <  6:87711/149‑1 (MQ=255)
                                      gcagcaGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACAcc                                                                                                                                      <  4:126782/117‑1 (MQ=255)
                                      gcagcaGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACAcc                                                                                                                                      >  3:126782/1‑117 (MQ=255)
                                             aaGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTccc                                                                                              >  8:16810/1‑149 (MQ=255)
                                                ggTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCttt                                                                                           >  1:60321/1‑149 (MQ=255)
                                                  tGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGAGCACGCAACGCGCGTTGCATCCCCCttttc                                                                                         <  6:129688/149‑2 (MQ=38)
                                                  tGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCttttt                                                                                         <  2:60321/149‑1 (MQ=255)
                                                  tGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCttttt                                                                                         >  5:7723/1‑149 (MQ=255)
                                                     catccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAg                                                                                                                           <  7:98582/113‑1 (MQ=255)
                                                     catccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAg                                                                                                                           >  8:98582/1‑113 (MQ=255)
                                                        ccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGa                                                                                                                     <  6:104456/115‑1 (MQ=255)
                                                        ccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGa                                                                                                                     >  5:104456/1‑115 (MQ=255)
                                                             tCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCgg                                                                                >  6:4117/1‑147 (MQ=255)
                                                                     aaTGAGTTTAATAGCCATGAGTTCCCCATGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTcc                                                                      <  5:77890/149‑1 (MQ=255)
                                                                     aaTGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTcc                                                                      >  1:74440/1‑149 (MQ=255)
                                                                      aTGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCtt                                                                                     <  3:66035/133‑1 (MQ=255)
                                                                      aTGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCtt                                                                                     >  4:66035/1‑133 (MQ=255)
                                                                      aTGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTccc                                                                     >  5:143451/1‑149 (MQ=255)
                                                                          gTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCttt                                                                  >  1:125860/1‑148 (MQ=255)
                                                                           tttAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTg                                                                                        >  6:83790/1‑125 (MQ=255)
                                                                           tttAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCttt                                                                  >  3:10938/1‑147 (MQ=255)
                                                                           tttAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTc                                                                 <  6:7723/148‑1 (MQ=255)
                                                                                   ccATGAGTTCCCCAGGTAGATTGGTTCCTGCCCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTCGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTc                                                         >  5:19491/1‑148 (MQ=255)
                                                                                   ccATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTcagc                                                              >  7:355/1‑143 (MQ=255)
                                                                                   ccATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTcagc                                                              <  8:355/143‑1 (MQ=255)
                                                                                   ccATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTc                                                         >  5:143516/1‑148 (MQ=255)
                                                                                     aTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCAt                                                      <  5:4117/149‑1 (MQ=255)
                                                                                     aTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCttt                                                                  >  2:1537/1‑136 (MQ=255)
                                                                                     aTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCttt                                                                  <  1:1537/136‑1 (MQ=255)
                                                                                       gAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAgg                                                                                                                          >  1:97570/1‑80 (MQ=255)
                                                                                       gAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAgg                                                                                                                          <  2:97570/80‑1 (MQ=255)
                                                                                         gTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCt                                                  >  7:53695/1‑149 (MQ=255)
                                                                                          ttCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTg                                                 <  2:125860/149‑1 (MQ=255)
                                                                                             cccAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAaga                                                                            >  3:61714/1‑119 (MQ=255)
                                                                                             cccAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAaga                                                                            <  4:61714/119‑1 (MQ=255)
                                                                                              ccAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCgg                                              >  7:38674/1‑148 (MQ=255)
                                                                                                   tAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCa                                        <  7:16810/149‑1 (MQ=255)
                                                                                                         ggTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCAt                                                      <  7:53788/129‑1 (MQ=255)
                                                                                                         ggTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCAt                                                      >  8:53788/1‑129 (MQ=255)
                                                                                                          gTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGt                                  >  2:113393/1‑148 (MQ=255)
                                                                                                           ttCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCa                                        <  3:125728/141‑1 (MQ=255)
                                                                                                           ttCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCa                                        >  4:125728/1‑141 (MQ=255)
                                                                                                             ccTGACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTAtt                              >  4:117464/1‑149 (MQ=255)
                                                                                                                gACCCATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGc                           >  5:71951/1‑149 (MQ=255)
                                                                                                                  cccATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCtt                         <  6:143516/149‑1 (MQ=255)
                                                                                                                    cATGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTtgt                       <  6:71951/149‑1 (MQ=255)
                                                                                                                     aTGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCt                                                                                                 >  7:69348/1‑74 (MQ=255)
                                                                                                                     aTGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCt                                                                                                 <  8:69348/74‑1 (MQ=255)
                                                                                                                     aTGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTtgtg                      <  2:74440/149‑1 (MQ=255)
                                                                                                                     aTGCTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTtgtg                      <  8:53695/149‑1 (MQ=255)
                                                                                                                        cTAACGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTgtggt                    <  3:117464/148‑1 (MQ=255)
                                                                                                                          aaCGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGa                  >  1:89699/1‑148 (MQ=255)
                                                                                                                          aaCGCGATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGa                  >  3:118020/1‑148 (MQ=255)
                                                                                                                             gcgATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTtgtg                      >  7:7903/1‑141 (MQ=255)
                                                                                                                             gcgATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTtgtg                      <  8:7903/141‑1 (MQ=255)
                                                                                                                              cgATTCCGCTCAAAAATCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTCGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACtt             <  6:19491/149‑1 (MQ=255)
                                                                                                                                         ataaaTCACTTCATCACCCGAGGGTAAAGGGGGATGCAACGCGCGTTGCGTGCTCCGTTTTTGGTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGGTTGTGGTGACACTTCACCTTGCGGt  <  7:12965/147‑1 (MQ=255)
                                                                                                                                         aaaaaTCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTtgtg                      <  2:19884/129‑1 (MQ=255)
                                                                                                                                         aaaaaTCAATACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTtgtg                      >  1:19884/1‑129 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
AATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGT  >  NC_000913/3876866‑3877151

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