Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,950,559 (C)5→6 intergenic (‑37/‑273) aspS ← / → yecD aspartyl‑tRNA synthetase/isochorismatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,5541.C97.9% 338.9 / ‑1.2 97intergenic (‑32/‑278)aspS/yecDaspartyl‑tRNA synthetase/isochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (48/47);  total (49/48)
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

Only 100 of 104 total aligned reads displayed.

ATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAC  >  NC_000913/1950422‑1950683
                                                                                                                                      |                                                                                                                                 
aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGcccccc                                                                                                                             >  2:683161/1‑138 (MQ=255)
aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGcccccc                                                                                                                             <  2:688164/139‑2 (MQ=255)
  cAGGCTACCAAGATCACGACGACGGTTGACCCAATCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCGATACGCATGAGATATCCCTTAACTTAGCTGCTGGCGTAATCCCCCctg                                                                                                                            >  1:2044530/1‑139 (MQ=255)
     gCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                         >  1:1241857/1‑139 (MQ=255)
       tACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTgcg                                                                                                                       <  1:1881585/139‑1 (MQ=255)
        aCCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTgcgc                                                                                                                      >  2:1728700/1‑139 (MQ=255)
               tCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGAc                                                                                                               >  1:1355818/1‑139 (MQ=255)
                   gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAg                                                                                                           >  1:968066/1‑139 (MQ=255)
                     cgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTc                                                                                                         <  1:1728700/139‑1 (MQ=255)
                          gTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                    >  1:1697250/1‑139 (MQ=255)
                            tGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                  <  2:1697250/139‑1 (MQ=255)
                             gACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGTTGCGCAGGTGACCAAGTc                                                                                                         <  1:647494/131‑1 (MQ=255)
                             gACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGTTGCGCAGGTGACCAAGTc                                                                                                         >  2:647494/1‑131 (MQ=255)
                                   aCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                          <  2:567549/139‑1 (MQ=255)
                                       cacaGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTgcg                                                                                                                       >  1:982711/1‑107 (MQ=255)
                                       cacaGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTgcg                                                                                                                       <  2:982711/107‑1 (MQ=255)
                                            agTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            >  2:182534/1‑129 (MQ=255)
                                            agTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  1:182534/129‑1 (MQ=255)
                                            agTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                  <  2:1241857/139‑1 (MQ=255)
                                            agTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                  >  2:1452778/1‑139 (MQ=255)
                                               cACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                       <  1:1429180/115‑1 (MQ=255)
                                               cACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                       >  2:1429180/1‑115 (MQ=255)
                                                 ccTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGccc                                                                                                                                >  1:2051257/1‑87 (MQ=255)
                                                 ccTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGccc                                                                                                                                <  2:2051257/87‑1 (MQ=255)
                                                    gctgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                         >  1:530820/1‑92 (MQ=255)
                                                    gctgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                         <  2:530820/92‑1 (MQ=255)
                                                      tgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGTCCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGaaa                                                                        >  2:1045133/1‑139 (MQ=255)
                                                        ccccACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGccc                                                                                                                                <  2:940861/80‑1 (MQ=255)
                                                        ccccACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGccc                                                                                                                                >  1:940861/1‑80 (MQ=255)
                                                          ccACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAg                                                                    >  2:650868/1‑139 (MQ=255)
                                                            aCGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  1:683161/139‑1 (MQ=255)
                                                            aCGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:968066/139‑1 (MQ=255)
                                                            aCGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:1355818/139‑1 (MQ=255)
                                                             cGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                       <  2:1576137/101‑1 (MQ=255)
                                                             cGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                       >  1:1576137/1‑101 (MQ=255)
                                                                            gCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCg                                                                                                        >  2:1103721/1‑85 (MQ=255)
                                                                            gCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCg                                                                                                        <  1:1103721/85‑1 (MQ=255)
                                                                                tCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:1480123/119‑1 (MQ=255)
                                                                                tCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  1:1480123/1‑119 (MQ=255)
                                                                                 ccACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAGGCGGCTATTATACTGGAAATTCTg                                            >  1:7663/1‑139 (MQ=255)
                                                                                   acaATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGACACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGc                                           <  1:1452778/139‑1 (MQ=255)
                                                                                       tatTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGTCCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCa                                       <  1:1045133/139‑1 (MQ=255)
                                                                                       tatTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  1:462257/1‑112 (MQ=255)
                                                                                       tatTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:462257/112‑1 (MQ=255)
                                                                                          tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  1:573849/83‑1 (MQ=255)
                                                                                          tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            >  2:573849/1‑83 (MQ=255)
                                                                                          tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTACGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCg                                    >  2:1845083/1‑139 (MQ=255)
                                                                                            tGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  1:744981/1‑107 (MQ=255)
                                                                                            tGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:744981/107‑1 (MQ=255)
                                                                                               aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                         >  1:947121/1‑49 (MQ=255)
                                                                                               aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                         <  2:947121/49‑1 (MQ=255)
                                                                                               aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                  >  2:1273812/1‑72 (MQ=255)
                                                                                               aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                  <  1:1273812/72‑1 (MQ=255)
                                                                                               aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                             >  2:324753/1‑109 (MQ=255)
                                                                                               aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                             <  1:324753/109‑1 (MQ=255)
                                                                                               aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAaga                               >  2:181429/1‑139 (MQ=255)
                                                                                                 gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGt                                                                                                                  >  2:500545/1‑54 (MQ=255)
                                                                                                 gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGt                                                                                                                  <  1:500545/54‑1 (MQ=255)
                                                                                                   aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                  >  2:1021025/1‑132 (MQ=255)
                                                                                                   aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                  <  1:1021025/132‑1 (MQ=255)
                                                                                                   aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCt                           >  2:1185624/1‑139 (MQ=255)
                                                                                                          atCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtat                                                            <  2:951586/99‑1 (MQ=255)
                                                                                                          atCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtat                                                            >  1:951586/1‑99 (MQ=255)
                                                                                                           tCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            >  2:1654615/1‑66 (MQ=255)
                                                                                                           tCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  1:1654615/66‑1 (MQ=255)
                                                                                                             ccTTAACTTAGCTGCCGGCGGATGCCCGCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCCGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCATGGCCCGCGCTg                >  1:1229086/1‑139 (MQ=255)
                                                                                                               ttAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAg                                                                    <  2:1455861/86‑1 (MQ=255)
                                                                                                               ttAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAg                                                                    >  1:1455861/1‑86 (MQ=255)
                                                                                                                tAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            >  2:64487/1‑61 (MQ=255)
                                                                                                                tAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  1:64487/61‑1 (MQ=255)
                                                                                                                 aaCTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:1187214/86‑1 (MQ=255)
                                                                                                                 aaCTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGc                        >  2:99647/1‑128 (MQ=255)
                                                                                                                 aaCTTAGCTGCCGGCGGATGCCCCCCTGCCGCGCAGGTGACCAAGCCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  1:1187214/1‑86 (MQ=255)
                                                                                                                  aCTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                       >  1:461302/1‑48 (MQ=38)
                                                                                                                  aCTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                       <  2:461302/48‑1 (MQ=38)
                                                                                                                  aCTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAATCACAGCCTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCGGGAAc            <  1:384688/139‑1 (MQ=255)
                                                                                                                     tAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCg                                                                                                   <  1:1755640/49‑1 (MQ=38)
                                                                                                                     tAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCg                                                                                                   >  2:1755640/1‑49 (MQ=38)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTc                                                                                                         <  1:655517/42‑1 (MQ=37)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTc                                                                                                         >  2:655517/1‑42 (MQ=37)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               <  1:461076/52‑1 (MQ=38)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               >  2:461076/1‑52 (MQ=38)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                          <  2:1050776/57‑1 (MQ=39)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                          >  1:1050776/1‑57 (MQ=39)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggcgg                                                                <  2:2040844/83‑1 (MQ=255)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggcgg                                                                >  1:2040844/1‑83 (MQ=255)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                  >  1:1563808/1‑113 (MQ=255)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                  <  2:1563808/113‑1 (MQ=255)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCt        >  2:1498609/1‑139 (MQ=255)
                                                                                                                      aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCt        >  2:55917/1‑139 (MQ=255)
                                                                                                                       gCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTc                                                                                                         <  1:768979/41‑1 (MQ=37)
                                                                                                                       gCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTc                                                                                                         >  2:768979/1‑41 (MQ=37)
                                                                                                                       gCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                  <  1:1509812/48‑1 (MQ=38)
                                                                                                                       gCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                  >  2:1509812/1‑48 (MQ=38)
                                                                                                                       gCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTc       >  1:370437/1‑139 (MQ=255)
                                                                                                                         tGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAGGCGGCTATTATACTGGAAATTCTGCCG‑‑CCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAc  <  2:7663/139‑1 (MQ=255)
                                                                                                                           ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  1:1498609/139‑1 (MQ=255)
                                                                                                                           ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   >  1:145112/1‑139 (MQ=255)
                                                                                                                           ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  1:55917/139‑1 (MQ=255)
                                                                                                                           ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGACGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  1:1455044/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                 
ATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAC  >  NC_000913/1950422‑1950683

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