Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,131	0	T	C	21.2%	186.1 / 13.3	87	intergenic (+114/‑84)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (44/23);  new base C (0/18);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 1.68e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGT  >  NC_000913/3781034‑3781228                                                                                                  |                                                                                                   aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTcc                                                                                                  <  2:1090843/100‑1 (MQ=255) aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTc                                                                                                   <  1:1542038/99‑1 (MQ=255)   ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                                                <  2:107462/100‑1 (MQ=255)    cTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                               >  2:1201881/1‑100 (MQ=255)     tAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                                                                              >  1:54144/1‑100 (MQ=255)        cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                                           >  2:881392/1‑100 (MQ=255)        cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                                           >  1:1500258/1‑100 (MQ=255)        cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                                           >  1:1778128/1‑100 (MQ=255)         ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                          >  1:176631/1‑100 (MQ=255)          ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                         <  1:2244179/100‑1 (MQ=255)          ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                         <  2:2232682/100‑1 (MQ=255)            ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                                                                        <  1:716147/99‑1 (MQ=255)            ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                                                       <  1:1318202/100‑1 (MQ=255)             tcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                                                      >  2:728510/1‑100 (MQ=255)                cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑ggg                                                                                  <  1:221433/100‑1 (MQ=255)                 tCAGGGAGAGGGGACAGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑gggg                                                                                 <  2:2559553/100‑1 (MQ=255)                  cAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGt                                                                                >  2:1342067/1‑100 (MQ=255)                    gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGa                                                                              >  1:2029282/1‑100 (MQ=255)                     ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAg                                                                             >  2:792840/1‑100 (MQ=255)                      gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAgg                                                                            <  2:54144/100‑1 (MQ=255)                      gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAgg                                                                            <  1:1996939/100‑1 (MQ=255)                      gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAgg                                                                            >  2:87797/1‑100 (MQ=255)                      gagaGGGGACCGTTCGGCGCTGTATGCACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAgg                                                                            <  1:694730/100‑1 (MQ=255)                       agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAggg                                                                           <  1:277156/100‑1 (MQ=255)                        gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAgggg                                                                          >  2:150613/1‑100 (MQ=255)                        gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAgggg                                                                          <  2:1551568/100‑1 (MQ=255)                         agGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGa                                                                         >  1:229139/1‑100 (MQ=255)                         agGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGa                                                                         <  1:2446392/100‑1 (MQ=255)                            ggACCGTTTGGCGCTGTATGTACTCCCTCACTCGGAACCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGGGTGCC‑TGGGGGAGGGGaaaa                                                                      >  2:649006/1‑100 (MQ=255)                            ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGaaaa                                                                      >  1:210146/1‑100 (MQ=255)                            ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGaaaa                                                                      >  2:1848726/1‑100 (MQ=255)                             gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAg                                                                     <  2:1955627/100‑1 (MQ=255)                               ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGc                                                                   <  2:2259523/100‑1 (MQ=255)                                 gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCg                                                                 >  1:1901418/1‑100 (MQ=255)                                       cgcTGTATGTCCTCCCTCCCCCTCAAACGACACCGCACTCTTTTTTTCGCCCTCGCCCGTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCAc                                                              <  1:1091021/97‑1 (MQ=37)                                        gcTGTATGTCTCCCCCCCCTCTGAACCGCCACCGCAGTCTTTTTTTCGCCGCCCCCCCCCGGGGGAGGGGGCG‑GGGGGGGGGGGCAAAGGCCGCACTGCt                                                          <  1:1201877/100‑1 (MQ=255)                                          tgtatgtcCTCCCCCCCTCTGAACCCACACCGCCCTCTTTTTTTCCCCCTCGCCCCCCGGGGGAGAGGGCA‑GGGGTAGGGGGAAAAGGCCGCACTGCTCa                                                        <  1:1470095/100‑1 (MQ=255)                                            tatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                       >  2:1682742/1‑99 (MQ=255)                                              tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                    >  1:17808/1‑100 (MQ=255)                                              tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                    >  2:1065985/1‑100 (MQ=255)                                                taCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtt                                                  >  1:605839/1‑100 (MQ=255)                                                 aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                 >  1:2121318/1‑100 (MQ=255)                                                  cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCAcagatcgg                                                >  1:289334/1‑93 (MQ=255)                                                  cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                >  2:2081619/1‑100 (MQ=255)                                                   tccctcccTCGGAACCGCCCCCGCCCTCTTTTTTTCCCCCTCGCCCCCCGGGGGGGGGGCCC‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                               <  1:881389/100‑1 (MQ=255)                                                   tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGAGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                               >  2:694012/1‑100 (MQ=255)                                                    cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTTTCTTTTCCCTCCGGGGAGAGGGCC‑GGGGTGATGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                              >  1:140160/1‑100 (MQ=255)                                                    cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                              >  2:1189263/1‑100 (MQ=255)                                                      cTCCCTCTGAACCGACCCCGCCCTTTTTTTTTCTCCCCCGCCCCCCGGGGGGGGGGGCG‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                            <  1:87797/100‑1 (MQ=255)                                                      cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                            >  2:163673/1‑100 (MQ=255)                                                        ccctctGAAACCCCCCCCCCCTTTTTTTTTCCCCCTCCCCCCCCCGGGGGGAGGCCC‑GGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTa                                          <  1:762404/98‑1 (MQ=255)                                                        ccctcgGAACCGCCACCGCTCTCTTTTTTCCCCCCCCCCCCCCCGGGGGAGGGGCCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTa                                          <  2:1778138/98‑1 (MQ=255)                                                          ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                        >  2:196329/1‑100 (MQ=255)                                                          ctctAAAACGACCCCGCTCTTTTTTTTTCTCCCCCGCCCCGCGGGGGAGGGGGCG‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                        <  2:210146/100‑1 (MQ=255)                                                            ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                      >  2:819463/1‑100 (MQ=255)                                                             taaaaGGACGCCACACTCTTTTTTTCCCCCTTCCCCCCCCGGGGAGAGTCCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                            <  1:511774/91‑1 (MQ=255)                                                             tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                     >  1:1016111/1‑100 (MQ=255)                                                               aaaCGCCCCCGCACTCTTTTTTTCCCCCGCCCCCCCCGGGGGGGGGGGCC‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                       <  1:1682733/96‑1 (MQ=255)                                                               aaaCGACACCGCACTTTTTTTTTCCCCCCCGCCCCCCGGGGGGGGGGGCG‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                     <  1:163673/98‑1 (MQ=255)                                                                accGACCCCGCACTTTTTTTTTCTCCCCCCCCCCCCGGGGGAAAGCCCC‑GGGTGAAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                  <  1:819462/98‑1 (MQ=255)                                                                aaCGCCCCCGCCCTCTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGGCC‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                       <  2:2062650/95‑1 (MQ=255)                                                                  cGACACCGCACTCTTTTTTTCCCCCCCGCCCCCCGGGGGAGGGGGCC‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                   <  1:1928908/97‑1 (MQ=255)                                                                    acacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGAGAGGGGAAAAGGCCGCCCTGCTCACCAACTTTCTGCTATCCTACCCCAGc                              >  2:477007/1‑100 (MQ=38)                                                                      ccctcatttttttttcCCCCCTCCCCCCCCGGGGGAGGGGGCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCa                             <  2:216315/90‑1 (MQ=255)                                                                        cGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGCCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                          >  2:1563711/1‑100 (MQ=255)                                                                         gccctctTTTTTTCTCCCTCGCCCCCCCGGGGAGAGGGCG‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                         <  1:1708814/97‑1 (MQ=255)                                                                         cccctttttttttCCCCCCCCCCCCCGCGGGGGGGGGGCCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                          <  2:17808/94‑1 (MQ=255)                                                                             tctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                     >  2:2190423/1‑100 (MQ=255)                                                                               tttttttCCCCCTCGCCCCTCCGGGGGGGGGGCC‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                       <  1:93242/96‑1 (MQ=39)                                                                               ttttttcCTCCCCCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                   <  1:2289645/100‑1 (MQ=255)                                                                                ttttttCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTTCTATCCTGCCCCCCCACAAAGGGGGCa                  >  2:672936/1‑100 (MQ=255)                                                                                 tttttCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                 >  2:2324206/1‑100 (MQ=255)                                                                                 tttttCTCCCTCGCCCCCCGGGGGAGAGGCCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                 <  1:881542/100‑1 (MQ=255)                                                                                 ttttcCCCCCTCCCCCCCCCGGGGAGAGGCCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                 <  1:1987743/100‑1 (MQ=255)                                                                                   tttCTCCCCCGCCCCCCCGGGGGGAGGGCC‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa               <  1:165619/100‑1 (MQ=255)                                                                                        ccccccccccccGGGGGAGGGGGCG‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa          <  2:2029293/94‑1 (MQ=255)                                                                                        cccTCTCCCCGGAGGGGCGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa          >  2:809176/1‑100 (MQ=255)                                                                                        atctgcccccTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa          <  2:2146886/94‑1 (MQ=255)                                                                                         ccTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa         >  2:943028/1‑100 (MQ=255)                                                                                          cTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa        >  2:2203392/1‑100 (MQ=255)                                                                                            gcccccTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaagat      >  1:2146874/3‑97 (MQ=255)                                                                                            cGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa      >  1:394921/1‑100 (MQ=255)                                                                                             gCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt     >  2:1723890/1‑100 (MQ=255)                                                                                              ccccccGGGGGAGAGGGCG‑GGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt     <  1:792839/99‑1 (MQ=255)                                                                                              ccccTCGGGGGGGAGGCCT‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt     <  2:2329134/99‑1 (MQ=255)                                                                                              ccccTCGGGGGAGAGGGCC‑GGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt     <  1:5395/99‑1 (MQ=255)                                                                                              ccccTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc    >  2:2121514/1‑100 (MQ=255)                                                                                              ccccTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc    >  2:2119906/1‑100 (MQ=255)                                                                                               ccttCCGGGGAGGGGCCC‑GGGGTGAGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg   <  2:568593/97‑1 (MQ=255)                                                                                                ccTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt  >  2:245594/1‑100 (MQ=255)                                                                                                  |                                                                                                   ACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGT  >  NC_000913/3781034‑3781228

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