Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,261,055	0	A	C	41.2%	109.4 / 13.6	68	V455G (GTG→GGG)	tdcE	pyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (33/7);  new base C (0/28);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.98e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATT  >  NC_000913/3260907‑3261203                                                                                                                                                      |                                                                                                                                                     gTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCAc                                                                                                                                                     <  2:1174909/150‑1 (MQ=255)      gAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgcc                                                                                                                                               >  1:213611/1‑151 (MQ=255)           atatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTa                                                                                                                                           >  1:1244784/1‑150 (MQ=255)                         gATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagcagc                                                                                                                            >  2:981803/1‑151 (MQ=255)                           tttccgccACCCCCAGCCAGTCCATGAAGTGACCGGGGCTCTCCATCCCTTTGTCGTAACCCACCACGTCGTCCACCAGCGGTGCTTTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTagag                                                                                                                                <  2:548384/138‑1 (MQ=255)                           tGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGt                                                                                                                          >  2:769108/1‑151 (MQ=255)                            ttACTGCACCCCCAGCCATTCCATGAAGTGATCGGGGCTGTCCATCTCTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTCCTGTTTT‑CGCCCCGCCCTGATTCTTCGCCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGagcag                                                                                                                             <  1:928894/146‑1 (MQ=255)                             tACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGttt                                                                                                                        >  2:449680/1‑151 (MQ=255)                                tGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCg                                                                                                                     >  2:1153436/1‑151 (MQ=255)                                  ccccccaagccatTCAATGAAGTGACCGAGGCTTTCCATCACTTTTTCGTAACCCAGCACGTCGCCCACCGGGGGTGTTGTTTG‑GGCCCCGACCTGAATTTTCAGTTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                      <  2:1283912/141‑1 (MQ=255)                                    ccgccgGCCAGTCCATAAGGTGACCGAGGCGGTCCACCTCTTGGTCGTAACCCAGCACGTCGTCCATCAGCGGTGCTTTTTT‑CGCCCCGACCTGAATCTTCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                                                                                   <  1:1260760/149‑1 (MQ=255)                                     ccccacccagCCCAGAAGGTGACCGGGGCTGTCCATCTCTTTGTCGAACCCCCGCCCGTCGTCCATCAGGGGTTTTTTTTC‑CGCCCCGCCCTAATTCTTCTCCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                                <  2:133420/149‑1 (MQ=255)                                      gccagccaGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                                                                                               >  1:1003644/1‑151 (MQ=255)                                        cagccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                                                                              >  2:706106/1‑150 (MQ=255)                                           ccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                           >  1:1048992/1‑150 (MQ=255)                                                     aaGTGATCGAGGCTGCCCATCTTTTTCTCGTACCCCAGCACGTCGTCCATCAGGGGTGTTTTTTC‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCa                                                                                                    <  1:716622/147‑1 (MQ=255)                                                      aGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa                                                                                                >  2:162848/1‑150 (MQ=255)                                                         gATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGTCCGAACTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGa                                                                                             >  2:1108918/1‑150 (MQ=255)                                                            cGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCTTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACt                                                                                          >  2:385995/1‑150 (MQ=255)                                                              aGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                                                        >  1:597128/1‑150 (MQ=255)                                                                     ccATCTCTTTGTCGTAATCCAGCACGTCGTCCACCAGCGGTGCTTTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCt                                                                                 <  1:259095/150‑1 (MQ=255)                                                                     ccATCCCTTTGTCGTAACCCGGCCCGTCTTCCATCAGCGGTGCGGTTTT‑CGCCCCGGCCGGATTCTTCCCCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAAc                                                                                           <  2:714733/140‑1 (MQ=255)                                                                     ccATCATTTTGTCGTAACCCACCACGTCGCCCACCGGCGGTG‑TGTTTTGCGGCCCGCCCTGATTCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGc                                                                                  <  2:907981/149‑1 (MQ=255)                                                                      cATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt                                                                                >  1:856216/1‑150 (MQ=255)                                                                         cccTTTGTCGAAACCCAGCCCGTCGTCCATCGGCGGTGCTTTTTT‑CGGCCCGGCCGGAATCTGCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCa                                                                            <  1:133256/149‑1 (MQ=255)                                                                         cattttGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTTTTTT‑CGCCCCGACCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCa                                                                            <  1:1000902/148‑1 (MQ=255)                                                                            tttGTCGTAACCCGGCACGTCGTCCACCAGCGGTGCTGTTTT‑CGCCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                                                             <  1:1220469/147‑1 (MQ=255)                                                                              tGTCGTAACCCCGCCCGTCCTCCACCAGCGTTGCTGTTTT‑CGGCCCGCCCTGAATCTTCACCTCCCCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCa                                                                        <  2:268137/150‑1 (MQ=255)                                                                                  gTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAtt                                                                   >  2:855278/1‑151 (MQ=255)                                                                                  gTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAtt                                                                   >  2:668266/1‑151 (MQ=255)                                                                                  gTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAtt                                                                   >  1:381678/1‑151 (MQ=255)                                                                                       ccacctcttcgtcCATCAGCGGTGCTTTTTT‑CGGCCCGACCTAAATCTTCATCTTCGCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTg                                                                  <  2:940958/139‑1 (MQ=255)                                                                                       ccAGCCCGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACTTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                               <  1:1103016/150‑1 (MQ=255)                                                                                       ccAGCACGTCGTCCATCAGCGGTTCTTTTTG‑CGCCCCGCCCTAAATTTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCt                                                              <  2:1068756/151‑1 (MQ=255)                                                                                          gCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGAc                                                           >  1:602996/1‑151 (MQ=255)                                                                                          ccaccccccccAACAGTGGTGTTGTTTT‑CGCCCCGCCCTGATCCTCCATCTTCTCCCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCt                                                              <  1:906633/139‑1 (MQ=255)                                                                                           cACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGAc                                                           <  1:445377/150‑1 (MQ=255)                                                                                           cACGTCGCCCACCAGCGGTGCTTTTTC‑CGGCCCGCCCTCATTCTTCGGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTg                                                             <  2:359908/148‑1 (MQ=255)                                                                                             cgtcgtCCATCGGCGGTGCTGTTTT‑CGGCCCGACCGGATTCTTCAGTTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGAc                                                           <  2:840546/148‑1 (MQ=255)                                                                                                cgtcCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                                      >  1:596423/1‑150 (MQ=255)                                                                                                 gtcCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                                      >  1:1109705/1‑149 (MQ=255)                                                                                                  tcCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAgg                                                   >  1:1160209/1‑151 (MQ=255)                                                                                                   ccATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGcagcag                                                    <  1:1172047/149‑1 (MQ=255)                                                                                                        aGCGGTGCTGTTTT‑CGCCCCGACCTGACTCTTAGGCTTCCCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAt                                               <  2:251810/149‑1 (MQ=255)                                                                                                         gCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                              <  2:703043/149‑1 (MQ=255)                                                                                                                 gTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgt                                    >  1:929515/1‑151 (MQ=255)                                                                                                                    tt‑CGGCCCGACCTGAATCTTCCGCTCCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcg                                  <  2:711777/150‑1 (MQ=255)                                                                                                                         gCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGtt                             >  1:770700/1‑151 (MQ=255)                                                                                                                         ccccGCCCTGATTCTCCGGTTCCTCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGtt                             <  2:123182/150‑1 (MQ=255)                                                                                                                          cccGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTg                            >  1:1029581/1‑151 (MQ=255)                                                                                                                             gACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGa                           >  1:1267976/1‑149 (MQ=255)                                                                                                                               ccTGAATCTTCAGCTTCTCTTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                       <  1:583191/151‑1 (MQ=255)                                                                                                                                  gAATCTTCAGTTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCa                      <  1:149616/149‑1 (MQ=255)                                                                                                                                   aaTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGt                    <  2:458797/150‑1 (MQ=255)                                                                                                                                     tCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                 >  2:1265909/1‑151 (MQ=255)                                                                                                                                      cTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGc                >  1:389130/1‑151 (MQ=255)                                                                                                                                      cTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGc                >  1:273573/1‑151 (MQ=255)                                                                                                                                        tCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCAt              >  1:481078/1‑151 (MQ=255)                                                                                                                                         cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc             >  1:906130/1‑151 (MQ=255)                                                                                                                                         cAGCTTCTCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCAt              <  1:1205061/150‑1 (MQ=255)                                                                                                                                           gCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCa            <  1:449609/150‑1 (MQ=255)                                                                                                                                               ctcctcCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCa            <  2:1048921/146‑1 (MQ=255)                                                                                                                                               ctcGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc             >  2:1108824/1‑145 (MQ=255)                                                                                                                                               ctcGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAg           <  1:449697/147‑1 (MQ=255)                                                                                                                                               cccGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGa          <  2:884315/146‑1 (MQ=255)                                                                                                                                                  gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTc     >  2:620581/1‑150 (MQ=255)                                                                                                                                                    ccACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt  >  1:255247/1‑151 (MQ=255)                                                                                                                                                    ccACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAt   <  1:456313/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTT‑CGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATT  >  NC_000913/3260907‑3261203

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