Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,172,619	C→T	9.4%	G30G (GGC→GGT)	mdaB →	NADPH quinone reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,172,619	0	C	T	9.4%	215.9 / 10.2	74	G30G (GGC→GGT)	mdaB	NADPH quinone reductase
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base T (4/3);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

ACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAA  >  NC_000913/3172473‑3172766                                                                                                                                                   |                                                                                                                                                    aCCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGca                                                                                                                                                    >  1:397877/1‑148 (MQ=255)        cGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACCCCCTGACCGAAGTCGCGGATGGCACACTgcgc                                                                                                                                            >  1:451001/1‑149 (MQ=255)          ctctGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGa                                                                                                                                          >  4:244958/1‑149 (MQ=255)               tgACGTACTATAGTAGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCtt                                                                                                                                      <  4:100003/148‑1 (MQ=255)                      cTATAGTCGGCAAACCTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCatg                                                                                                                               <  7:410741/148‑1 (MQ=255)                         taGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGATAAAATTCGCACACTCCCATGGTCAACTGAACGACACCCTGCCCGAAGTCCCGGATGGCACACTGCGCGACCATGGGCATGAAGt                                                                                                                           >  6:190852/1‑149 (MQ=255)                         taGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCatgat                                                                                                                             <  4:299314/147‑1 (MQ=255)                           gTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGt                                                                                                                           >  2:361413/1‑147 (MQ=255)                           gTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGt                                                                                                                           <  1:361413/147‑1 (MQ=255)                            tCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCg                                                                                                                        <  2:451001/149‑1 (MQ=255)                            tCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCg                                                                                                                        >  3:13411/1‑149 (MQ=255)                              ggCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGc                                                                                                                       <  7:222379/148‑1 (MQ=255)                                cAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATc                                                                                                                    >  6:306117/1‑149 (MQ=255)                                 aaaCGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGTACACTGCGCGACCTTGGGCATGATGTCCGCATCg                                                                                                                   >  2:19128/1‑149 (MQ=255)                                 aaaCGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGTACACTGCGCGACCTTGGGCATGATGTCCGCATCg                                                                                                                   >  2:26211/1‑149 (MQ=255)                                  aaCGTCTCACCTTGAGGTTAAAAAATGAGCAACCTCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTgg                                                                                                                                    <  5:307021/131‑1 (MQ=255)                                  aaCGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTgg                                                                                                                                    >  6:307021/1‑131 (MQ=255)                                       ctcACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTcgcg                                                                                                             >  7:188278/1‑149 (MQ=255)                                                   ttAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGAc                                                                                                  >  1:443697/1‑148 (MQ=255)                                                   ttAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGAc                                                                                                  >  2:326527/1‑148 (MQ=255)                                                     aaaaaaTGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTAc                                                                                               >  1:290857/1‑149 (MQ=255)                                                      aaaaaTGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTAc                                                                                               <  4:13411/148‑1 (MQ=255)                                                        aaaTGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGAt                                                                                            <  5:448945/149‑1 (MQ=255)                                                        aaaTGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGAt                                                                                            <  2:397877/149‑1 (MQ=255)                                                           tGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACACTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGAAAGCGACTACGATGt                                                                                          >  8:99779/1‑148 (MQ=255)                                                               cAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAg                                                                                     >  2:405098/1‑149 (MQ=255)                                                                    tCCTGATTATCACCGGCGCGAAAAACTTCTCCCACTCCAATGGTCACCTGAAAGCCACCCTGACCGAAGTCGCGGATGGCACCCTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCTGaa                                                                                <  5:31855/149‑1 (MQ=255)                                                                    tCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGa                                                                                             <  6:443701/136‑1 (MQ=255)                                                                    tCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGa                                                                                             >  5:443701/1‑136 (MQ=255)                                                                    tCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGaa                                                                                <  2:290857/149‑1 (MQ=255)                                                                      cTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAg                                                                               >  1:66163/1‑148 (MQ=255)                                                                         aTTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACTCCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTAc                                                                            <  5:190852/148‑1 (MQ=255)                                                                         aTTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGTACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTa                                                                             <  1:26211/147‑1 (MQ=255)                                                                         aTTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGTACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTa                                                                             <  1:19128/147‑1 (MQ=255)                                                                         aTTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTa                                                                             <  6:413508/147‑1 (MQ=255)                                                                         aTTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACa                                                                           <  3:244958/149‑1 (MQ=255)                                                                               aaCGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACt                                                                      <  8:188278/148‑1 (MQ=255)                                                                                    cgcgAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTctct                                                                <  5:306117/149‑1 (MQ=255)                                                                                      cgAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTgg                                                              >  8:6761/1‑149 (MQ=255)                                                                                       gAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGca                                                                                                                                                    >  4:322886/1‑62 (MQ=255)                                                                                       gAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGca                                                                                                                                                    <  3:322886/62‑1 (MQ=255)                                                                                       gAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTggg                                                             <  2:66163/149‑1 (MQ=255)                                                                                                    cACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGa                                                 >  8:174943/1‑148 (MQ=255)                                                                                                       tCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCgg                                                                                  <  5:291545/112‑1 (MQ=255)                                                                                                       tCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCgg                                                                                  >  6:291545/1‑112 (MQ=255)                                                                                                         cAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTgg                                            >  3:288357/1‑148 (MQ=255)                                                                                                          aaTGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCa                                          <  1:326527/149‑1 (MQ=255)                                                                                                             ggTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGAt                                       >  2:291963/1‑149 (MQ=255)                                                                                                                 aaCTGAACGACACCCTGACCGAAGTCGCGGATGGTACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGc                                     >  5:273634/1‑147 (MQ=255)                                                                                                                    tGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGcc                                    >  2:199197/1‑145 (MQ=255)                                                                                                                    tGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGcc                                    <  1:199197/145‑1 (MQ=255)                                                                                                                      aaCGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAACCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAgacgt                              >  6:292900/1‑145 (MQ=255)                                                                                                                      aaCGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGc                                <  1:291963/147‑1 (MQ=255)                                                                                                                      aaCGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCtg                              <  7:174943/149‑1 (MQ=255)                                                                                                                            acCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTctct                                                                >  7:426726/1‑109 (MQ=255)                                                                                                                            acCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTctct                                                                <  8:426726/109‑1 (MQ=255)                                                                                                                              ccTGACCGAAGTCGCGGATGGTACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGCCAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATgg                      >  1:470513/1‑149 (MQ=255)                                                                                                                              ccTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTAGGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATgg                      <  1:405098/149‑1 (MQ=255)                                                                                                                              ccTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGGAGTACAAAACTTTCTCTGGGc                                                            >  5:172491/1‑111 (MQ=255)                                                                                                                              ccTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGGAGTACAAAACTTTCTCTGGGc                                                            <  6:172491/111‑1 (MQ=255)                                                                                                                              ccTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATgtggtg                                                  >  3:182646/1‑121 (MQ=255)                                                                                                                              ccTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATgtggtg                                                  <  4:182646/121‑1 (MQ=255)                                                                                                                              ccTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATgg                      >  4:130889/1‑149 (MQ=255)                                                                                                                              ccTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATg                       >  8:382845/1‑148 (MQ=255)                                                                                                                              ccTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATg                       <  7:382845/148‑1 (MQ=255)                                                                                                                              ccTGACCGAAATCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATgg                      >  3:277808/1‑149 (MQ=255)                                                                                                                               atGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTTGTGATCTGGCAGATGCCAGGCTGGTGGATggg                     >  8:143770/2‑149 (MQ=255)                                                                                                                                  aCCGAAATCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGgc                    <  4:277808/147‑1 (MQ=255)                                                                                                                                    cGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATg                       >  7:407867/1‑142 (MQ=255)                                                                                                                                    cGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATg                       <  8:407867/142‑1 (MQ=255)                                                                                                                                    cGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGgcgcgc                >  1:416893/1‑149 (MQ=255)                                                                                                                                      aaGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAg                                  <  7:52046/129‑1 (MQ=255)                                                                                                                                      aaGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAg                                  >  8:52046/1‑129 (MQ=255)                                                                                                                                         tCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGgc                    >  8:126896/1‑140 (MQ=255)                                                                                                                                         tCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGgc                    <  7:126896/140‑1 (MQ=255)                                                                                                                                           gcgGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGa          <  2:443697/148‑1 (MQ=255)                                                                                                                                              gATGGTACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGt      <  6:273634/149‑1 (MQ=255)                                                                                                                                                   cacaCTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGaaa  >  6:227006/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAACGACACCCTGACCGAAGTCGCGGATGGCACACTGCGCGACCTTGGGCATGATGTCCGCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAA  >  NC_000913/3172473‑3172766

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: ‑