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

AGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATAC  >  NC_000913/3443129‑3443419                                                                                                                                                        |                                                                                                                                               aGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTAcc                                                                                                                                              >  2:339392/1‑151 (MQ=255)    cGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTAccccc                                                                                                                                           <  1:566563/151‑1 (MQ=255)     gCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCg                                                                                                                                          >  2:1302873/1‑151 (MQ=255)         agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCat                                                                                                                                      >  1:606392/1‑151 (MQ=255)         agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCat                                                                                                                                      >  2:633831/1‑151 (MQ=255)             tAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTc                                                                                                                                  >  1:619734/1‑151 (MQ=255)                   aTAAAGCGGTTGCCCACGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTc                                                                                                                            >  2:533296/1‑151 (MQ=255)                    tAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTc                                                                                                                            >  1:26583/1‑150 (MQ=255)                        gCGGTTTCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCgg                                                                                                                       >  1:1056387/1‑151 (MQ=255)                         cGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGt                                                                                                                      >  1:337249/1‑151 (MQ=255)                              gCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGGTTACCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                 >  1:416713/1‑151 (MQ=255)                               ccggggctaaaaaTACAGCAAATTTTTTGCAACCCT‑TTCCCAA‑CCACCCCCCCCCCCGACCCATGCCGCGATGGCCGCCGGAAAAAAAATAATCCTGGAAGCGAAGTTGGCC‑GGTATCCCCCCCGCCATATTCACTTTCAGCGGt                                                                                                                      <  1:1073300/135‑1 (MQ=255)                               cccGGGCTGAAAATACGGCAAATTTGTTGCCAGCCG‑GTCCAAA‑CCAGTCCCCCCCCCAAACCAGGACGCGATGGTCGCCGGGAACAAAATAATACTGGAAGGAAAGATGGCC‑GGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                  <  2:1121796/149‑1 (MQ=255)                               ccagggctaaaaaTACAGCAAATTTGTTGCCACCCA‑‑TTCCAACCCAGTCCCCCCCCCGAACCATGCCGCGAGGGTCCCCGGAAAAAAAATAATACGGGAAGCAAAGATGGCG‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                  <  2:257272/139‑1 (MQ=255)                               caagggctaaaaaTACAGCAAATTTGTTGTCAGCCG‑TTCCCAC‑CCAGTCCCCCCCCCGACCCATGACGCGATGGCCGCCGGAAAAAAAATAATACTGGAAGGAAAATTTGCG‑GGGTTTACCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                  <  1:969763/139‑1 (MQ=255)                                   ggcccaaaaTACGGCAAATTTTTTGCCAGCCA‑‑TTCCAACCCAGTCCCCCCCCCGACCCAGGACGCGATGGCCGCGGAAAACAAAATAATACGGAAAGGAAAGTTGGCG‑GGGATACCCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                  <  1:231943/139‑1 (MQ=255)                                       acAAATACAG‑AAAAATGGTGGCAGCCAGGTCCCAA‑CAAGTACCCCCCCCAAACAAAGACGGGAGGGTCGCCGGAAACAAAATAATATTGAAGGAGAAGTTGGCC‑GGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                         <  1:713911/149‑1 (MQ=255)                                             aCAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                                  >  2:1166769/1‑151 (MQ=255)                                             aCAACGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                                  >  2:1166710/1‑151 (MQ=255)                                                gCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                >  1:945866/1‑150 (MQ=255)                                                    aaTTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                            >  1:162260/1‑150 (MQ=255)                                                      ttgttgTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacg                                                                                         >  2:791140/1‑151 (MQ=255)                                                       tgttgtCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacg                                                                                         >  2:1037986/1‑150 (MQ=255)                                                         ttgtCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg                                                                                      >  1:1065285/1‑151 (MQ=255)                                                          tgtCAGCCA‑‑TTCCAACCCAGTCCCCCCCCCGACCCATGCCGCGAGGGTCGCGGGAAAAAAAAAAATCGGGGAAGCGAAGATGGCC‑GGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg                                                                                      <  2:394779/150‑1 (MQ=255)                                                             cAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt                                                                                  >  1:281896/1‑151 (MQ=255)                                                                 cA‑GTCCCAA‑CAAGTCCCCCCCCCGAACCAGGACGCGAGGGTCGCCGGAAAAAAAAAAATACGGGAAGCGAAGATTGCC‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                              <  2:483575/151‑1 (MQ=255)                                                                     tcccAA‑CCAGTCCCGCCCCCAAACCATGCCGCGATGGTCGCCGGAAACAAAATAATACGGGAAGCGAAGTTTGCG‑GGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                 <  2:151632/143‑1 (MQ=255)                                                                       ccaacccggtcccccccccGAACCAGGAGCTGAGGCTCGCGGAAAACAAAAAAATGCGGGAAGGGAAGTTTGCG‑GGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                           <  1:762778/136‑1 (MQ=255)                                                                       ccAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                         >  1:1327757/1‑151 (MQ=255)                                                                       ccAA‑CCAGCCCCCCCCCCAACCCAGGCCGCGTTGGCCGCGGGAAAAAAAAAAATACTGAAAGGGAAGATGGCC‑GGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACg                                                                          <  1:1097775/150‑1 (MQ=255)                                                                          a‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                      >  2:34669/1‑151 (MQ=255)                                                                          a‑CCAGCACCGCCCCCAACCAAGGGCGCGATGGTCCCGGGAAAAAAAAAAATCAGGGAAGAG‑AGATTGCCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                         <  2:932385/148‑1 (MQ=255)                                                                             cAGTCCCCCCCCAGAACCATGACGTGATGGTCGCGGGGAAAAGAATAATACTGGAAGGAAATTTTGCG‑GGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGttt                                                                       <  2:246025/148‑1 (MQ=255)                                                                             cAGCACCGCCCCCGAACCATGACGCGAGGGTCGCGGGGAAAAGAATAATACTGAAGGCGAAGATTGCC‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGttt                                                                       <  2:669171/148‑1 (MQ=255)                                                                               gccccGCCCCCGACCAATGCCGCGATGGCCGCCGGAAAAAGAATAATACTGGAAGCGAAGTTTGCC‑GGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                  <  2:1290653/148‑1 (MQ=255)                                                                                 cccccccccAACCCATGACGCGATGGTCGCCGGAAAAAAAATAATACGGGAAGCGAAGATTGCG‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAg                                                                <  2:677419/147‑1 (MQ=255)                                                                                  ccccccccAACCCATGCCGCGATGGCCGCGGGAAAAAAAATAATACGGGAAGGAAAGATGGCG‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                  <  2:1162937/145‑1 (MQ=255)                                                                                  ccGCCCCCGAACCATGCCGGGATGGTCGCGGGAAAAAAAATAATACTGAAAGCGAAGATTGCG‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTa                                                                 <  2:12926/149‑1 (MQ=255)                                                                                     cccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                             >  2:917820/1‑150 (MQ=255)                                                                                       cccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc                                                          >  2:501521/1‑151 (MQ=255)                                                                                           aaCCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                       >  2:811767/1‑150 (MQ=255)                                                                                               aTGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                  >  1:485447/1‑151 (MQ=255)                                                                                                   cgcgATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGt                                              >  1:494628/1‑151 (MQ=255)                                                                                                    gggttGGTCGCGGGGAACAAAATAATACTGGAGGCGAAGTTTGCG‑GGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                               <  1:540086/145‑1 (MQ=255)                                                                                                      gATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                            >  2:1039088/1‑150 (MQ=255)                                                                                                        tGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAcc                                         >  2:942651/1‑151 (MQ=255)                                                                                                         ggTCGCGGGGAAAAAAATAATACGGGAAGAGAAGTTTGCG‑GGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCa                                        <  1:829741/151‑1 (MQ=255)                                                                                                            cgcgGGAAACAAAAAAATCCGGGAAGGGAAGATTGCG‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                      <  2:1213273/150‑1 (MQ=255)                                                                                                               cGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                   >  1:939198/1‑150 (MQ=255)                                                                                                               cGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                   >  2:717825/1‑150 (MQ=255)                                                                                                                ggaaaaaaaataatACGGGAAGCGAAGATGGCG‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCt                                    <  2:977339/140‑1 (MQ=255)                                                                                                                      aaaataatACGGGAAGAGAAGATTGCC‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                            <  2:478063/148‑1 (MQ=255)                                                                                                                               cTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                     >  2:893561/1‑148 (MQ=255)                                                                                                                                 ggaggagAAGTTTGCG‑GGCTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc                <  1:963262/145‑1 (MQ=255)                                                                                                                                 ggAGGCGAAGATTGCG‑GGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACaaa                    <  2:528721/147‑1 (MQ=255)                                                                                                                                   aaGCTAAGGTTCCC‑GAGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt              <  2:902927/151‑1 (MQ=255)                                                                                                                                   aaGCGAAGATGGCG‑GGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg               <  2:152795/150‑1 (MQ=255)                                                                                                                                     gCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc             >  2:738837/1‑150 (MQ=255)                                                                                                                                     gCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc             >  2:251356/1‑150 (MQ=255)                                                                                                                                     gCGAAGATGGCC‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc             <  1:199808/150‑1 (MQ=255)                                                                                                                                      cGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCAc           >  2:1246856/1‑151 (MQ=255)                                                                                                                                       gAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt          >  1:1127092/1‑151 (MQ=255)                                                                                                                                         agtttGCC‑GGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt          <  2:174002/146‑1 (MQ=255)                                                                                                                                              gCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat    >  1:968555/1‑150 (MQ=255)                                                                                                                                              gCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat    >  2:79458/1‑150 (MQ=255)                                                                                                                                                  gggATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa      <  2:619742/145‑1 (MQ=255)                                                                                                                                                  gggATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatac  <  1:974886/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                               AGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCA‑GTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCC‑GGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATAC  >  NC_000913/3443129‑3443419

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