Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,261,842	0	T	G	46.7%	91.6 / 12.0	60	T12P (ACC→CCC)	prs	phosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (8/24);  new base G (28/0);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 3.80e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACA  >  NC_000913/1261696‑1261990                                                                                                                                                    |                                                                                                                                                     acaAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCg                                                                                                                                                  <  2:999958/151‑1 (MQ=255)    aGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTTGGAATACGTTGTGCTAG‑TTTCCGGGGGGCggt                                                                                                                                                >  1:846076/1‑148 (MQ=255)       ggACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCa                                                                                                                                            <  1:523286/151‑1 (MQ=255)       ggACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAAACGTTGTGCTAG‑TTCCGGGGGGGGGTTAcc                                                                                                                                             >  1:268357/1‑150 (MQ=255)       ggACTGGATGATGAAAATATCACCACCGCGTACAATTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGGTGTGCTTG‑TTTCGGGGGGGGGTTTcc                                                                                                                                             >  2:59571/1‑150 (MQ=255)         aCTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAAACGTTGTGCTAA‑TTCCGGGGGGGGGGTACCAGc                                                                                                                                          >  1:935973/1‑151 (MQ=255)            ggATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCaa                                                                                                                                        <  2:51512/150‑1 (MQ=255)               tgatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATAAGGTGGGCCAG‑TTCCGGGGGGGCGGTACCACCaaaaa                                                                                                                                     >  1:869685/1‑150 (MQ=255)                gatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAg                                                                                                                                    <  2:899716/150‑1 (MQ=255)                 atgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTTG‑TTCCGGGGTGGGGGTTACAGCcaaaa                                                                                                                                     >  2:826540/1‑148 (MQ=255)                   gaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGGTACCAGCAAAAAACtt                                                                                                                                 >  2:802609/1‑150 (MQ=255)                   gaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATAAGTTGTGCTAG‑TTCCGGGGGGGCGTTAACAGCAAAAAGCtt                                                                                                                                 >  1:356524/1‑150 (MQ=255)                      aaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCata                                                                                                                             <  2:732322/151‑1 (MQ=255)                        tatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAACGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGGTACCAGCAAAAAGCTTCttat                                                                                                                            >  1:957622/1‑150 (MQ=255)                           caccacCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCTGCGTCGCCGAGTGAAGTGCACAGGCGGCAGGGCAGACGGTGCGGTAGTTTCCGCGGGGGCGTTA‑CAGCTAAAAGTGTCat                                                                                                                              >  1:42582/1‑145 (MQ=255)                           caccacCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAgg                                                                                                                        <  2:544759/151‑1 (MQ=255)                                ccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGGGGTACCCGCAAAAAACTTCATATCCGGCAc                                                                                                                     >  1:1030842/1‑149 (MQ=255)                                     tACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGTTACCAGCAAAAAACTTCATATCAGGCACgaga                                                                                                                 >  2:879199/1‑148 (MQ=255)                                      aCATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGGGTTACCAGCAAAAAGCTTCATATTAGGCCCGagaag                                                                                                               >  2:1011547/1‑149 (MQ=255)                                        aTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAAcc                                                                                                           <  1:754644/151‑1 (MQ=255)                                        aTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGCTGTGCTAG‑TTCCCGGGGGGCGGTACCAGCAAAAAACTTCATATTAAGCACGAgaagaa                                                                                                             >  1:792305/1‑149 (MQ=255)                                            tCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCa                                                                                                        <  2:702316/150‑1 (MQ=255)                                             cATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAg                                                                                                       <  2:761878/150‑1 (MQ=255)                                             cATTAATTTGTACACTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTTCGGGGTGGGGTTACCAGCAAAAAACTTCCTATCAGGCACGAGAAGAAACTCAgg                                                                                                      >  1:490330/1‑151 (MQ=255)                                              aTTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCCATACGGTGTGCTAG‑TTCCGGGGGGGCGGTACCAGCAAAAAACTTCATATCAGGCACGAGAAGAAc                                                                                                            >  2:609648/1‑144 (MQ=255)                                                 aaTTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATAAGTTGTGCTTG‑GTCCGGGGGGGAGTTGGCCACAAAAAGCTTCATATAAGGCa                                                                                                                      >  2:302108/1‑131 (MQ=255)                                                         cGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCAt                                                                                           <  1:73736/150‑1 (MQ=255)                                                           cTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTTGCAATACGGTGGGCCAG‑TTCCGGGGGGGGGGTACCCGCAAAAAGCTTCATATCAGGGACGAGAAGAAACTCAGGGATGCGTCCATtgg                                                                                        >  1:127403/1‑151 (MQ=255)                                                               cTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACCTTGTTGTAG‑TTTCTGGGGGGCGGTACCAGCAAAAAACCTCATATCAGGCACGAGAAGAAAATCAGGCCTGCGTCCATtggggg                                                                                     >  1:709909/1‑150 (MQ=255)                                                                 tCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGa                                                                                    <  1:285280/149‑1 (MQ=255)                                                                 tCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGGTACCAGCAAAAAACTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCCTtggtgg                                                                                     >  1:369050/1‑148 (MQ=255)                                                                  cGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAg                                                                                 <  2:68074/151‑1 (MQ=255)                                                                       tCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑GTCCGGGGGGGGGGTACCAGCAAAAAACTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCCTTGGGGGAAAGaa                                                                               >  1:254121/1‑148 (MQ=255)                                                                           tAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGGTAG‑TTCCGGGGGGGCGTTTCCAGCAAAAAGCTTCATATAAGGCACGAGAAGAACCCCAGGCAAGCGCCCCTTGGTGGGAAGAATTTGcc                                                                         >  1:634766/1‑150 (MQ=255)                                                                             aaGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGATGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGCGAAGAAACTCAGGCATGCGTCCATTTGGGGGAAGAATCTGCCg                                                                        >  2:985618/1‑149 (MQ=255)                                                                               gCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGaaaa                                                                    <  2:804112/151‑1 (MQ=255)                                                                                 gACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATAAGTTGTGCTAG‑TTCCGGGGTGGCGGTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAAACTCCGGGATGCGTCCATTTGTGGGAAGAAACTGCCGGAAAc                                                                   >  2:500071/1‑150 (MQ=255)                                                                                  aCCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAAAAGTTGGGCTAG‑TTCCGGGGGGGGGGTAACAACAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGGATGGGTCCATTGGGGGGAAAAATATGCCGaaaa                                                                    >  1:680998/1‑148 (MQ=255)                                                                                   ccTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACtgt                                                                <  1:878505/151‑1 (MQ=255)                                                                                        aGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGGTACCAGCAAAAAACTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGGAAGAATCTGCCGAAAACTGTGCggg                                                           >  2:748671/1‑151 (MQ=255)                                                                                            gcgTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAg                                                       <  1:825866/151‑1 (MQ=255)                                                                                            gcgTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAg                                                       <  1:825509/151‑1 (MQ=255)                                                                                             cgTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAgg                                                      <  1:503176/151‑1 (MQ=255)                                                                                              gTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGc                                                     <  1:169654/151‑1 (MQ=255)                                                                                                cGCCGAGTGAAGTGTACAGGCGGTTTGCAATACGTTGTGCTAG‑TTCCGGGGGGGGGTTACCAGCAAAAAACTTCATATCCGGCACGAGAAGAACCTCAGGCATGCGTCCCTTGGTGGAAAGAATCTGCCGAAAAATGTGCGGGCCAGGCa                                                    >  2:220287/1‑150 (MQ=255)                                                                                                cGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCAt                                                   <  2:759501/151‑1 (MQ=255)                                                                                                 gCCGAGTGAAGTGTACAGGCGGTTGGCAATAAGTTGTGCTAG‑GTTCGGGGGGGCGTTACCCGCAAAAAACTTCATATCAAGCCCGGGAAGAAACTCAAGGAAGAGTCCCTTTGGGGGAAGAAACTGCCGAAAACTTTGCGGGGCAgg                                                      >  2:660615/1‑147 (MQ=255)                                                                                                   cGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGGTACCAGCAAAAAACTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGGAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGa                                                 >  2:807384/1‑150 (MQ=255)                                                                                                    gAGTGAAGTGTACAGGCGGTTGGCAATACGGTGTGCTAG‑GTCCGGGGGGGCGTTACCAGCAAAAAACTTCATATTAGGCACGAGAAGAAACTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAAACGTGCGGGCCAGGCAt                                                   >  1:226003/1‑147 (MQ=255)                                                                                                      gTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCt                                             <  1:222134/151‑1 (MQ=255)                                                                                                           gtgtACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCa                                        <  1:406846/151‑1 (MQ=255)                                                                                                           gtgtACAGGCGGTTGGCAAAACGTTGTGCTAG‑TTTCGGGGGGGCGGTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCCGGCCTGATCCttt                                           >  1:506975/1‑148 (MQ=255)                                                                                                             gtACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAg                                       <  2:639233/150‑1 (MQ=255)                                                                                                                         ggCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTaa                           <  2:273285/150‑1 (MQ=255)                                                                                                                               aCGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAgcgc                    >  1:669623/1‑151 (MQ=255)                                                                                                                                cGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGcgcg                   >  2:814444/1‑151 (MQ=255)                                                                                                                                      gCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCa              <  2:254688/150‑1 (MQ=255)                                                                                                                                      gCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCa              >  2:967781/1‑150 (MQ=255)                                                                                                                                          g‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAAAGt          >  1:969726/1‑150 (MQ=255)                                                                                                                                             tCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTg       >  2:122316/1‑151 (MQ=255)                                                                                                                                                    tGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACa  >  2:622590/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACA  >  NC_000913/1261696‑1261990

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