Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	41.7%	86.0 / 20.1	60	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (7/28);  new base G (25/0);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 7.75e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508511‑3508806                                                                                                                                                   |                                                                                                                                                      tCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGgaaa                                                                                                                                                    >  1:140675/1‑150 (MQ=255)    aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGgaaaa                                                                                                                                                   >  2:633286/1‑148 (MQ=255)    aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGgaaaa                                                                                                                                                   >  1:492731/1‑148 (MQ=255)         gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcacaca                                                                                                                                           <  2:576360/151‑1 (MQ=255)          tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGgaaaaa                                                                                                                                                  >  2:877699/1‑142 (MQ=255)            tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACg                                                                                                                                         <  2:774229/150‑1 (MQ=255)            tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCaaacac                                                                                                                                          >  1:428251/1‑149 (MQ=255)               aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                                                                     <  1:343870/151‑1 (MQ=255)                aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATcaaaaa                                                                                                                                           >  1:503282/1‑140 (MQ=255)                   ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACCGTAGgccc                                                                                                                                 >  1:993842/1‑151 (MQ=255)                   ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAACACACGGGAGcccc                                                                                                                                 >  2:1082615/1‑151 (MQ=255)                    cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGgcccc                                                                                                                                >  1:246522/1‑151 (MQ=255)                    cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGAGCGCCCGGGAAAAAGCAAACACAGTaaggccc                                                                                                                                >  2:136008/1‑145 (MQ=255)                       agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCcaca                                                                                                                             <  2:105481/151‑1 (MQ=255)                           aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATcacacaacggaaccccc                                                                                                                                >  1:55394/1‑133 (MQ=255)                           aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCACATc                                                                                                                           >  1:962893/1‑149 (MQ=255)                             cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATACACACACGGTAGCCCCCACATCAGa                                                                                                                        <  2:1159882/150‑1 (MQ=255)                               atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACACGGTAGCCCCCCCATCAGaaa                                                                                                                      >  1:1215170/1‑150 (MQ=255)                                     aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACgggg                                                                                                                >  1:366488/1‑150 (MQ=255)                                            cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGcac                                                                                                         <  2:706636/150‑1 (MQ=255)                                                         gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGGAGCCCCCCCAACAGAAAGCGTGGTGGGCACACCCAGGTGCGt                                                                                             >  1:336767/1‑149 (MQ=255)                                                         gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                            <  1:727925/150‑1 (MQ=255)                                                         gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                                                           <  2:5397/151‑1 (MQ=255)                                                           aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                                                         <  1:828003/151‑1 (MQ=255)                                                           aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                                                         <  2:242565/151‑1 (MQ=255)                                                              ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                      <  1:155810/151‑1 (MQ=255)                                                                gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcat                                                                                    <  1:444067/151‑1 (MQ=255)                                                                 tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACAAACGGTAGCCCCCACATCAGAAAGCGGGGTGCGCACACCGAGCTGCGTAGCCatca                                                                                     >  1:272025/1‑149 (MQ=255)                                                                 tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGCGGAAAAAATCACACAACGGAGCCCCCCCATCAGAAAACGTGGGGCGGACACCCAGCTGCGGAGGGatcat                                                                                    >  1:75204/1‑150 (MQ=255)                                                                   ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                                                                 <  2:82387/151‑1 (MQ=255)                                                                   ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGTGCGCGCGGGAAAAATTACACACGGGAGGCCCCCCATCAGAAAGGGGGGGGGGCACCCCGCGCGGCGGAGCGGTCATccc                                                                                 >  1:457908/1‑151 (MQ=255)                                                                               aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                                                                     <  1:581294/151‑1 (MQ=255)                                                                               aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                                                                     <  1:692938/151‑1 (MQ=255)                                                                                      gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                               <  1:918731/150‑1 (MQ=255)                                                                                         tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCCGGGAAATTCACACACGGTAGCCCCCACCACACAAAGCGTGGGGCGCACACCGAGCTGCGGATCGGTCATTCCCTCCTGATAAACCAGCCCTGGc                                                           >  1:140465/1‑151 (MQ=255)                                                                                          ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                          <  1:83136/151‑1 (MQ=255)                                                                                          ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                          <  1:424520/151‑1 (MQ=255)                                                                                             gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAGCCCCCCCATCAGAAAGCGTGGGGGGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAACACCGGcacc                                                        >  1:1218072/1‑150 (MQ=255)                                                                                               aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGGGGGAAAAATTAAACACGGGAGCCCCCCCATCAAAAAAGGGGGGGGGGACACCCAGCGGCGTGGCGATCACCCCCTTCTGGTAAAACAgacatggcccc                                                        >  2:1078968/1‑138 (MQ=255)                                                                                               aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACAAGGGAGCCCCCCCCTCAGAAAGCGTGGGGCGCCCACCGCGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACCGGcacc                                                        >  2:254883/1‑148 (MQ=255)                                                                                                   cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCGCGGAAAAAATCACACACCGGAGCCCCCCCATCAGAAAGCGGGGGGGGCACACCGAGCTGCGCAGCGATCAACCCTTGCTGATAAAACAGCACTGGGCCCACtt                                                    >  1:141687/1‑148 (MQ=255)                                                                                                   cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTCACACACGGTAGCCCCCCCCACCGAAAGCGTGGTGCGCCCACCGCGCCGCGGAGCGATCAACCCTTTCTGATaaaaaaccacc                                                              >  2:1026620/1‑130 (MQ=255)                                                                                                    ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                <  2:428239/151‑1 (MQ=255)                                                                                                      agggaggAAAGAAAAGAGCGACCATTTTGTGGCAGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                              >  1:446760/1‑151 (MQ=255)                                                                                                       gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAACACCGGCCCCACTTCGCggg                                              >  1:1217930/1‑150 (MQ=255)                                                                                                       gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                             <  1:39594/151‑1 (MQ=255)                                                                                                       gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                             <  2:250167/151‑1 (MQ=255)                                                                                                             aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                       <  2:366477/151‑1 (MQ=255)                                                                                                               agaaaagaGCGACCATTTTGTGGCGGTGGGGGCCGGGAAAAATCACAAACCGGAGCCCCCACCTCAAAAAAGGGGGTGGGGACACCCAGCCGGGCAGCGGACATTCCTTTCTGGTAAACCAGCACCGGGACCACtttcccgggatgggggg                                     >  1:66538/1‑136 (MQ=255)                                                                                                               agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGc                                     <  2:503263/151‑1 (MQ=255)                                                                                                               agaaaagaGCGACCATTTTGTGGCGGATGGGCGCGGGAAAAATTAAAAACGGGAGCCCCCCCACCAGAAAAGGGGGGGGCCGCACCCAGCTGCGAAACGCTCAACCCTTTCCGATAAAACAACACTGGCACCACTCCGCgggggttggtg                                      >  2:532654/1‑142 (MQ=255)                                                                                                                gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCt                                    <  1:836345/151‑1 (MQ=255)                                                                                                                 aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACCGGAGCCCCCCCATACGAAAACCTGGGGCGCACACCCAGCCGCGGAGCGGTCATCCCCTTCTGATAAAACAACCCTGGCCCCCCCTCGCGGGGAtgg                                         >  1:133642/1‑145 (MQ=255)                                                                                                                   aagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGTAGCCCCCCCCTCAGAAAGCGGGGGGGGCAAACCCAGCTGCGTAGAGGTCATCCCTTTCTGCTAAACCAACACCGGCCCCACCTCCCGGGGATGGGGGCtgt                                  >  2:89631/1‑150 (MQ=255)                                                                                                                   aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                 <  1:416829/151‑1 (MQ=255)                                                                                                                    agagCGACCATTTTGTGGCGGGGGGGCGCGGAAAAATTTACACACCGTAACCCCCCCATCAGAAAACGGGGTGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCCCCACCTCCCCGGGATGGGGGGtg                                   >  2:273039/1‑148 (MQ=255)                                                                                                                     gagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAc                               <  2:367222/151‑1 (MQ=255)                                                                                                                      agCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAGCGGGGTGGGCACACCCCGCTGCGTAGCGGTCATCCCTTTCTGGTAAACCAGCACCGGCCCCACTTCCCGGGGATGGTGGCggtggc                               >  2:163797/1‑148 (MQ=255)                                                                                                                       gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                             <  1:136001/151‑1 (MQ=255)                                                                                                                         gACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAACGATCATCCCTTTCTGATAAACCAGCACCGGCACCACTTCGCGGGGATGGTGGGTGTGGcc                              >  1:983574/1‑148 (MQ=255)                                                                                                                                   tggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                 >  1:432190/1‑151 (MQ=255)                                                                                                                                   tggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                 >  1:313200/1‑151 (MQ=255)                                                                                                                                               gcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc     >  2:1225612/1‑151 (MQ=255)                                                                                                                                                cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg    <  1:877704/151‑1 (MQ=255)                                                                                                                                                   aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGCATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  >  2:831724/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                      TCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508511‑3508806

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