Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	345,328	0	C	A	41.9%	137.5 / 58.9	86	intergenic (+337/‑76)	yahL/yahM	uncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (45/5);  new base A (0/36);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 3.50e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345181‑345477                                                                                                                                                      |                                                                                                                                                      aaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCaa                                                                                                                                                     >  2:290748/1‑150 (MQ=255)  aaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCaaa                                                                                                                                                    >  2:1176785/1‑150 (MQ=255)     tGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAAcc                                                                                                                                                >  2:305664/1‑151 (MQ=255)      gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAAcc                                                                                                                                                <  2:113869/150‑1 (MQ=255)         ctcATTCAGATGTTTTGTGAATGTGTTGGTTGACATACATATGAAAAAAATAATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGAAAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAAcc                                                                                                                                                >  1:218942/1‑147 (MQ=255)           cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCa                                                                                                                                            <  2:635038/149‑1 (MQ=255)           cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCa                                                                                                                                            <  2:303843/149‑1 (MQ=255)           cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCAtt                                                                                                                                          >  2:1189641/1‑151 (MQ=255)               cAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCt                                                                                                                                      >  2:1240612/1‑151 (MQ=255)                aGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCtc                                                                                                                                     >  1:787305/1‑151 (MQ=255)                    gTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGAAAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCtcatc                                                                                                                                  >  1:739590/1‑150 (MQ=255)                      tttGTGAATGGGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAgg                                                                                                                               >  1:642268/1‑151 (MQ=255)                         gtgAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGa                                                                                                                              >  1:650361/1‑149 (MQ=255)                           gattttgtTGGTTTCCTTTAATATAAAAAAAAAAATAATTCCATAATGTTTGGGTGGCTTCTTTTCTTTTTCTTCCCCGCCACCCCCACCCCTGAAAGAAAATGTGCTAG‑GA‑GTTAAATAAAAAAACCCGCCATTAGCTCATCAGGAag                                                                                                                            <  1:1267156/144‑1 (MQ=255)                            aaTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                         >  2:199966/1‑151 (MQ=255)                              tgtgTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                         >  2:1330856/1‑149 (MQ=255)                                ttttggttgAAATTCATAAAAAAAAAATCATAATTCCATCTTGTTTGTGGGGGTTTTTTTTTTTATCTTACGCACCCCCCCCAAGCCGAAAAAAAAATATGCTGA‑AG‑GGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAg                                                                                                                       <  1:748438/147‑1 (MQ=255)                                tgtTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                                                     >  1:25602/1‑151 (MQ=255)                                    gtttgACATTCATATAAAAAAAATAAAAATCCCATCATTTTTGTGTGGGTTTTTTTCTTTATCTTCCGCTTCGCCCTCAAGCCTGACAAAAAATTTGCTGA‑AA‑GCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                                  <  2:1215859/148‑1 (MQ=255)                                    gtttgACATTCAAAAAAAAAAAAACATAATTCCATCATGTTTGGGTGGGGTTTTTTCTGTATTTTACCCTCCCCACTCAACCCGGAAAAAAAATATTCTTT‑AA‑GGCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACa                                                                                                                 <  2:989074/149‑1 (MQ=255)                                    ggttgACATTAATATAAAAAAAATAAAAATCCCATCATTTTTGGGGGGGTTTTTTTTTTTATCTTCCCCACCGCCCCCACCCCTAACAAAAAATTTGCTGA‑AG‑GCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACa                                                                                                                 <  1:1022208/151‑1 (MQ=255)                                            ttCATAAAAAAAAAATCATAATCCCATCATGTTTGTGGGGGTTTTTTTCTTTATCTTACGCACCGCACTAAACCCTAAAAAAAAATATGCTGA‑AG‑GCCTCATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCa                                                                                                             <  1:461735/147‑1 (MQ=255)                                             tCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                                                        >  2:849273/1‑151 (MQ=255)                                              cATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttg                                                                                                       >  2:978013/1‑151 (MQ=255)                                               atatAAAAAAAATCATAATCCCATCTTGTTTGTGTGGGGTCTTTTCTGTATCTTACGCACCGCACTCAAGCCTGAAAAAAAATATGCTGA‑AA‑GGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                                                        <  2:31905/149‑1 (MQ=255)                                                   aaaaaaaaTAATAATCCCATTATTTTTTTGTGGGTTTTTTTCTTTATTTTACGCATCGCACTAAACCCTGACAGAAAATATGCTGA‑AA‑GGCTAATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTa                                                                                                  <  2:248881/150‑1 (MQ=255)                                                     aaatcattaTAATTCCATCATGTTTGTGTGGGGGCATTTGTGTATCTTATGCATCGTTTTAAAGACTGACTGAAAATTTGCTGT‑AA‑GGCTCATATCAACAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAGTGTTGTATg                                                                                                >  1:418799/9‑151 (MQ=255)                                                      aaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                               >  1:1307233/1‑151 (MQ=255)                                                        aaaTCATAATCCCATCATTTTTGTGGGGGTTTTTTTCTGTATCTAACGCACCGCCCTAAAGCCTAAAAAAAAATATGCTGA‑AA‑GGCTCAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTa                                                                                             <  1:726186/151‑1 (MQ=255)                                                         aataatATTTCCATCATGTTTGTGTGGGTTTTTTTCTTTCTCTTACGCATCCCCCACAACCTTAAAAAAAAATATTTTGA‑AA‑GGCACATATAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAc                                                                                            <  1:461181/151‑1 (MQ=255)                                                         aaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAc                                                                                            >  2:445414/1‑151 (MQ=255)                                                         aaTCAAAATTCCATCATGTTTGTGGGGGTTTTTTTCTGTATCTTACGCCTCGCACCCAACCCTGACAAAAAAAATGCTGA‑AA‑GGTTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAc                                                                                            <  1:637947/151‑1 (MQ=255)                                                           tCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgg                                                                                          >  1:813695/1‑151 (MQ=255)                                                           tCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgg                                                                                          >  1:756018/1‑151 (MQ=255)                                                             aTAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                                                        >  1:1104891/1‑151 (MQ=255)                                                              tAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGt                                                                                       >  2:466362/1‑151 (MQ=255)                                                               aaTTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                                                      >  1:209189/1‑151 (MQ=255)                                                                 ttCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                                    >  2:361490/1‑151 (MQ=255)                                                                  tCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                                    >  1:74042/1‑150 (MQ=255)                                                                   ccATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                                    >  1:175182/1‑149 (MQ=255)                                                                   ccATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                                    >  1:1121182/1‑149 (MQ=255)                                                                     atcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGc                                                                                >  2:1228361/1‑151 (MQ=255)                                                                        atGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                             >  1:257282/1‑151 (MQ=255)                                                                         tttttGTGTGGGGTTTTTTCTGTATTTAACGCACCGCCCCCAACCCAAACAAAAAATATGCTGA‑AA‑GGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAg                                                                                  <  1:1136277/143‑1 (MQ=255)                                                                         ttttggtgtgGGTTTTTTTTTGTATCTTCCACATCGCACACAAGCCTGAAAAAAAATTTGTTGT‑AA‑GGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg                                                                                 <  1:458576/141‑1 (MQ=255)                                                                            ttGTGTGGGTTTTTTTCTGTTTCTTACGCATCGCACTCAAGCCTGAAAAAAAATATGCTGA‑AA‑GGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGAt                                                                          <  1:796195/150‑1 (MQ=255)                                                                             tgtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                                        >  2:1257864/1‑151 (MQ=255)                                                                              gtgggGGGTCTTTTCTGTATCTTACGCACCGCACTCAAGCCTGAAAAAAAATATGCTGT‑AA‑GGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                                         <  1:5524/149‑1 (MQ=255)                                                                                gtgGGTTCTTTTCTTTATCTTCCGCACCGCACCCAAGCCTGAAAAAAAATATGCTGT‑AA‑GGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGc                                                                       <  2:929329/149‑1 (MQ=255)                                                                                gtgGGTTCTTTTCTGTATCTTCCGCATCGCAATCACGCCTGAAAAAAATTATGCTGA‑AG‑GGTTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                                         <  2:110431/147‑1 (MQ=255)                                                                                 tgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt                                                                    >  2:773721/1‑151 (MQ=255)                                                                                 tgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt                                                                    >  1:92181/1‑151 (MQ=255)                                                                                  gggtttttttctttttctAACGCACCGCACCAAAGCCGAAAAAAAAATATGCTGA‑AG‑GGCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCgg                                                                     <  1:421496/143‑1 (MQ=255)                                                                                  gggtttttttCTGTTTTTAACCCACCGCCCTCAAGCCTAAAAAAAAATATGCTGA‑AA‑GGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGc                                                                       <  1:1172614/141‑1 (MQ=255)                                                                                  ggggTTTTTTCTGTATCTTACGCACCGCACTCAAGCCTGAAAAAAAATATGCTGT‑AA‑GGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc                                                                   <  2:113258/151‑1 (MQ=255)                                                                                        ttttCTGTATCTTACGCATCGCACTAAAGCCGGAAAAAAAATATGCTGA‑AA‑GCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACt                                                              <  2:82345/150‑1 (MQ=255)                                                                                          ttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtttt                                                           >  2:1368086/1‑151 (MQ=255)                                                                                           tCTGTATTTTACCCACCGCCCCAAACCCCAAAAAAAAATATGTTAT‑AA‑GTCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtttt                                                           <  1:162965/150‑1 (MQ=255)                                                                                           tCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTc                                                          >  1:888331/1‑151 (MQ=255)                                                                                            cTGTATCTTACGCCCCGCAATCAAGCCGAAAAAAAAATATG‑TGTAAA‑GTCCCATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtttt                                                           <  2:521552/149‑1 (MQ=255)                                                                                            cTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTc                                                          >  1:1106386/1‑150 (MQ=255)                                                                                              gTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaa                                                        >  1:914027/1‑150 (MQ=255)                                                                                                   ttACGCACCGCACTCAAGCCTGACAAAAAATATGCTGT‑AA‑GGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaa                                                       <  2:369475/146‑1 (MQ=255)                                                                                                      cGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtcttct                                               >  2:636753/1‑151 (MQ=255)                                                                                                       gCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTg                                              >  2:244836/1‑151 (MQ=255)                                                                                                        cATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTg                                              >  2:1201526/1‑150 (MQ=255)                                                                                                             cACTAAAGCGTAAAAAAAAATATGTTGT‑AG‑GGCAAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCa                                         <  1:811234/150‑1 (MQ=255)                                                                                                              aCTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          >  1:981521/1‑148 (MQ=255)                                                                                                                  aaGCCTAACAAAAAATGT‑CTGG‑AAGGGCTAATACCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATccc                                   <  1:632082/151‑1 (MQ=255)                                                                                                                  aaGCCTAAAAAAAAATATGCTGT‑AA‑GGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATc                                     <  2:1216707/149‑1 (MQ=255)                                                                                                                   aGCCTGACAGAAAATGTGTTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATccc                                   <  2:148323/150‑1 (MQ=255)                                                                                                                   aGCCTGAAAAAAAATATGTTGA‑AA‑GGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATccc                                   <  2:428310/150‑1 (MQ=255)                                                                                                                     ccTGAAAAAAAATATGCTGA‑AA‑GGCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACAGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                                <  2:323893/151‑1 (MQ=255)                                                                                                                         aaaaaaaaTATGTTGT‑AA‑GGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtta                               <  1:1176776/144‑1 (MQ=255)                                                                                                                           aaaaaaTAGGCTGA‑AG‑GGCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatc                          <  1:290747/149‑1 (MQ=255)                                                                                                                           aGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatc                          >  1:1207127/1‑151 (MQ=255)                                                                                                                              aaaTATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTATTCTGAATCAGATCCCCTTATTatcatc                       >  2:30441/1‑151 (MQ=255)                                                                                                                                   tGCTGT‑AG‑GCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTg                   <  1:383864/150‑1 (MQ=255)                                                                                                                                    gCTGA‑AG‑GGTTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa                 <  2:787312/151‑1 (MQ=255)                                                                                                                                     cTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa                 >  1:886221/1‑150 (MQ=255)                                                                                                                                      tGT‑AA‑GCCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt                    <  1:1190463/146‑1 (MQ=255)                                                                                                                                        t‑AA‑GGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa                <  1:199966/148‑1 (MQ=255)                                                                                                                                          aa‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg            >  2:1179808/1‑151 (MQ=255)                                                                                                                                          aa‑GGCTAATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAg              <  1:1263247/149‑1 (MQ=255)                                                                                                                                          aa‑GCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGc                  <  2:837774/145‑1 (MQ=255)                                                                                                                                           a‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg           >  2:42372/1‑151 (MQ=255)                                                                                                                                                 caaaccAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGaa         <  2:1214481/143‑1 (MQ=255)                                                                                                                                                      cAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGcc  >  2:1226673/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      AAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGT‑AA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345181‑345477

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