Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,443,277	0	A	C	50.7%	74.3 / 120.1	67	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (31/2);  new base C (0/34);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 4.43e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-10

AGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTA  >  NC_000913/3443129‑3443421                                                                                                                                                        |                                                                                                                                                 aGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑Acc                                                                                                                                                >  1:1042879/1‑151 (MQ=255)   aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑Accc                                                                                                                                               >  2:328615/1‑150 (MQ=255)      cATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGc                                                                                                                                           >  1:823811/1‑151 (MQ=255)      cATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGc                                                                                                                                           <  2:691903/151‑1 (MQ=255)       aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGcc                                                                                                                                          <  1:865669/151‑1 (MQ=255)        tAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGcc                                                                                                                                          >  1:651523/1‑150 (MQ=255)             tAACACATACAGCGGTTGCCCAGTCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATTACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGCTTGCCGG‑GAGT‑ACCCCCGCCata                                                                                                                                       >  1:673056/1‑148 (MQ=255)                cacaTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACt                                                                                                                                 >  1:1056026/1‑151 (MQ=255)                   aTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTc                                                                                                                              >  1:207250/1‑151 (MQ=255)                    aaaagcgtggccccAGGCTGAAAAAACGGAGAATTTGTTGTCACCCAGTCCCAACCAGCCCCCCCCCCAAACCAGGCCGCGTTGGCCG‑CGGGAAAAAAAAAAAGCCGGAAGGCGAAAATGGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCa                                                                                                                             <  1:213264/140‑1 (MQ=255)                     aaaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCTTTTTCACTTTc                                                                                                                              >  1:766925/1‑149 (MQ=255)                          ggTTGCCCGGTCAGAAAAAACAGCAAATTTTTTTCCAGCCTTTCCCACCCAGCCCCCCCCCCAAACAAGGCCGCGTGGGTCGCCGGAAAAAAAATAATACGGGAA‑GCGAAGATGGCCGG‑GATC‑CCCCCCGCCATATTCACTTTCAGCg                                                                                                                          <  1:795389/148‑1 (MQ=255)                               cccAGGCTGAAAATACAGCAAAATTGTTGCCAGCCAGTCCCAACCAGTACCGCCCCCGAACCATGACGCGAGGGCCGCGGGAAAAAAAAAAATACGGGAA‑GCGAAGATGGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                                  <  2:252011/151‑1 (MQ=255)                               cccAGGCTAAAAATACAGCAAATTTGTTGCCAGCCAGTCCCAACCAGTCCCGCCCCCGAACCAGGACGGGATGTTCGCCGGGAACAAAATAATCCTGGAA‑GCGAAGATGGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                    <  2:192510/149‑1 (MQ=255)                               ccagtgccaaaaaTACAGAGATTTTGTGGCCTGCCTATCCCAACCAGCCCCCCCCCCAAACAATGACGCGATGGCCGCCGAAAAAAAAATAATGGAGAAC‑GGGAAG‑TTGCCGGTTATT‑CCCCCCGCCATATTCACTTTCAGCGGTaaa                                                                                                                     <  1:95461/138‑1 (MQ=255)                                  aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                               >  2:734078/1‑151 (MQ=255)                                   ggcggaaaaTACAGCAAAATTGTTGTCAGCCATTTCCAACCAGTCCCCCCCCCAAACCATGCCGCGATGGTCGCGGGAAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTaaa                                                                                                                     <  2:608800/138‑1 (MQ=255)                                         aaaTACAGAAAATTTGTTGCCGGCCGGTCCCACCCGGTCCCCCCCCCGAACCATGACGGGTGGGTCGCGGGGAAAAAAATAATGCGGGAA‑GCGAAGATGGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                           <  2:502211/148‑1 (MQ=255)                                                 cGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                  >  2:139103/1‑149 (MQ=255)                                                      ttgttgTCAGCCATTCCCACCCTGTCCCCCCCCCGACCCGTGCCGGGATGGTCGCGGGAAAAAAAAAAATCGTGGAA‑GCGAAGATTGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAg                                                                                                <  1:623506/146‑1 (MQ=255)                                                        gttgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                         >  2:8962/1‑151 (MQ=255)                                                         ttgcCACCCAGTTCCAACCAGTCCCCCCCCCAAACCATGACGCGATGCTCGCGGGAAACAAAATAATGCTGGAA‑GCGAAGATTCCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                         <  2:539349/150‑1 (MQ=255)                                                                ccAGTCCCACCCAGCCCCGCCCCCCAACCAGGACGCGAGGGCCGCCGGAAACAAAAAAATCCTGGAA‑GAGAAGATTGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                   <  1:918261/149‑1 (MQ=255)                                                                        aaCCAGTCCCCCCCCCGAACCATGACGCGAGGGTCGCCGGAAACAAAATAATACTGGAG‑GCGAAGATTGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                             <  2:555825/147‑1 (MQ=255)                                                                        aaCCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGttt                                                                         >  2:684015/1‑151 (MQ=255)                                                                            agccccGCCCCAAACCCAGGACGGGAGGGTCGCCGGAAAAAGAATAATACGGGAA‑GCGAAGATTGCCGG‑GTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                        <  2:197671/144‑1 (MQ=255)                                                                                ccccccccGAACCATGACGCGATGGTCGCCGGAAACAAAATAATACGGGAA‑GCGAAGATGGCCGG‑GA‑T‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                                <  2:148450/148‑1 (MQ=255)                                                                                ccGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGt                                                                 >  2:408673/1‑151 (MQ=255)                                                                                     cccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc                                                            >  2:927506/1‑151 (MQ=255)                                                                                         aaCCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                         >  2:255833/1‑150 (MQ=255)                                                                                              tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcgg                                                   >  1:435708/1‑151 (MQ=255)                                                                                                 cgggATGGTCGCGGAAAAAAAAAAAATACGGGAG‑GCGAAG‑TTGGCGG‑GTTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                                 <  1:459322/147‑1 (MQ=255)                                                                                                 cgcgATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGt                                                >  1:122160/1‑151 (MQ=255)                                                                                                  gggagggcccccGGGAAAATAAAAATTCTGGAA‑GCAAATTTGGCCGG‑TATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACGCGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                               <  1:210902/141‑1 (MQ=255)                                                                                                  gcgATGGCGGCGGGAAACAAAATAATCCTGAAA‑GCGAAGTTTGCGGG‑GTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                                 <  1:961677/149‑1 (MQ=255)                                                                                                       ggTCGCCGGGAAAAGAATAATACTGGAA‑GCGAAGATTGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                              <  2:773378/147‑1 (MQ=255)                                                                                                        ggcgcgggaaaaaaaataatACTGGGG‑GCGAGGATGGGCGG‑ATTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc                                         <  1:771382/137‑1 (MQ=255)                                                                                                         tCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                        >  1:505486/1‑151 (MQ=255)                                                                                                             cGGGAAAAAAATAATACTGGAA‑GCGAAGTTTGCCGG‑GTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                     <  2:843308/150‑1 (MQ=255)                                                                                                               ggAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTca                                    >  1:15349/1‑149 (MQ=255)                                                                                                                 aaaaaaataatAGTGGAG‑GAAAAGATTGCGGG‑GTTC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTca                                    <  1:824952/142‑1 (MQ=255)                                                                                                                 aaaaaaataatAGGGGAA‑GCGAAGATTGCCGG‑GATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                                  <  1:748102/144‑1 (MQ=255)                                                                                                                  aCAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                                >  2:352308/1‑150 (MQ=255)                                                                                                                    aGAATAATACTGGAA‑GCAAAGATTGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                             <  2:186351/151‑1 (MQ=255)                                                                                                                     aaaaaatcCTGGAA‑GCGAAGTTTGCCGG‑GATA‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACaaa                               <  1:517569/144‑1 (MQ=255)                                                                                                                         aataCTGGAG‑GCGAAGATTGCCGG‑AATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                         <  2:119806/150‑1 (MQ=255)                                                                                                                             cggaaa‑GCGAAGATGGCCGG‑TTTC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAg                     <  2:381501/146‑1 (MQ=255)                                                                                                                             cTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                       >  2:779695/1‑148 (MQ=255)                                                                                                                               ggAG‑GCGAAGATGGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                    <  1:253989/149‑1 (MQ=255)                                                                                                                               ggAA‑GCAAAGATTGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGCTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc                  <  2:721460/151‑1 (MQ=255)                                                                                                                                gAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg                 >  1:281472/1‑151 (MQ=255)                                                                                                                                 aa‑GAGAAGTTTGCGGG‑AATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc                  <  2:59697/149‑1 (MQ=255)                                                                                                                                       aagttggCCGG‑GATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt                <  1:994555/140‑1 (MQ=255)                                                                                                                                        aGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc          >  2:1046033/1‑151 (MQ=255)                                                                                                                                        aGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc          >  1:778597/1‑151 (MQ=255)                                                                                                                                          aTTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa         >  2:368208/1‑150 (MQ=255)                                                                                                                                           ttGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa         <  1:212978/149‑1 (MQ=255)                                                                                                                                           ttGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa        >  2:875772/1‑150 (MQ=255)                                                                                                                                              cggg‑gATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat      <  1:8962/147‑1 (MQ=255)                                                                                                                                              ccGG‑TATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat      <  2:72601/149‑1 (MQ=255)                                                                                                                                              ccGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatac    >  2:143642/1‑151 (MQ=255)                                                                                                                                              ccGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatac    >  2:440767/1‑151 (MQ=255)                                                                                                                                               cGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact   >  1:294246/1‑151 (MQ=255)                                                                                                                                               cGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact   >  2:967761/1‑151 (MQ=255)                                                                                                                                                gg‑gTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact   <  2:346656/150‑1 (MQ=255)                                                                                                                                                gg‑gTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact   <  2:933068/150‑1 (MQ=255)                                                                                                                                                 g‑gATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatacta  >  1:24451/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                 AGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTA  >  NC_000913/3443129‑3443421

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