Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	446,976	0	A	C	40.3%	112.6 / 17.5	62	W244G (TGG→GGG)	cyoE	protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/3);  new base C (0/25);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 2.22e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAAT  >  NC_000913/446826‑447121                                                                                                                                                       |                                                                                                                                                  cagccagcagcGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAacca                                                                                                                                                   <  2:177189/151‑1 (MQ=255) cagccagcagcGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAacca                                                                                                                                                   <  1:545737/151‑1 (MQ=255)      agcagcGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAaccacc                                                                                                                                                 >  1:1191251/1‑148 (MQ=255)        cagcGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGAc                                                                                                                                             >  1:876208/1‑150 (MQ=255)         agcGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGc                                                                                                                                           >  2:290966/1‑151 (MQ=255)            gTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTa                                                                                                                                         >  2:767406/1‑150 (MQ=255)                                aaaaTCAACGGACATCATCACGGAAAGGGAAGTGATGGCGATGATAGAGAAGCCGAACAGTTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCCCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGa                                                                                                                    <  2:221950/151‑1 (MQ=255)                                aaaaTCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGaccacc                                                                                                                       >  2:609204/1‑148 (MQ=255)                                         ggACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATAc                                                                                                           >  1:978393/1‑151 (MQ=255)                                             atcatcACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATAccc                                                                                                         >  1:238596/1‑149 (MQ=255)                                             atcatcACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCa                                                                                                        <  2:542168/150‑1 (MQ=255)                                               catcaCCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCg                                                                                                     >  2:872959/1‑151 (MQ=255)                                                     cgggagGGCAGTGATGGCGATGATAGAGAAGCCGAACGGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATACCCGCGCAGAGCCATTCCTAACCCCCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCg                                                                                               <  1:672341/148‑1 (MQ=255)                                                     cGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCg                                                                                               >  1:361457/1‑151 (MQ=255)                                                     cGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCAATACCTAACCCCCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCg                                                                                               <  1:225118/151‑1 (MQ=255)                                                      gagaGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGc                                                                                              >  2:1197383/1‑151 (MQ=255)                                                        gagGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGcc                                                                                             >  2:119646/1‑150 (MQ=255)                                                          ggggAGTGAGGGAGATGAAAGAAAAGCCGAACACCTGGCCCCCCCATATTCTTTCATAACCAACTTTAAACCCGCCCAGAGCAATCCCAACCCCCCAGCCCCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAag                                                                                          <  1:405905/151‑1 (MQ=255)                                                          ggacAGTGAGGGCGATGATAGAAAAGCCAAACCGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATCCCTACCCCCCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCaa                                                                                           <  1:687270/147‑1 (MQ=255)                                                           gccAGTGAGGGCGATGATAGAGAACCCAAACAGTTGCCGCCCCCAGATTCTGTCACCACAACCTTTATACCCGCGCAGCCCCATACCAACCCCCCAGCCGCTACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAaga                                                                                         <  2:78947/149‑1 (MQ=255)                                                               gTGAGGGCGTTAAAAGAAACGCCAAACCGTTGGCGCGCCCAGATTCTGTCATCAGCAATTTTAAACCCGCGCAGACCAATACCTACCCCCCAGACGCAACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAagaga                                                                                       <  1:743103/149‑1 (MQ=255)                                                                   tGGCTATGATAGAGAAGCCAAACAGCTTGCGCCCCCAGATTTTGTCATCAACACCTTTAGAACCGCGCAGAGCCTTACCAACCCCCCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGAAACCGCCAAGAGAGAGc                                                                                   <  1:683907/149‑1 (MQ=255)                                                                    ggCGTTGATAGAGAACCCAAACAGTTGGCGCCCCCAAATTCTGTCACCAAAAACTTTAAACCCGCACAGAGCCATCCCAACCCCCCAGACGCTAACCGCCGCGGCGACCCCCAGATATTTATACCCAGCGTAACCGCCAagagag                                                                                      <  1:32377/145‑1 (MQ=255)                                                                    ggCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATc                                                                                >  1:483888/1‑151 (MQ=255)                                                                    ggCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATc                                                                                >  2:300279/1‑151 (MQ=255)                                                                      cGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATc                                                                                >  1:268821/1‑149 (MQ=255)                                                                        atgatAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGCCATCAGCAACTTTATACCCGCGCAGAGCCATCCCTACCCCCCAGACGCTACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGc                                                                             <  2:347219/150‑1 (MQ=255)                                                                           atAGAGAAGCCAAACAGCTTGCGCGCCCAGATTCTGTCACCAGCAACTTTATAACCGCGCAGAGCCATCCCTAACCCCCAGACGCTACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGc                                                                             <  2:66830/147‑1 (MQ=255)                                                                             agagAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGc                                                                        >  2:871129/1‑150 (MQ=255)                                                                               agaAGCCGAACAGCTTGCGCCCCCAGATTCTGTCATCAGCAACTTTATACCCGCGCAGAGCCATCCCAACCCCCCAGCCGCTACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGc                                                                             <  2:686653/143‑1 (MQ=255)                                                                                gaAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAAc                                                                     >  1:530962/1‑150 (MQ=255)                                                                                 aaGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACgg                                                                   >  1:621441/1‑151 (MQ=255)                                                                                  aGCCAAACGGCTGGCGCGCCAATTTTCTTTAATCACCAATTTAATACCCGCGAAAACCCATCCCAACCCCCCAGACGCTAACCGCCGCGGCGCCCACCAGATATTTATCCCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGc                                                                  <  1:634861/151‑1 (MQ=255)                                                                                         cacctggcgcgcCCAAATTCTGTCACAAGAAACTTTATACCCGCGCAGAGCCATCCCAACCCCCCAGCCGCTAACCGCCGCGGCGACCCCCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGc                                                            <  2:1182402/144‑1 (MQ=255)                                                                                         cAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCg                                                           >  1:1205796/1‑151 (MQ=255)                                                                                             ttGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCgatg                                                        >  1:330556/1‑150 (MQ=255)                                                                                              tGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCgatgat                                                      >  2:868038/1‑151 (MQ=255)                                                                                              tGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCgatgat                                                      >  1:447180/1‑151 (MQ=255)                                                                                                 gcccccAGATTTTGTCATGAGCAACTTTAGAACCGCGCAGGCCCATGCCTACCCCCCAGCCGTTAACCGCCGCGGCGACCCCCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCgatg                                                        <  1:97569/143‑1 (MQ=255)                                                                                                  cccccAGATTCTGTCATCACCATTTTTATACCCCCGCAGAGCTATCCCTACCCCCCAGCCTTTAACCCCCGCGGCGCCCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCgatga                                                       <  2:49290/144‑1 (MQ=255)                                                                                                    cccAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATCATATAc                                                  >  1:516858/1‑149 (MQ=255)                                                                                                    cccAAATCCTGTCATCAGCACCTTTAAACCCGCGCAGAGCCATACCTACCCCCCAGACGCAAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAg                                                <  1:531835/151‑1 (MQ=255)                                                                                                    ccaaaaTTCTGCAATCAGCAACTTTAAACCCGCACAGACCCATACCTACCCCCCAGCCGCAAACCGCCGCGGCGACCCCCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGatat                                                    <  1:1037502/142‑1 (MQ=255)                                                                                                     ccAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAg                                                >  1:152007/1‑150 (MQ=255)                                                                                                                  tcaGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGCTATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATCTACAGCGTGATGTGAttc                                   >  1:577212/1‑150 (MQ=255)                                                                                                                   caGCACTTTTATACCCGCGGGGACCCTTCCCTAACCCCCAGCCCCAACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGAttctt                                 <  2:1142500/151‑1 (MQ=255)                                                                                                                   caGCAACTTTATACCCGCCCAGACCCATCCCTAACCCCCAGACCCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCgtgatgtga                                      <  2:705477/146‑1 (MQ=255)                                                                                                                       accTTTAAAACCGCGCAGAGCCATACCTACCCCCCAGCCGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCg                               <  2:1191603/147‑1 (MQ=255)                                                                                                                       aattttATACCCGCGCAGACCCATACCTACCCCCCAGCCGCTACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGc                              <  1:609045/147‑1 (MQ=255)                                                                                                                       aaCTTTATACCCGCGCAGACCCATCCCAACCCCCCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGc                              <  2:1107446/150‑1 (MQ=255)                                                                                                                         cTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCAc                           >  2:765625/1‑151 (MQ=255)                                                                                                                          tttATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCAcc                          >  1:117833/1‑151 (MQ=255)                                                                                                                          ttaaaacccGCGCAGAGCCATCCCTACCCACCAGCCGCAACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCa                            <  2:502554/142‑1 (MQ=255)                                                                                                                           ttATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCg                         >  2:240784/1‑151 (MQ=255)                                                                                                                           ttAAAACCGCGCAGAGCCATCCCTACCCCCCAGCCGCTACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCAcc                          <  2:71691/150‑1 (MQ=255)                                                                                                                            tataACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCg                         >  1:523157/1‑150 (MQ=255)                                                                                                                             ataACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGa                        >  2:593442/1‑150 (MQ=255)                                                                                                                              taACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGaa                       >  2:856872/1‑150 (MQ=255)                                                                                                                               acccGCGCAGAGCCATCCCTACCCCCCAGACGCAAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGaa                       <  2:356862/147‑1 (MQ=255)                                                                                                                                  cgcgCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGcc                  >  1:819022/1‑151 (MQ=255)                                                                                                                                     gcAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCttt               >  2:1147035/1‑151 (MQ=255)                                                                                                                                         agCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTAcc           >  1:587025/1‑151 (MQ=255)                                                                                                                                            cATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACt        >  2:312448/1‑151 (MQ=255)                                                                                                                                              tcccaaaCCCCGAGGCGCTACCCGCCGCGGCGACCACCAGATATTTATACCCGGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTg       <  1:1133263/145‑1 (MQ=255)                                                                                                                                                   aaCCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAAt  >  1:827390/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                  CAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAAT  >  NC_000913/446826‑447121

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