Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	160,736	A→T	100%	intergenic (‑132/+46)	dksA ← / ← sfsA	transcriptional regulator of rRNA transcription; DnaK suppressor protein/sugar fermentation stimulation protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	160,736	0	A	T	93.5%	263.9 / 6.7	91	intergenic (‑132/+46)	dksA/sfsA	transcriptional regulator of rRNA transcription; DnaK suppressor protein/sugar fermentation stimulation protein A
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base T (40/43);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAT  >  NC_000913/160639‑160832                                                                                                  |                                                                                                 tGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGtttt                                                                                                >  1:977727/1‑100 (MQ=255)   ttCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTg                                                                                              <  1:905523/100‑1 (MQ=255)    tCGGGGGAAAAGTGAGGCCGCTATAAATAGCAGATGGTTTTTCGGAGAGCGATGATCTCTACGTAACACTTGGCAGCTGTGTGAGCAAAAGCGTCTTTg                                                                                              <  2:459539/99‑1 (MQ=255)    tCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTgc                                                                                             >  1:1293962/1‑100 (MQ=255)     cGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTgcg                                                                                            <  2:1551172/100‑1 (MQ=255)       ggggAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCa                                                                                          >  2:1031828/1‑100 (MQ=255)        gggAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCAc                                                                                         >  2:1986551/1‑100 (MQ=255)        gggAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCAc                                                                                         >  1:1090496/1‑100 (MQ=255)           aaaaaTGAGGCCGCTATAAATAGCAGATGCTTTTCCCGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCAcgcg                                                                                      <  1:1051612/100‑1 (MQ=255)            aaaaTGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGa                                                                                     >  1:1752570/1‑100 (MQ=255)                  aGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTCTTTGCGCACGCGACCAGaa                                                                               >  2:2073357/1‑100 (MQ=255)                   ggCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCATGCGACCAGAAt                                                                              >  2:1145330/1‑100 (MQ=255)                          ataAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAAtt                                                                       <  2:901601/100‑1 (MQ=255)                           taAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTa                                                                      <  1:888669/100‑1 (MQ=255)                             aaTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTacc                                                                    >  1:357439/1‑98 (MQ=255)                             aaTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTacc                                                                    >  1:355071/1‑98 (MQ=255)                              aTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAAcc                                                                   <  2:2420294/100‑1 (MQ=255)                              aTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGAAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAAcc                                                                   <  2:2508237/100‑1 (MQ=255)                               tAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAGCACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCa                                                                  <  2:1392165/100‑1 (MQ=255)                               tAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCa                                                                  >  2:2391114/1‑100 (MQ=255)                                aGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAg                                                                 <  2:732329/100‑1 (MQ=255)                                    gATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGttac                                                             >  2:1611888/1‑100 (MQ=255)                                    gATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGttac                                                             >  1:455173/1‑100 (MQ=255)                                     ttGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGttact                                                            <  2:2402054/99‑1 (MQ=255)                                      tGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGttactt                                                           <  1:323539/100‑1 (MQ=255)                                      tGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGttacta                                                           >  1:2402068/1‑99 (MQ=255)                                        cTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGttacttac                                                         <  1:73015/100‑1 (MQ=255)                                        cTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGttacttac                                                         <  2:156354/100‑1 (MQ=255)                                        cTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGttacttac                                                         >  1:819081/1‑100 (MQ=255)                                         ttttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTtacttact                                                        <  2:2543695/100‑1 (MQ=255)                                         ttttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTtacttact                                                        <  1:1696356/100‑1 (MQ=255)                                         ttttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTtacttact                                                        >  2:660005/1‑100 (MQ=255)                                         ttttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTtacttact                                                        <  1:2352422/100‑1 (MQ=255)                                          tttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCCTAATGTAAATTAACCAGTTacttactt                                                       >  2:1572110/1‑100 (MQ=255)                                           ttCCGGATAGCAATTATCTAAACGTAACCCTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTacttactt                                                       <  2:2342695/99‑1 (MQ=255)                                           ttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACttacttt                                                      <  2:2504666/100‑1 (MQ=255)                                             ccGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACttactttac                                                    <  2:987891/100‑1 (MQ=255)                                              cGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGAACGCGACCAGAATGTAAATTAACCAGTTACttactttact                                                   <  2:757879/100‑1 (MQ=255)                                                gATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTtactac                                                 <  1:1247173/100‑1 (MQ=255)                                                gATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTtactac                                                 <  2:820355/100‑1 (MQ=255)                                                  tAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTGTACTACan                                               <  2:2369462/100‑2 (MQ=255)                                                    gCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATg                                             >  1:2222340/1‑100 (MQ=255)                                                    gCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATg                                             >  1:2283403/1‑100 (MQ=255)                                                    gCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATg                                             <  2:227626/100‑1 (MQ=255)                                                    gCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATg                                             >  2:1874373/1‑100 (MQ=255)                                                     cAATTATCTAAACGTAACACTTGACGACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGt                                            >  1:899411/1‑100 (MQ=255)                                                     cAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGt                                            <  2:1352796/100‑1 (MQ=255)                                                     cAATTATCTAAACGTAACACTTGACAACTGTGTAAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTGACTACAATGt                                            >  2:1682947/1‑100 (MQ=255)                                                      aaTTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTa                                           >  2:2034716/1‑100 (MQ=255)                                                      aaTTATCTAAACGTAACACTTAACAACTTTGTGAGAAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGt                                            <  2:1414761/99‑1 (MQ=255)                                                         tATCTAAACGTAACACTTGACAACTGTGTGATGAAAATCGTTTTTGCGCACGCTACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAAcc                                        <  2:878150/100‑1 (MQ=255)                                                          aTCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATATCAATTACCCACTTACTTACTTTACTACAATGTAACCg                                       >  2:1088131/1‑100 (MQ=255)                                                            cTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGc                                     <  2:1615653/100‑1 (MQ=255)                                                             tAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACTAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAAcgtgca                                    >  2:1438567/1‑95 (MQ=255)                                                             tAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCa                                    >  1:240615/1‑100 (MQ=255)                                                             tAAACGCAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCa                                    <  1:1131929/100‑1 (MQ=255)                                                              aaaCGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAg                                   >  2:368679/1‑100 (MQ=255)                                                              aaaCGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAg                                   <  1:2305838/100‑1 (MQ=255)                                                                 cGTAACACTTGACAACTGTGTGAGGAGAAGCGTTTTTGCGCACGCGAACAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACAGGAAGTGa                                <  2:721452/100‑1 (MQ=255)                                                                  gTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAt                               >  2:953705/1‑100 (MQ=255)                                                                   tAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAtt                              <  1:488838/100‑1 (MQ=255)                                                                   tAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAt                               >  2:176042/1‑99 (MQ=255)                                                                    aaCACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGACAGTGAttt                             >  2:301467/1‑100 (MQ=255)                                                                      cacTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAACTTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAttttt                           <  2:1424081/100‑1 (MQ=255)                                                                       acTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTAACTTACTTTACTACAATGTAACCGGCAGTGAtttttt                          <  2:1366781/100‑1 (MQ=255)                                                                        cTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTa                         >  1:751815/1‑100 (MQ=255)                                                                         ttGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTaa                        <  2:787243/100‑1 (MQ=255)                                                                         ttGACAACTGTGTAATGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTACCCGGTTACTTACTTTACTACAATGTAACCGGCAGTGAttttt                           <  2:1065504/97‑1 (MQ=255)                                                                           taCAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAag                       >  2:508457/2‑99 (MQ=255)                                                                           gACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAaga                      >  2:1599520/1‑100 (MQ=255)                                                                           gACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAaga                      >  1:2000354/1‑100 (MQ=255)                                                                              cttccgatctggAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                   <  1:638074/90‑1 (MQ=255)                                                                                cTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCAt                 >  2:343111/1‑100 (MQ=255)                                                                                    gtgAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCt             >  2:73293/1‑100 (MQ=255)                                                                                    gtgAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCt             >  1:501432/1‑100 (MQ=255)                                                                                      gAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTc           >  1:1690856/1‑100 (MQ=255)                                                                                      gAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTc           <  1:1273304/100‑1 (MQ=255)                                                                                       aGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTca          <  2:777171/100‑1 (MQ=255)                                                                                       aGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTca          <  1:918086/100‑1 (MQ=255)                                                                                        ggAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCagatcggaag         >  2:638067/1‑91 (MQ=255)                                                                                        ggAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag         >  1:1259292/1‑100 (MQ=255)                                                                                        ggAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag         >  2:1426119/1‑100 (MQ=255)                                                                                        ggAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag         <  2:2484994/100‑1 (MQ=255)                                                                                        ggAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag         >  1:2093884/1‑100 (MQ=255)                                                                                         gAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagc        >  2:1756494/1‑100 (MQ=255)                                                                                         gAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagc        <  1:2323008/100‑1 (MQ=255)                                                                                          aaaaGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagca       <  2:501426/100‑1 (MQ=255)                                                                                           aaaGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagcag      <  1:2438838/100‑1 (MQ=255)                                                                                           aaaGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagcag      >  1:69732/1‑100 (MQ=255)                                                                                             aGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGaa    >  1:1871115/1‑100 (MQ=255)                                                                                              gCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGaaa   <  1:1652704/100‑1 (MQ=255)                                                                                               cGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCATTCAGCAGAAAt  >  1:1263813/1‑100 (MQ=255)                                                                                                  |                                                                                                 TGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAT  >  NC_000913/160639‑160832

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