Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,443,277	0	A	C	53.2%	68.5 / 100.6	62	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (26/3);  new base C (0/33);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.74e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-13

GCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATAC  >  NC_000913/3443127‑3443425                                                                                                                                                        |                                                                                                                                                     gCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTa                                                                                                                                                      <  1:763199/151‑1 (MQ=255)    gACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAccc                                                                                                                                                   >  2:914724/1‑151 (MQ=255)     aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc                                                                                                                                                  <  2:1020229/151‑1 (MQ=255)      cGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAccccc                                                                                                                                                 >  1:626365/1‑151 (MQ=255)           agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCACCat                                                                                                                                            <  2:1058512/151‑1 (MQ=255)           agagTAACAAATAAAGCGGTTCCCCAGGCTACAAATACAGCAAATTTGTTATCACCCAGTTCCAC‑CCAGTACCGCCCCCGAACCATGCCGCGATGGTCGCCGGAAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCat                                                                                                                                            <  2:86577/151‑1 (MQ=255)            gagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCat                                                                                                                                            >  2:905798/1‑150 (MQ=255)                aaCACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTc                                                                                                                                        >  1:959392/1‑150 (MQ=255)                     aTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTc                                                                                                                                  >  2:929722/1‑151 (MQ=255)                                 ccggggctaaaaaaACAGAGAAATGTTTGCCAGCAAGTCCAAA‑CCAGTCCCCCCCCCAACCCAGGACGCGAGGGTCCCCGGAAAAAAAATAATCCGGGAGGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                                      <  1:678425/141‑1 (MQ=255)                                 cccGGGCTGAAAATACAGCAAATTTTTTGTCGGCCAGTCCCAA‑CCAGTCCCCCCCCCGAACCATGACGCGTTGGTCGCGGGAAACAAAAAAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTaa                                                                                                                          <  1:921585/147‑1 (MQ=255)                                       cTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                                                                                                >  2:1075730/1‑151 (MQ=255)                                           aaaTCCAGCAAATTTTTTGCCACCCAGTCCCAC‑CCAGCACCCCCCCCAACCCATGACGCGATGGCCGCCGGGAACATAATAATCCTGAAAGAGAAATTGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                             <  1:723208/150‑1 (MQ=255)                                           aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                            >  1:631134/1‑151 (MQ=255)                                           aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                            >  1:1006732/1‑151 (MQ=255)                                              tACAGCAAAATTGTTGCCAGCCAGTCCAAA‑CCAGTACCCCCCCCAACCCAGGACGCGATGGTCGCGGGAAAAAAAATAATACGGAAAGGGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgca                                                                                                          <  1:425524/150‑1 (MQ=255)                                                 aGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                      >  1:210684/1‑151 (MQ=255)                                                 aGCAAAATTTTTGTCAGCCATTCCCAA‑CCAGTCCCCCCCCCGACCCATGCCGCGATGGTCGCCGGAAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagca                                                                                                       <  1:846107/150‑1 (MQ=255)                                                   cGAATTTGTTGTCGGCCAGTCCCAC‑CCGGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAAAATAATACTGGAAGCGAAGATTGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcag                                                                                                         <  1:523986/146‑1 (MQ=255)                                                    gAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                      >  2:684610/1‑148 (MQ=255)                                                            tgccacccAGTCCCAA‑CCAGTCCCCCCCCCGACCCAGGACGCGATGGTCGCCGGAAAAAAAATAATACTGGAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacga                                                                                              <  2:441455/142‑1 (MQ=255)                                                                    cagtCCAACCCAGTACCGCCCCCAAACCATGCCGCGTTGGTCGCGGGAAAAAAAATAATACTGGAAGCGATGATTGCCGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                       <  2:659498/145‑1 (MQ=255)                                                                      ttCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                                  >  1:734394/1‑150 (MQ=255)                                                                      ttCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                                 >  2:1077697/1‑151 (MQ=255)                                                                         cAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                              >  2:187021/1‑151 (MQ=255)                                                                            ccagcccccccccccGAACCATGACGGGGTGGGGGCAGGGAACAAAATAATTGAGGAAGGAAAGAGTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg                                                                                   <  2:325341/134‑1 (MQ=255)                                                                              cAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgc                                                                          >  2:728317/1‑151 (MQ=255)                                                                               aGTCCCCCCCCCAAACCATGACGCGATGGTCGCGGGAAAAAAAAAAATACTGGAAGCAAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgc                                                                          <  1:337762/150‑1 (MQ=255)                                                                               aGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                                         >  2:73569/1‑151 (MQ=255)                                                                               aGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                                         >  1:175790/1‑151 (MQ=255)                                                                                  aCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                        >  2:484468/1‑149 (MQ=255)                                                                                   ccccccccGACCCATGCCGCGAGGGTCGCGGGGAAAAAAAAAATCCTGGAAGCGAAGATTGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAg                                                                      <  1:370441/147‑1 (MQ=255)                                                                                   ccccccccGACCCAGGACGGGATGGCCGCGGGGAAAAAAATAATACGGGAAGCGAAGTTTGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTa                                                                       <  2:409721/146‑1 (MQ=255)                                                                                    cccccccGACCCATGCCGCGATGGTCGCCGGAAAAAAAAAAATACTGGAAGGGAAGTTTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                                    <  1:281616/149‑1 (MQ=255)                                                                                    cccccccGACCAAGGACGCGATGGCCGCGGGAAAAAAAAAAATACGGGAAGCGAAGATGGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                                    <  1:409835/149‑1 (MQ=255)                                                                                            acccATGACGCGTGGGTCGCCGGAAAAAAAATAATACTGGAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                              <  1:597851/147‑1 (MQ=255)                                                                                            aaCCAGCACGCGTTGGTCGCCGGAAACAAAATAATCTTGGAAGCGAAGATTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                             <  2:260626/150‑1 (MQ=255)                                                                                              ccATGACGCGATGGCCGCCGGAAACAAAATAATACTGAAAGCGAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                           <  1:913200/150‑1 (MQ=255)                                                                                                aTGACGCGATGGTCGCCGGAAAAAAAATAATACTGGAAGCGAAGTTTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                        <  2:765687/151‑1 (MQ=255)                                                                                                    cgcgATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGt                                                    >  2:39465/1‑151 (MQ=255)                                                                                                     gcgAGGGCCGCGGGGAAAAAAATAATACGGGAAGGGAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcgg                                                       <  1:542085/147‑1 (MQ=255)                                                                                                          ggTCGCGGGAAACAAAATAATACTGGAAGCGAAGTTTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                                <  2:868463/149‑1 (MQ=255)                                                                                                          ggTCGCGGGAAAAAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGa                                                 <  1:177089/148‑1 (MQ=255)                                                                                                            tCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                            >  2:292352/1‑151 (MQ=255)                                                                                                               cgggaaaaaaaataatCCGGGAAGGAAAGATGGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                                     <  1:1000957/129‑1 (MQ=255)                                                                                                                        aaaaaatccgggAAGAGAATATTGCCGGTATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                <  1:593805/141‑1 (MQ=255)                                                                                                                         aataagACTGGAAGGAAAGTTTGCCGGTATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                <  2:597492/150‑1 (MQ=255)                                                                                                                          ataataCGGGAAGAGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaa                              <  2:11316/151‑1 (MQ=255)                                                                                                                          aaaatcgtGGAGGCGAAGTTGGCCGGAGTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaa                              <  1:1063180/144‑1 (MQ=255)                                                                                                                               aCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                           >  2:1061708/1‑149 (MQ=255)                                                                                                                                cTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                           >  2:889681/1‑148 (MQ=255)                                                                                                                                    aaGCAAAGATTGCCGGCATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg                     <  1:504956/150‑1 (MQ=255)                                                                                                                                     aGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc                   >  2:1081042/1‑151 (MQ=255)                                                                                                                                         aattttGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCa                  <  2:600146/144‑1 (MQ=255)                                                                                                                                         aaGTTTGCCGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCa                  <  2:720080/148‑1 (MQ=255)                                                                                                                                         aaGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg               >  1:820792/1‑151 (MQ=255)                                                                                                                                            aTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa             <  1:929751/150‑1 (MQ=255)                                                                                                                                             ttGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa             <  1:895044/149‑1 (MQ=255)                                                                                                                                                  gggATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaata         <  1:282736/148‑1 (MQ=255)                                                                                                                                                  gggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact       >  2:424440/1‑150 (MQ=255)                                                                                                                                                      ttACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatac  >  2:83656/1‑151 (MQ=255)                                                                                                                                                      ttACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatac  >  1:562073/1‑151 (MQ=255)                                                                                                                                                      ttACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaa     <  1:674223/148‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                     GCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATAC  >  NC_000913/3443127‑3443425

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: ‑