Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,486,160	G→A	5.6%	G1034S (GGC→AGC)	hrpA →	putative ATP‑dependent helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,486,160	0	G	A	5.6%	234.1 / 2.0	71	G1034S (GGC→AGC)	hrpA	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base A (2/2);  total (35/36)
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.82e-01

CGCTCGCTACAAGATCTGAAAGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTCTACTGCTGCTGAATATTCCATC  >  NC_000913/1486019‑1486299                                                                                                                                              |                                                                                                                                            cgctcgctACAAGATCTGAAAGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACt                                                                                                                                       <  3:52887/148‑1 (MQ=255)     cgctACAAGATCTGAAAGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACa                                                                                                                                    >  7:212825/1‑147 (MQ=255)       ctACAAGATCTGAAAGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAg                                                                                                                                 <  8:282317/148‑1 (MQ=255)       ctACAAGATCTGAAAGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAg                                                                                                                                 <  7:75581/148‑1 (MQ=255)       ctACAAGATCTGAAAGACGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTAGCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAg                                                                                                                                 <  2:237004/148‑1 (MQ=255)          cAAGATCTGAAAGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTg                                                                                                                               <  7:95383/147‑1 (MQ=255)                   aaaGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGc                                                                                                                          <  6:75527/143‑1 (MQ=255)                   aaaGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCg                                                                                                                         <  5:69114/144‑1 (MQ=255)                   aaaGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGc                                                                                                                     >  3:274232/1‑148 (MQ=255)                       aTGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCTTTGGCGGATGAAGGTATAGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGgc                                                                                                                 <  6:21609/148‑1 (MQ=255)                            cTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCtgg                                                                                                            >  7:257074/1‑148 (MQ=255)                               accggCAAACTCCAGCAAACCCCCCCTGCGGCGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCtggtgg                                                                                                         <  2:141579/145‑1 (MQ=255)                                    cAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAg                                                                                                    >  6:49004/1‑148 (MQ=255)                                    cAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAg                                                                                                    <  4:274232/148‑1 (MQ=255)                                    cAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGGTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAg                                                                                                    <  3:78247/148‑1 (MQ=255)                                        gTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAgcgc                                                                                                 >  5:36620/1‑147 (MQ=255)                                         tGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCtggtgg                                                                                                         <  5:153799/138‑1 (MQ=255)                                         tGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCtggtgg                                                                                                         >  6:153799/1‑138 (MQ=255)                                             ggAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGa                                                                                                                              >  7:48378/1‑113 (MQ=255)                                             ggAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGa                                                                                                                              <  8:48378/113‑1 (MQ=255)                                             ggAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCAGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAgt                                                                                           >  4:144343/1‑148 (MQ=255)                                                aaCGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGtgtg                                                                                        <  7:323476/148‑1 (MQ=255)                                                 aCGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACaa                                                                                                                                   <  4:212447/104‑1 (MQ=255)                                                 aCGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACaa                                                                                                                                   >  3:212447/1‑104 (MQ=255)                                                 aCGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTgg                                                                                       >  6:46122/1‑148 (MQ=255)                                                           cGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCaa                                                                                <  5:49004/145‑1 (MQ=255)                                                             gTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGGTGCCGGAAAGCTACGAACAGAAGCGTGGCATCTAGAAAGTGAAGGCGTGGCAGGCGCTGGTGGATGAGCGCGACAGGGTGGCGATCAAACTGt                                                                           <  3:158382/148‑1 (MQ=255)                                                              tGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGt                                                                           >  4:154234/1‑147 (MQ=255)                                                              tGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGt                                                                           <  3:154234/147‑1 (MQ=255)                                                              tGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCCCGACAGTGTGGCGATCAAACTggt                                                                          >  4:158382/1‑146 (MQ=255)                                                                 cGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGa                                                                       <  5:46122/148‑1 (MQ=255)                                                                    aTGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGATCAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATaa                                                                    <  8:257074/148‑1 (MQ=255)                                                                    aTGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATaa                                                                    <  2:326275/148‑1 (MQ=255)                                                                    aTGACGGTATCGAGCAGAGCGGCTT‑CATATCTGGAGGTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAAc                                                                   <  6:106560/148‑1 (MQ=255)                                                                              cGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGAc                                                                                              <  4:92128/112‑1 (MQ=255)                                                                              cGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGAc                                                                                              >  3:92128/1‑112 (MQ=255)                                                                                aGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGc                                                         >  7:12947/1‑147 (MQ=255)                                                                                 gCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGCACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCa                                                        >  1:326736/1‑147 (MQ=255)                                                                                       cGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCa                                                  <  8:212825/147‑1 (MQ=255)                                                                                          cTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGc                                               >  4:105740/1‑147 (MQ=255)                                                                                          cTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTAGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTgg                                                                                       <  8:206028/107‑1 (MQ=255)                                                                                          cTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTAGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTgg                                                                                       <  8:201045/107‑1 (MQ=255)                                                                                          cTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTAGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTgg                                                                                       >  7:206028/1‑107 (MQ=255)                                                                                          cTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTAGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTgg                                                                                       >  7:201045/1‑107 (MQ=255)                                                                                           ttACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCaa                                             >  6:57451/1‑148 (MQ=255)                                                                                              cATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAAtgt                                          <  3:144343/148‑1 (MQ=255)                                                                                                tatCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAAtgtg                                         >  4:63202/1‑147 (MQ=255)                                                                                                 atCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAAtg                                           <  2:54604/144‑1 (MQ=255)                                                                                                 atCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAAtg                                           >  1:54604/1‑144 (MQ=255)                                                                                                 atCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGAGCAAACTGTTTGATACCCCGATGGAGCAAAAGCAGGCAATGTgg                                        >  2:284702/1‑147 (MQ=255)                                                                                                  tCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGaa                                      >  3:138687/1‑148 (MQ=255)                                                                                                  tCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGa                                       <  3:105740/147‑1 (MQ=255)                                                                                                       aGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTc                                 >  3:152275/1‑148 (MQ=255)                                                                                                         ttttGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCtt                               <  6:36620/148‑1 (MQ=255)                                                                                                          tttGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTc                              >  5:60702/1‑148 (MQ=255)                                                                                                            tGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGc                            >  5:276628/1‑148 (MQ=255)                                                                                                               tCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGt                         <  4:138687/148‑1 (MQ=255)                                                                                                                 aGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTc                        <  6:60702/147‑1 (MQ=255)                                                                                                                  gCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTCTa                      >  3:231786/1‑148 (MQ=255)                                                                                                                     gCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGc                                                   >  1:107410/1‑116 (MQ=255)                                                                                                                     gCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGc                                                   <  2:107410/116‑1 (MQ=255)                                                                                                                     gCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTCTActg                   >  3:65798/1‑148 (MQ=255)                                                                                                                       cGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGCGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTCTActgct                 >  2:151032/1‑148 (MQ=255)                                                                                                                             gCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAgtg                                                                                          <  8:100829/69‑1 (MQ=255)                                                                                                                             gCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAgtg                                                                                          >  7:100829/1‑69 (MQ=255)                                                                                                                              cTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGGTAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTCTACTGCTGCTGaa           >  2:275575/1‑147 (MQ=255)                                                                                                                              cTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTCTACTGCTGCTGAat          <  1:284702/148‑1 (MQ=255)                                                                                                                               tACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTCTACTGCTGCTGAata         >  6:29098/1‑148 (MQ=255)                                                                                                                                aCGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAAGGTGGAACGGTCTGCGCCGTCTACTGCTGCTGAatat        >  1:284242/1‑148 (MQ=255)                                                                                                                                    aCAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTCTACTGCTGCTGAATATTCCa    <  8:12947/148‑1 (MQ=255)                                                                                                                                      aGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACAGTCTTAGCCGTCTACTGCTGCTGAATATTCCATc  >  5:127328/1‑148 (MQ=255)                                                                                                                                              |                                                                                                                                            CGCTCGCTACAAGATCTGAAAGATGCGCTGAAAGGCAAAGTGCAGGAAACGCTATCTGCGGTGGCGGATGACGGTATCGAGCAGAGCGGCTTACATATCTGGAGTTTTGGTCAGCTGCCGGAAAGCTACGAACAGAAGCGTGGCAACTACAAAGTGAAGGCGTGGCCGGCGCTGGTGGATGAGCGCGACAGTGTGGCGATCAAACTGTTTGATAACCCGCTGGAGCAAAAGCAGGCAATGTGGAACGGTCTTCGCCGTCTACTGCTGCTGAATATTCCATC  >  NC_000913/1486019‑1486299

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