Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,211,722	T→C	6.1%	F206L (TTC→CTC)	dnaG →	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,211,722	0	T	C	6.1%	194.3 / 3.8	65	F206L (TTC→CTC)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (2/2);  total (33/32)
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.97e-01

TTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGAT  >  NC_000913/3211579‑3211868                                                                                                                                                |                                                                                                                                                   ttGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTccc                                                                                                                                                >  5:129857/1‑148 (MQ=255)  tGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTccc                                                                                                                                                >  2:178908/1‑147 (MQ=255)   ggTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTcccc                                                                                                                                               <  6:129857/147‑1 (MQ=255)         gcgcCCCCCGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTcgc                                                                                                                                         >  2:328308/1‑147 (MQ=255)             cccccGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATaa                                                                                                                                    <  2:169060/148‑1 (MQ=255)             cccccGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATaa                                                                                                                                    <  3:172527/148‑1 (MQ=255)                ccGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTcgcg                                                                                                                                        <  1:298630/141‑1 (MQ=255)                ccGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTcgcg                                                                                                                                        >  2:298630/1‑141 (MQ=255)                 cGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAAc                                                                                                                                  >  6:347207/1‑146 (MQ=255)                       ggACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCgg                                                                                                                          <  2:251263/148‑1 (MQ=255)                        gACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCggg                                                                                                                         <  1:178908/148‑1 (MQ=255)                                ccTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGt                                                                                                                 >  3:189445/1‑148 (MQ=255)                                 ttGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTAATTGATTGATGCGGGGATGTTGGTCACTAACTATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTTATTGGtt                                                                                                                <  5:123744/147‑1 (MQ=255)                                    aaGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTg                                                                                                             <  1:328308/148‑1 (MQ=255)                                       cGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTggcg                                                                                                          <  8:334690/148‑1 (MQ=255)                                        ggTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGt                                                                                                                                                    >  4:119682/1‑105 (MQ=255)                                        ggTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGt                                                                                                                                                    <  3:119682/105‑1 (MQ=255)                                              gcggcAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGt                                                                                                   >  1:83017/1‑148 (MQ=255)                                              gcggcAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGt                                                                                                   <  5:347207/148‑1 (MQ=255)                                                     tCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTggg                                                                                             >  1:158563/1‑147 (MQ=255)                                                      ccAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGc                                                                                            >  2:385556/1‑147 (MQ=255)                                                      ccAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGc                                                                                            >  1:5260/1‑147 (MQ=255)                                                         gAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACg                                                                                        <  7:387086/148‑1 (MQ=255)                                                           aaaTCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGcgcg                                                                                                    >  2:191195/1‑134 (MQ=255)                                                           aaaTCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGa                                                                                       >  2:293914/1‑147 (MQ=255)                                                           aaaTCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGAt                                                                                      <  4:189445/148‑1 (MQ=255)                                                           aaaTCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCAATCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGcgcg                                                                                                    <  1:191195/134‑1 (MQ=255)                                                            aaTCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATa                                                                                     <  2:5260/148‑1 (MQ=255)                                                                      cATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATaa                                                                                                                                    <  1:26605/91‑1 (MQ=255)                                                                      cATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATaa                                                                                                                                    >  2:26605/1‑91 (MQ=255)                                                                           attgatGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGa                                                                      <  1:293914/148‑1 (MQ=255)                                                                                  cGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGc                                                                <  2:158563/147‑1 (MQ=255)                                                                                     ggATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGTGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGttt                                                                                                               <  6:329464/97‑1 (MQ=255)                                                                                     ggATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGTGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGttt                                                                                                               >  5:329464/1‑97 (MQ=255)                                                                                               cACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACAtt                                                  >  7:395154/1‑148 (MQ=255)                                                                                                 cTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCaa                                                                              <  7:110699/118‑1 (MQ=255)                                                                                                 cTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCaa                                                                              >  8:110699/1‑118 (MQ=255)                                                                                                   aaCGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGTCAACGATACCCCCAAATACCTGACCTCGCCGGAAACAGACATTTTc                                               >  6:361277/1‑147 (MQ=255)                                                                                                   aaCGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAACCAGACATTTTc                                               >  7:237952/1‑147 (MQ=255)                                                                                                   aaCGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACAt                                                   <  7:387219/143‑1 (MQ=255)                                                                                                   aaCGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACAt                                                   >  8:387219/1‑143 (MQ=255)                                                                                                      gATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATa                                           <  1:385556/148‑1 (MQ=255)                                                                                                       aTCAGGGACGCACTGACGCTCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGAGAAACGCGGTCGGGTGAGTGGCTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACCGACACTTTCCaaaa                                          <  8:237952/148‑1 (MQ=255)                                                                                                            ggACGAAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATACGAGATAAACGCGTTCGGGTGAAAGGTTTTTGCAGGCGCGTCCTGGGCAACGATACCCCCAACTACCTGAACTCGCCGGAAACACACATTTTCCATAAAGgcc                                     >  2:71676/1‑148 (MQ=255)                                                                                                                 cAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAg                                >  6:94349/1‑148 (MQ=255)                                                                                                                   gTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCt                              >  7:74681/1‑148 (MQ=255)                                                                                                                    ttACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCtt                             >  3:106844/1‑148 (MQ=255)                                                                                                                    ttACGATCGTTTCAGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCaa                                                                              <  6:278745/99‑1 (MQ=255)                                                                                                                     tACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCtt                             >  7:119147/1‑147 (MQ=255)                                                                                                                     tACGATCGTTTCCGCGAGCGGCTCATGTTCCCCATACGCGATAACCGCGGTCGGGTCATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAACTACCTGAACTCGCCGGAAACAGACATTATCCATAAAGGCCGCCAGCtt                             >  5:34227/1‑147 (MQ=255)                                                                                                                      aCGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTa                           >  1:379516/1‑148 (MQ=255)                                                                                                                      aCGATCGATTCCGCGAGCGGGTGATGTTACCCATTCGCGATTAACGCGGTCGGGAGATTGCTTTTGGCGGGCGCGTGCTGGGCAACGAGACACACAAATACCTGATCTAGCCGGATGCAGACAGTTTCGATAAAGTCCGCGAGCtttt                           >  1:91260/1‑147 (MQ=255)                                                                                                                           cGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACAttt                                                 >  2:37733/1‑121 (MQ=255)                                                                                                                           cGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACAttt                                                 <  1:37733/121‑1 (MQ=255)                                                                                                                           cGTTTCCGCGAGCGGGTGATGTTCCCCATTAGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGAGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTc                      >  1:391997/1‑148 (MQ=255)                                                                                                                           cGTTTCCGCGAGCGGGTGATGCTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTc                      >  7:298380/1‑148 (MQ=255)                                                                                                                           cGTTTCCGCGAGCGGGTGATGCTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTc                      >  7:303668/1‑148 (MQ=255)                                                                                                                              ttCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCtt                    >  3:190433/1‑147 (MQ=255)                                                                                                                               tCCGCGAGCGGGTGATGTTCCCCATTCTCGATAAACGCGGTCGGGTTATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAATTACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCttt                   <  2:391997/147‑1 (MQ=255)                                                                                                                               tCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCttt                   <  4:106844/147‑1 (MQ=255)                                                                                                                               tCCGCGAGCGGGTGATGCTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCttt                   <  8:298380/147‑1 (MQ=255)                                                                                                                               tCCGCGAGCGGGTGATGCTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCttt                   <  8:303668/147‑1 (MQ=255)                                                                                                                                  gcgAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGa               <  5:94349/148‑1 (MQ=255)                                                                                                                                      gCGGTTGATGTTCCCCATACGCGATAAACGCGGTCTGGTGATTGGCTATGGCGGGCGCGTGCTGGGCAACGTTACCCCCAAATAGCTGAACTCGCCGGAAAGAGACATATTCCATAAACGCCGCCAGCTTTACGGTCTTTATGAAgcg           <  2:91260/148‑1 (MQ=255)                                                                                                                                       cGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAgcgc          <  8:395154/148‑1 (MQ=255)                                                                                                                                            gATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGcagcag     <  4:190433/148‑1 (MQ=255)                                                                                                                                               gTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGa   >  5:104486/1‑147 (MQ=255)                                                                                                                                                ttCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGAt  >  6:226158/1‑147 (MQ=255)                                                                                                                                                |                                                                                                                                                   TTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCGGTTTGGCGGCAATCCAGAAAATCGCCAGTCATTGATTGATGCGGGGATGTTGGTCACTAACGATCAGGGACGCAGTTACGATCGTTTCCGCGAGCGGGTGATGTTCCCCATTCGCGATAAACGCGGTCGGGTGATTGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGAT  >  NC_000913/3211579‑3211868

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