Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,096,938	0	G	T	25.6%	174.9 / 59.2	82	A86A (GCG→GCT)	rdgB	dITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base T (9/12);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

AGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGGCGCGTTTCC  >  NC_000913/3096839‑3097035                                                                                                    |                                                                                                  aGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgcg                                                                                                   <  1:402257/100‑1 (MQ=255)   cgcgcCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgcgcg                                                                                                 >  1:509050/1‑100 (MQ=255)    gcgcCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGt                                                                                                <  1:1613903/100‑1 (MQ=255)     cgcCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGtt                                                                                               >  1:502098/1‑100 (MQ=255)         aTGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTc                                                                                           <  2:1354201/100‑1 (MQ=255)         aTGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTc                                                                                           <  1:195290/100‑1 (MQ=255)           gCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTc                                                                                           <  1:593956/98‑1 (MQ=255)           gCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCg                                                                                         >  2:632076/1‑100 (MQ=255)            cGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCgg                                                                                        <  2:435500/100‑1 (MQ=255)            cGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCgg                                                                                        <  1:1289831/100‑1 (MQ=255)             ggCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGt                                                                                       >  2:819593/1‑100 (MQ=255)                aaaaGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGaa                                                                                    <  2:783250/100‑1 (MQ=255)                     tGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgc                                                                               <  2:1676746/100‑1 (MQ=255)                      gACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgcg                                                                              <  1:548828/100‑1 (MQ=255)                      gACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgcg                                                                              <  1:548750/100‑1 (MQ=255)                       aCCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGa                                                                             >  2:117462/1‑100 (MQ=255)                        ccGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAc                                                                            >  1:1565592/1‑100 (MQ=255)                         cGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAcc                                                                           >  1:52171/1‑100 (MQ=255)                         cGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAcc                                                                           <  2:949449/100‑1 (MQ=255)                             ttACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATc                                                                       >  2:102585/1‑100 (MQ=255)                             ttACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATc                                                                       <  1:746075/100‑1 (MQ=255)                             ttACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATc                                                                       <  1:1649025/100‑1 (MQ=255)                             ttACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATc                                                                       <  2:20415/100‑1 (MQ=255)                                ccGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCaaa                                                                    >  2:551842/1‑100 (MQ=255)                                ccGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaa                                                                    <  2:757082/100‑1 (MQ=255)                                ccGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaa                                                                    <  2:204338/100‑1 (MQ=255)                                ccGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaa                                                                    <  1:445117/100‑1 (MQ=255)                                ccGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaa                                                                    >  2:1383686/1‑100 (MQ=255)                                 cGGCAATTGCCGACGGCTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCaaa                                                                    <  1:1896371/99‑1 (MQ=255)                                 cGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaan                                                                   <  1:176830/100‑2 (MQ=255)                                 cGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaaa                                                                   >  1:538029/1‑100 (MQ=255)                                 cGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaaa                                                                   <  2:873916/100‑1 (MQ=255)                                  ggCAATTGTCGTCGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaaat                                                                  <  1:459184/100‑2 (MQ=255)                                  ggCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAg                                                                  >  1:418572/1‑100 (MQ=255)                                  ggCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAg                                                                  <  2:1773778/100‑1 (MQ=255)                                  ggCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCACTCGTTATTCCGGTGAAGACGCGACCGATCaaaa                                                                   <  1:1886458/99‑1 (MQ=255)                                     aaTTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAAt                                                               >  2:1272886/1‑100 (MQ=255)                                         gCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGc                                                           <  2:234710/100‑1 (MQ=255)                                          ccGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTTCGGTGAAGACGCGACCGATCAAAAGAATCTGCa                                                          >  1:1511708/1‑100 (MQ=255)                                            gacgacTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaaa                                                        >  1:171784/1‑100 (MQ=255)                                            gacgacTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaan                                                        <  2:1641310/100‑2 (MQ=255)                                               gacTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaaaaa                                                      <  1:807176/99‑1 (MQ=255)                                                      gtctggCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGa                                              >  2:1163332/1‑100 (MQ=255)                                                       tctggCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGaa                                             >  1:1932360/1‑100 (MQ=255)                                                       tctggCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGaa                                             <  1:415132/100‑1 (MQ=255)                                                        ctggCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGaaa                                            >  2:1840460/1‑100 (MQ=255)                                                          ggCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGCTCAAAAGAATCTGCAAAAACTGCTGGAAACa                                          >  1:590555/1‑100 (MQ=255)                                                          ggCGGTACATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACa                                          >  2:1019709/1‑100 (MQ=255)                                                           gCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACaa                                         >  2:81708/1‑100 (MQ=255)                                                            cGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAAt                                        <  2:1864180/100‑1 (MQ=255)                                                             ggTAGTTGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATg                                       >  1:230969/1‑100 (MQ=255)                                                             ggTAGATGTGCTTGGCGGCGCGGCGGGCATTTACTCCGCGCGTTAGTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCCAGAACTGCTGGAGACAATg                                       >  1:1131562/1‑100 (MQ=255)                                                             ggTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATg                                       <  2:1565603/100‑1 (MQ=255)                                                             ggTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATg                                       <  2:502100/100‑1 (MQ=255)                                                             ggTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATg                                       >  1:932614/1‑100 (MQ=255)                                                                aGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGaaa                                    <  2:1548377/100‑1 (MQ=255)                                                                 gATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAg                                   <  2:1574882/100‑1 (MQ=255)                                                                 gATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAg                                   >  2:1644182/1‑100 (MQ=255)                                                                  aTGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGa                                  <  2:1235854/100‑1 (MQ=255)                                                                    gtgCTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACg                                >  1:336322/1‑100 (MQ=255)                                                                    gtgCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACg                                <  1:1329544/100‑1 (MQ=255)                                                                      gCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTa                              >  2:20833/1‑100 (MQ=255)                                                                      gCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTa                              <  2:489136/100‑1 (MQ=255)                                                                       cTTGGCGGCGCGCCGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACAGATCAAAAGAATCTGCAAAAACTGCTGGACACAATGAAAGAAGTa                              >  2:1622656/1‑99 (MQ=255)                                                                       cTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTAc                             <  2:76235/100‑1 (MQ=255)                                                                       cTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTAc                             >  2:1399375/1‑100 (MQ=255)                                                                          ggcggcGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCgg                          <  1:117462/100‑1 (MQ=255)                                                                               cgcgcCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGgacgac                     >  2:1799891/1‑100 (MQ=255)                                                                               cgcgcCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGgacgac                     >  2:446409/1‑100 (MQ=255)                                                                                gcgcCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATATGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGAcc                    <  1:385411/100‑1 (MQ=255)                                                                                   ccGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAAc                 <  1:108123/100‑1 (MQ=255)                                                                                   ccGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAAc                 <  2:874850/100‑1 (MQ=255)                                                                                    cGGGCATTTACTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACg                <  1:1848793/100‑1 (MQ=255)                                                                                     gggCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGt               >  1:1372464/1‑100 (MQ=255)                                                                                        cATTTACACCGCGCGTTATTCCGGTGAAGACGCTACCGATCAAAAGAATCTGCAAAAACTTCTGGAAACAATGAAAGACGTACCGGACGACTAACGTCa             >  2:1021730/1‑99 (MQ=255)                                                                                          tttACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGgc          <  2:1064330/100‑1 (MQ=255)                                                                                            tACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGgcgc        <  2:1923071/100‑1 (MQ=255)                                                                                            tACTCCGAGCGTTATTCCGGTGAAGACGCGACCGGCCAAAAGAATCTGCAACAAATGCTGGAAACAACGAAAGACGTACCAGACGACCCACGTCAgg           >  1:639880/1‑97 (MQ=255)                                                                                             aCTCCGCTCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGGcgcg       <  2:171785/100‑1 (MQ=255)                                                                                             aCTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGGcgcg       <  2:1919071/100‑1 (MQ=255)                                                                                              cTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACGATGAAAGACGTACCGGCCGACCAACGTCAGGCGCGt      >  2:1365830/1‑100 (MQ=255)                                                                                              cTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGGCGCGt      >  1:746279/1‑100 (MQ=255)                                                                                              cTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGGCGCGt      >  1:868022/1‑100 (MQ=255)                                                                                               tCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGGCGCGtt     <  2:828839/100‑1 (MQ=255)                                                                                                 cgcgcgTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGGCGCGTTTc   >  2:173029/1‑100 (MQ=255)                                                                                                  gcgcgTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGGCGCGTTTcc  >  1:560603/1‑100 (MQ=255)                                                                                                    |                                                                                                  AGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGGCGCGTTTCC  >  NC_000913/3096839‑3097035

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