Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,574,051	0	G	T	44.6%	91.1 / 33.4	65	V384V (GTC→GTA)	pqqL	putative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base G (7/29);  new base T (29/0);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 3.43e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573908‑1574201                                                                                                                                                |                                                                                                                                                        ttCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAcatca                                                                                                                                                   <  2:386780/150‑1 (MQ=255) ttCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTacaaca                                                                                                                                                   >  2:993555/1‑150 (MQ=255)  tCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAcatcat                                                                                                                                                  <  1:550496/150‑1 (MQ=255)       aGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCCAggt                                                                                                                                            >  2:793146/1‑149 (MQ=255)        gTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCCAGtt                                                                                                                                            >  1:453432/1‑150 (MQ=255)          tGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGTTGCCAGGCGACTGGTCAGCATACGTAAATCACGTTGAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGTGTATATTTGACATCATCGAGttct                                                                                                                                          <  1:642093/150‑1 (MQ=255)                gTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGGTTTTcag                                                                                                                                      >  2:109464/1‑148 (MQ=255)                    ctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGaa                                                                                                                                 <  1:22699/149‑1 (MQ=255)                      tCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAGTTTACATCATCGAGTTTTTCAGCAGaaaaa                                                                                                                              >  1:149184/1‑150 (MQ=255)                           cGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGc                                                                                                                        <  1:922595/151‑1 (MQ=255)                            gACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTTTTCAGCAGaaaaa                                                                                                                              >  2:205780/1‑144 (MQ=255)                               aaGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGAt                                                                                                                    <  1:241029/151‑1 (MQ=255)                                aGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTATATCATCGAGTTCTTCAGCCGAAAAAACAAGCTGATc                                                                                                                   >  2:801116/1‑151 (MQ=255)                                aGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTTACATCATCGAGGTCTTCAGCAGAAAAAACATGCTGAt                                                                                                                    >  2:843362/1‑150 (MQ=255)                                 gAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCa                                                                                                                  <  1:204964/151‑1 (MQ=255)                                     ggCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTTTTTCGCAGAAAAAACAAGCTGATCAA‑Tgg                                                                                                             >  2:801812/1‑151 (MQ=255)                                     ggCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTGACATCATCGAGGTCTTCCGCAGAAAAAACATGCTGATCaa                                                                                                                 >  2:913355/1‑148 (MQ=255)                                       cGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCATCCAGTTTTTTAGCAGAAAAAACAAACTGATCAA‑TTGtt                                                                                                           >  1:971648/1‑151 (MQ=255)                                        gTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTg                                                                                                          <  2:600218/151‑1 (MQ=255)                                           ttattTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTGGATTTTAGAATATAGGGTTTTTTAGGAGAAAAAACACGATGGTTAA‑TTGTTGcc                                                                                                        >  2:754345/1‑150 (MQ=255)                                              tttAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTAAATTATCGAGGTCTTCAGCAGAAAAAACAAGCTGATCAA‑AGGGTGccca                                                                                                      >  1:709707/1‑147 (MQ=255)                                                 aaTGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCATCGAGGTCTTTAGCAGAAAAAACATGCTGATCAA‑TGGTTGCCAActc                                                                                                   >  2:299508/1‑149 (MQ=255)                                                 aaTGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTGCATCATCGAGTTCTTCAACCGAAAAAACAAGCTGATTAA‑Tgggt                                                                                                           >  2:769620/1‑139 (MQ=255)                                                  aTGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGc                                                                                                <  1:388286/151‑1 (MQ=255)                                                     aGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGAAGGGTTTCTAATATCCAGGTTTTTAGCACAAAAAACACGCTGATCAA‑TGGTTTCCaaaact                                                                                                  >  1:1033168/1‑142 (MQ=255)                                                       cTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTAAATCATCGAGTTCTTTCGCAGAAAAAACAAGATGATCAA‑TGGTTGGCAAATCTTCCCtt                                                                                            >  1:884678/1‑150 (MQ=255)                                                        tGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTCTTTAGCAGAAAAAACATGCTGATCAA‑TTGTTGCCAACTCTTCCAtt                                                                                            >  2:259975/1‑149 (MQ=255)                                                        tGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTTACATCATCGAGTTTTTTAGCAGAAAAACCATGCTGGTCAA‑TGGTTGGCCACTCTGCCCtt                                                                                            >  1:841342/1‑149 (MQ=255)                                                         ggATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTaa                                                                                          <  1:308441/150‑1 (MQ=255)                                                           aTGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGc                                                                                       <  1:993575/151‑1 (MQ=255)                                                                  gCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCgca                                                                                <  1:722736/151‑1 (MQ=255)                                                                  gCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCgca                                                                                <  1:656702/151‑1 (MQ=255)                                                                   cGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGTCATCATCCAGTTTTTTAGCAGAAAAAACAAGCTGATCAA‑TGGTTGCCCACTCTTCCCTTAATACact                                                                                    >  2:1150993/1‑144 (MQ=255)                                                                    gACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTATACTATCTAGGTTTTTTGCAGAAAAAACATGCTGATCAA‑ATGTTTCCCACTCTTCCCTTTAttc                                                                                       >  1:1026852/1‑140 (MQ=255)                                                                    gACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTAAATCATCGAGTTTTTCTTCAGAAAAAACATGCTGATCAA‑TGGTTTCCAACTCTTCCATTAATGCAttttca                                                                                >  1:1038816/1‑145 (MQ=255)                                                                        ggTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGTTTTTTAGCCGAAAAAACACGCTGATCAA‑TGGGTGCCAACTCTGCCCTTTATGCATTCGCAGCAccc                                                                          >  1:858828/1‑151 (MQ=255)                                                                             gCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCaa                                                                                                                 <  1:675594/108‑1 (MQ=255)                                                                              cATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCata                                                                    <  1:869033/151‑1 (MQ=255)                                                                               aTACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCatat                                                                   <  2:209167/151‑1 (MQ=255)                                                                               aTACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCatat                                                                   <  1:240702/151‑1 (MQ=255)                                                                               aTACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCatat                                                                   <  2:443802/151‑1 (MQ=255)                                                                                     aaaTCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGt                                                             <  1:664918/151‑1 (MQ=255)                                                                                     aaaTCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTAAATCATCGCGTTCTTTAGCAGAAAAAACATGCTGATCAA‑TGGTTGCCAACTCTTCCCTTTATTCATTCCCAGCATCCTGCATATTTTCGt                                                             >  2:610860/1‑151 (MQ=255)                                                                                        tCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATTATCGAGTTTTTTAGCAGAAAAAACAAGCTGATCAA‑TGGTTGCCAACTCTGCCCTTAATGCATTTGCAGCAACCTGCATATTATCGtct                                                           >  1:433421/1‑150 (MQ=255)                                                                                           cGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTTTTTAACAGAAAAACCACGCTTATTAA‑TGGTTGCCAACTCTTCCCTTAATGCATTCCCAGCATCCTGCCTATTATTGtct                                                           >  2:282573/1‑147 (MQ=255)                                                                                                   ttGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGCGTTTTTTCGCCGAAAAACCATGCTGATCCA‑TGGTTTCCAACTCTTCCCTTTATTCATTTGCAGCATCCCTCATATTTTTGTCTCGTGCCttt                                                  >  1:401530/1‑148 (MQ=255)                                                                                                    tGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTTTTTAGCAGAAAAAACAAGCTGATCAA‑AGGTTGCCAACTCTTCCCTTTAAGTATTCTCCGCATCCTGCATATTTTTGTCTCGTggcttt                                                  >  1:962469/1‑142 (MQ=255)                                                                                                    tGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGTTATAATCGAGTTTTTTAGCAGAAAAAACAAGCTGATCAA‑TGGTGTTCAACTCTGCCCTTAATACATTCGCAGCATCCTGCATATTATCGTCTCCTGCATTTtca                                               >  1:633060/1‑150 (MQ=255)                                                                                                    tGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTac                                                <  1:108781/149‑1 (MQ=255)                                                                                                     gCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGGCAACATCGAGGTCTTTTCCCGCAAAAACAAAATGATCAATTGGTTCCCAAC‑CTGCCCTTTATTCATTCTCAGCACCCTGCATATTTTCTCCTCCTGCAtttt                                                 >  2:1095739/1‑146 (MQ=255)                                                                                                      cTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGTGATTTTACATCATCGAGTTCTTCCGCAGAAAAAACAAGCTGATCAA‑TTGTTGCCAACTCTGCCCTTAATGCATTTGCAGCATCCTGCATATTATCGTCTCGTGCGttt                                                  >  2:250645/1‑145 (MQ=255)                                                                                                           cAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGaaaaa                                       <  2:344689/151‑1 (MQ=255)                                                                                                            aaCCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGTCCTCACCCCGGTCTTTAGCCGAAAAAACATGCTGATCAA‑AGGGTGTCAACTCTTCCCTTAATGCATTCTCAGCAACCTGCATATTATCGTCTCGCGGATTTTCACCGaaaaa                                       >  1:952816/1‑150 (MQ=255)                                                                                                              ccGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAg                                    <  1:295941/151‑1 (MQ=255)                                                                                                               cGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGc                                   <  1:743528/151‑1 (MQ=255)                                                                                                                        cAGCCAGGTGAGGCGGGTAGATTTTATACCATCCAGGTCTTCCCCAGAAAAACCAAGCTGGTCAA‑AGGGTGGCCACTCTTCCCTTTATTCCTTCTCCACATCCTCCATATTATCGTCTCGTGCACTTTCACCGaaaaaaaaccgctg                              >  1:638493/1‑138 (MQ=255)                                                                                                                             aGGTGAGGCGGGGAGATTTTAACTCATCGAGGTCTTTCGCAGGAAAAACAAGCTGGTCAA‑TTGGTGCCAACTCTGCCATTAATGCAATTTCAGCAACCCCCATATTATTGTCTCGTGCATTTTTAAGGAAAAAAAAAGACTGATAATTggg                     >  1:891829/1‑151 (MQ=255)                                                                                                                              ggTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGt                    <  1:869319/151‑1 (MQ=255)                                                                                                                               gTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTg                   <  1:205771/151‑1 (MQ=255)                                                                                                                                     cGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAttt             >  2:134414/1‑151 (MQ=255)                                                                                                                                     cGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAttt             >  1:1121106/1‑151 (MQ=255)                                                                                                                                     cGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAtt              <  1:856881/150‑1 (MQ=255)                                                                                                                                          aGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAAc         <  2:328753/150‑1 (MQ=255)                                                                                                                                           gATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGc       <  2:704553/151‑1 (MQ=255)                                                                                                                                               tGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgc    <  2:1104649/150‑1 (MQ=255)                                                                                                                                                gACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcgc  >  2:276310/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                        TTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAA‑TGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573908‑1574201

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