Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,098	0	A	C	45.1%	106.5 / 10.7	73	Y66S (TAC→TCC)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (29/10);  new base C (0/32);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 1.15e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCG  >  NC_000913/1490951‑1491243                                                                                                                                                    |                                                                                                                                                  ttttACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATa                                                                                                                                                   <  1:697857/148‑1 (MQ=255) ttttACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATAcc                                                                                                                                                 <  1:517329/150‑1 (MQ=255)           tATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCccgcc                                                                                                                                       <  2:483979/150‑1 (MQ=255)               gCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGAtt                                                                                                                                   >  1:468565/1‑150 (MQ=255)                cAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCACCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATtat                                                                                                                                 >  2:26434/1‑151 (MQ=255)                  gCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATAc                                                                                                                               <  1:604533/151‑1 (MQ=255)                         gcgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGcc                                                                                                                        >  2:782663/1‑151 (MQ=255)                         gcgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGcc                                                                                                                        >  1:156951/1‑151 (MQ=255)                          cgcAATGGAGTTTCGTGTTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCa                                                                                                                       <  1:1121190/151‑1 (MQ=255)                             aaTGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCCCGTGCTGATGGTCGTTCGTCTTTTGTAAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCaaaa                                                                                                                    <  2:274127/151‑1 (MQ=255)                               tGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAAcc                                                                                                                  >  1:148491/1‑151 (MQ=255)                                 gAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCg                                                                                                                 >  2:1058326/1‑150 (MQ=255)                                    tttCTTGGTTTCTCCCCCCGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGGGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTTTTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCgat                                                                                                               <  1:660489/149‑1 (MQ=255)                                    tttCGTGGTTTCTCCCCCCGTGGCGTTCGGCCACCCATCACCTTAATTCATGTTTCCTGGGGAACCTAAATCCTCGTCCTGATGGTCGTGCGCCTTTTTTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCgat                                                                                                               <  2:148422/149‑1 (MQ=255)                                            tttCTTCCCACGTAGCGAGCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTAAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACggg                                                                                                      <  2:250728/150‑1 (MQ=255)                                                  cccACGTGGCGATCGGCCACTCATCAACATGTTTCTTTTTTCCTGTGGCATCTCAATTCCCGTGCTGATGTTCGTTCGTCTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTgg                                                                                                 <  2:868309/149‑1 (MQ=255)                                                   ccACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGgcgc                                                                                              >  1:1026829/1‑151 (MQ=255)                                                      cGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCAtt                                                                                           >  1:149282/1‑151 (MQ=255)                                                        tggCGTTCGCCCACTCATCAACTTTATTCATGTTCCCTGTGGCATCTCAATTCTCGTGCTAAGGGTGATTCGTCTTTTTTTAAGGTTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCAttt                                                                                          <  2:766554/148‑1 (MQ=255)                                                                ggCCACTCATCAACATGATTCATGTTTCCTGGGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTTTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACAttt                                                                                 <  2:993094/151‑1 (MQ=255)                                                                   cACTCATCAACATGATTCATTTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGt                                                                              <  2:177573/151‑1 (MQ=255)                                                                    aCTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTg                                                                             >  1:1123861/1‑151 (MQ=255)                                                                     cTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTg                                                                             <  1:307577/150‑1 (MQ=255)                                                                      tcatcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGAt                                                                           >  2:190935/1‑151 (MQ=255)                                                                       catcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGAtt                                                                          >  2:229262/1‑151 (MQ=255)                                                                          caACATGTTTCATGTTTCCTGTGGCATCTCAATTCCCGTGCTGATGGTCGTTCGTCTTTTGTAAAGGTTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTa                                                                        <  2:997622/150‑1 (MQ=255)                                                                           aaCATGATTCATGTTTCCTGTGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGa                                                                            <  2:365874/145‑1 (MQ=255)                                                                            aCATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAtct                                                                     >  1:785727/1‑151 (MQ=255)                                                                               tGATTCATTTTTCCTGTGGCATCTCAATTCTCTGGCGGATGGTCGTTCGTCTTTTGTTAAGGCGGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTcc                                                                   <  2:212040/150‑1 (MQ=255)                                                                                gtttCATGTTTCCTGGGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTcc                                                                   <  2:940683/147‑1 (MQ=255)                                                                                gATTCATGTTTCTTGGGGAATATCAATTCTCGGGCTGATGGTCTTTCGTCTTTTGTTAGGGCTAAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCt                                                                  <  2:556525/150‑1 (MQ=255)                                                                                gATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAACTACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGAGCATTTGGGACATTTGGTGATTTATCTCCt                                                                  >  1:389354/1‑150 (MQ=255)                                                                                  ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCAGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGt                                                                <  1:402229/150‑1 (MQ=255)                                                                                   ttatttttCCTGGGGCATCTCAATTCTCGTGCGGATGGTCGTTCGTCTTTTGTAAAGGCAGAAATGCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCt                                                                  <  2:908274/142‑1 (MQ=255)                                                                                   tCATTTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTACTGt                                                                <  1:85460/149‑1 (MQ=255)                                                                                   tCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGtt                                                               <  1:1058378/150‑1 (MQ=255)                                                                                   tCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGttt                                                              >  1:813027/1‑151 (MQ=255)                                                                                      tttttCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGttt                                                              <  2:237150/146‑1 (MQ=255)                                                                                      tGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTAtt                                                           >  1:387614/1‑151 (MQ=255)                                                                                      tGTTTCCTGTGGCATCTAAATTCTCGTGCGGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTAt                                                            <  1:661010/150‑1 (MQ=255)                                                                                        tttCTTGTGGCATCTCATTTCTCGTGCTGATGGTCGTTGGTTTTTTGTTAGGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTAtt                                                           <  1:1164271/149‑1 (MQ=255)                                                                                        tttCCTGTGGCATCTCAATTCTCGTGCTGATGGTCTTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTg                                                          <  1:215783/150‑1 (MQ=255)                                                                                        ttcccTGTGCCTTCACAATTCCCGGCCGGATGGTCTTTCTTTTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTAtt                                                           <  2:997072/146‑1 (MQ=255)                                                                                            cTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTg                                                     >  1:940567/1‑151 (MQ=255)                                                                                             tgtgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                                    >  2:399498/1‑151 (MQ=255)                                                                                             tgtgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                                    >  1:112459/1‑151 (MQ=255)                                                                                               tgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGcc                                                   >  2:491724/1‑150 (MQ=255)                                                                                                 gCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGt                                                <  2:348886/151‑1 (MQ=255)                                                                                                    tctcAATTCTCGTTCTGATGTTCGTTCGTCTTTTGTGAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTg                                               <  2:191364/149‑1 (MQ=255)                                                                                                        aaTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttg                                           <  1:791098/149‑1 (MQ=255)                                                                                                        aaTTCTCGTGCTGATGGTCGTTCGTCTTTTGTAAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                                             <  2:556732/147‑1 (MQ=255)                                                                                                            ctcGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtttgg                                     >  1:780618/1‑151 (MQ=255)                                                                                                            ctcGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtttgg                                     >  2:658967/1‑151 (MQ=255)                                                                                                            ctcGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtttgg                                     >  2:581700/1‑151 (MQ=255)                                                                                                              cgggCTGATGGTCGTTCGTTTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtttg                                      <  2:226509/145‑1 (MQ=255)                                                                                                              cGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtttggt                                    <  1:90022/150‑1 (MQ=255)                                                                                                                      tGGTCGTTCTTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGatgatg                             <  1:573511/149‑1 (MQ=255)                                                                                                                       ggTCGTTCGTCTTTTGTAAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGta                           <  2:965011/150‑1 (MQ=255)                                                                                                                        gTCATTCTTTTTTTATTAAGGCTGAAATTCCTCACCCCGCCTATTATTCCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgat                                 <  2:335438/143‑1 (MQ=255)                                                                                                                          cgttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAAc                       >  1:114041/1‑151 (MQ=255)                                                                                                                           gttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAAcc                      >  2:189340/1‑151 (MQ=255)                                                                                                                            ttcgtCTTTTTTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggttt                                       <  2:177140/133‑1 (MQ=255)                                                                                                                            ttcgtCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAAc                       <  2:789187/149‑1 (MQ=255)                                                                                                                             tcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCgg                    >  1:471213/1‑151 (MQ=255)                                                                                                                              cgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCggg                   >  2:1025405/1‑151 (MQ=255)                                                                                                                                  ttttGTTAAGGCTAAAACCCCCACCTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCgggg                  <  1:113622/148‑1 (MQ=255)                                                                                                                                       ttAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAAccc            <  2:158039/149‑1 (MQ=255)                                                                                                                                         aaGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTg         <  1:491622/150‑1 (MQ=255)                                                                                                                                         aaGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTgg        >  2:248260/1‑151 (MQ=255)                                                                                                                                         aaGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTgg        >  1:868873/1‑151 (MQ=255)                                                                                                                                         aaGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCg           <  2:636617/148‑1 (MQ=255)                                                                                                                                           ggCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt      <  2:468551/151‑1 (MQ=255)                                                                                                                                             cTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGttt     <  1:992440/150‑1 (MQ=255)                                                                                                                                               gAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCg  >  1:88345/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCG  >  NC_000913/1490951‑1491243

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