Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,964,373	C→T	6.0%	P261P (CCG→CCA)	thyA ←	thymidylate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,964,373	0	C	T	6.0%	222.2 / 3.7	67	P261P (CCG→CCA)	thyA	thymidylate synthetase
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base T (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

TTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATAT  >  NC_000913/2964225‑2964521                                                                                                                                                     |                                                                                                                                                     tttCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCAcc                                                                                                                                                      >  2:1102581/1‑149 (MQ=255)      tcGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCggtgc                                                                                                                                                 >  3:62477/1‑146 (MQ=255)       cGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCt                                                                                                                                                <  3:200388/149‑1 (MQ=255)              tCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGc                                                                                                                                         <  4:806501/149‑1 (MQ=255)                ttCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGccc                                                                                                                                       >  2:884274/1‑149 (MQ=255)                   cGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGa                                                                                                                                    >  3:281835/1‑149 (MQ=255)                    ggACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGAt                                                                                                                                   <  4:1119934/149‑1 (MQ=255)                     gACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGTGCTTTAATGCCCGGATg                                                                                                                                  <  4:62477/149‑1 (MQ=255)                                ccGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGc                                                                                                                       <  3:950435/149‑1 (MQ=255)                                  gTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCCCCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCt                                                                                                                     <  1:233362/149‑1 (MQ=255)                                     gCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCa                                                                                                                  >  1:787503/1‑149 (MQ=255)                                       aGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAAt                                                                                                                >  1:1036127/1‑149 (MQ=255)                                               gCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACTGGCGCTTTAATGCCCGGATGCgg                                                                                                                               <  2:814329/126‑1 (MQ=255)                                               gCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACTGGCGCTTTAATGCCCGGATGCgg                                                                                                                               >  1:814329/1‑126 (MQ=255)                                               gCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGt                                                                                                        >  2:904453/1‑149 (MQ=255)                                               gCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGt                                                                                                        >  1:392347/1‑149 (MQ=255)                                                 aaaaTTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCt                                                                                                      >  1:705751/1‑149 (MQ=255)                                                  aaaTTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCtt                                                                                                     <  4:277405/149‑1 (MQ=255)                                                  aaaTTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCtt                                                                                                     <  4:487513/149‑1 (MQ=255)                                                     tttCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGa                                                                                                  >  2:861513/1‑149 (MQ=255)                                                       tCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGaaa                                                                                                <  2:787503/149‑1 (MQ=255)                                                         gggAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACg                                                                                              >  1:1055028/1‑149 (MQ=255)                                                           gAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGt                                                                                            >  1:364325/1‑149 (MQ=255)                                                            aaGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTa                                                                                           <  1:1102581/149‑1 (MQ=255)                                                                  tctcGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACTGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGaa                                                                                     >  4:322157/1‑149 (MQ=255)                                                                  tctcGAAGAATTGAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGaa                                                                                     >  3:84465/1‑149 (MQ=255)                                                                    tcGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGa                                                                                   <  2:1055028/149‑1 (MQ=255)                                                                     cGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGAt                                                                                  <  2:364325/149‑1 (MQ=255)                                                                      gaagaaTTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATg                                                                                 >  2:919620/1‑149 (MQ=255)                                                                        agaaTTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGa                                                                               <  4:84465/149‑1 (MQ=255)                                                                         gaaTTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGAt                                                                              <  4:281835/149‑1 (MQ=255)                                                                           aTTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTc                                                                            <  2:1036127/149‑1 (MQ=255)                                                                            tttAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCg                                                                           >  2:486677/1‑149 (MQ=255)                                                                            tttAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCg                                                                           >  3:230620/1‑149 (MQ=255)                                                                             ttAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCgg                                                                          <  1:904453/149‑1 (MQ=255)                                                                                      ggTAAAAAAACCGACGTACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGtt                                                                 >  3:523913/1‑149 (MQ=255)                                                                                         aaaaaaaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGa                                                              <  1:884274/149‑1 (MQ=255)                                                                                          aaaaaaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGAt                                                             >  3:136740/1‑149 (MQ=255)                                                                                           aaaaaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGc                                                                                                                                 >  3:520194/1‑80 (MQ=255)                                                                                           aaaaaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGc                                                                                                                                 <  4:520194/80‑1 (MQ=255)                                                                                              aaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATc                                                                                                                            >  2:926528/1‑82 (MQ=255)                                                                                              aaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATc                                                                                                                            <  1:926528/82‑1 (MQ=255)                                                                                              aaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATc                                                         >  4:810184/1‑149 (MQ=255)                                                                                                 cGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAAc                                                      >  1:820142/1‑149 (MQ=255)                                                                                                   aCGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACtt                                                    >  1:963169/1‑149 (MQ=255)                                                                                                     gCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCg                                                  <  1:486677/149‑1 (MQ=255)                                                                                                       acacTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGc                                                >  3:1146157/1‑149 (MQ=255)                                                                                                       acacTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGc                                                >  1:880353/1‑149 (MQ=255)                                                                                                          cTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTTAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGc                                             <  2:820142/149‑1 (MQ=255)                                                                                                                  ggCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCgg                                     <  4:1137216/149‑1 (MQ=255)                                                                                                                  ggCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCgg                                     >  1:242610/1‑149 (MQ=255)                                                                                                                   gCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGt                                    >  2:36131/1‑149 (MQ=255)                                                                                                                    ctctGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGtt                                   >  1:798459/1‑149 (MQ=255)                                                                                                                       tGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTcgc                                >  2:716005/1‑149 (MQ=255)                                                                                                                         gCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCgg                              >  1:253608/1‑149 (MQ=255)                                                                                                                             gATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTa                          >  3:887305/1‑149 (MQ=255)                                                                                                                                gTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAAtt                       <  4:1146157/149‑1 (MQ=255)                                                                                                                                   tCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCTGACGCGGTTCGCGGCTTAATTGCa                    >  2:750432/1‑149 (MQ=255)                                                                                                                                   tCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCa                    <  4:136740/149‑1 (MQ=255)                                                                                                                                      tAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTACAGAt                 <  2:392347/149‑1 (MQ=255)                                                                                                                                        aTTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGa               <  1:919620/149‑1 (MQ=255)                                                                                                                                        aTTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGCCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGa               <  2:290313/149‑1 (MQ=255)                                                                                                                                            gATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTT‑GAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTg          >  3:1066415/1‑149 (MQ=255)                                                                                                                                            gATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGCCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGGTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGttt           >  1:655549/1‑149 (MQ=255)                                                                                                                                               aGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGAt        <  2:880353/149‑1 (MQ=255)                                                                                                                                                gCCACTGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATc       <  3:322157/149‑1 (MQ=255)                                                                                                                                                gCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATc       <  1:861513/149‑1 (MQ=255)                                                                                                                                                   aCCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGTGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCat    >  3:704831/1‑149 (MQ=255)                                                                                                                                                     cGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGTGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCatat  <  4:704831/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATAT  >  NC_000913/2964225‑2964521

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