Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	46.2%	120.6 / 18.1	78	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (32/10);  new base C (0/36);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 1.90e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCT  >  NC_000913/786707‑787001                                                                                                                                                     |                                                                                                                                                    gTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAcc                                                                                                                                                   >  2:1145463/1‑151 (MQ=255)  tCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAccc                                                                                                                                                  >  2:374137/1‑151 (MQ=255)   cAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAccc                                                                                                                                                  <  2:432352/150‑1 (MQ=255)    agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg                                                                                                                                                >  1:243391/1‑151 (MQ=255)    agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg                                                                                                                                                <  2:633707/151‑1 (MQ=255)     gagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTgg                                                                                                                                               <  1:1248376/151‑1 (MQ=255)        tGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGttt                                                                                                                                            >  1:253854/1‑151 (MQ=255)        tGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGttt                                                                                                                                            >  1:953711/1‑151 (MQ=255)         ggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTg                                                                                                                                           >  2:5284/1‑151 (MQ=255)         ggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTg                                                                                                                                           <  2:287050/151‑1 (MQ=255)              aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑cc                                                                                                                                     <  1:571448/151‑1 (MQ=255)                ccGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑ccgc                                                                                                                                   >  1:234751/1‑151 (MQ=255)                ccGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑ccg                                                                                                                                    <  1:1145345/150‑1 (MQ=255)                   aTGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑Ccgccgc                                                                                                                                >  1:276560/1‑151 (MQ=255)                               aGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGc                                                                                                                    >  1:214772/1‑151 (MQ=255)                                    tGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGAtct                                                                                                               >  1:400094/1‑151 (MQ=255)                                                          tCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                          >  1:1030850/1‑150 (MQ=255)                                                           cgggCAATGTGCTCCATTTTTAGCAACAAAAAAGCCGACTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                          <  2:941364/147‑1 (MQ=255)                                                           caggaatggtgCTTCATTGTAAGCAAAAAAAAGGCCGTCTCCCTTGCAGTCGGTTTTTTCTTTTAAAACAAAGGACTTTTACTCCGCTTTACCCTGTTTTGCAT‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGat                                                                                             <  1:600945/138‑1 (MQ=255)                                                            aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                       >  2:401208/1‑151 (MQ=255)                                                              gCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                     >  2:1078560/1‑151 (MQ=255)                                                                   gtgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                  >  2:1252148/1‑149 (MQ=255)                                                                      cTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                             >  2:330560/1‑151 (MQ=255)                                                                      cTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                             >  1:541164/1‑151 (MQ=255)                                                                      cTCCATTGTTAGCAAAAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCTCATTTAAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCg                                                                              <  1:75921/150‑1 (MQ=255)                                                                        ccATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGt                                                                            <  2:143941/150‑1 (MQ=255)                                                                         cATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                          >  1:738723/1‑151 (MQ=255)                                                                         cATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                          >  1:443050/1‑151 (MQ=255)                                                                           ttGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTTAGCGGTTTg                                                                         >  1:295154/1‑150 (MQ=255)                                                                           ttGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                         <  2:1147919/150‑1 (MQ=255)                                                                           ttGTTAGAAACAAAAAAGCCGCCTCACTTGCAGTCGGCTTTCTCATTTAAAACAAATGACTTTTACTTCTCTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                         <  1:467777/150‑1 (MQ=255)                                                                              ttAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaaa                                                                      <  2:191643/150‑1 (MQ=255)                                                                              ttAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaa                                                                       <  1:1105848/149‑1 (MQ=255)                                                                              ttAACAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCTTTTTAAACGAATGACTTTTACTCCGCTTTCCCCTGTTTTGAAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaa                                                                       <  1:307906/149‑1 (MQ=255)                                                                                aGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtc                                                                   >  1:854324/1‑151 (MQ=255)                                                                                   aacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGat                                                                                             >  2:800638/1‑122 (MQ=255)                                                                                       aaaaaCCCGCCTCCCTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGAAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGa                                                            <  1:237491/151‑1 (MQ=255)                                                                                          aaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACt                                                         >  1:880849/1‑151 (MQ=255)                                                                                              cGACTCAGTTGCAGGCGGCTTTCTCATTTTAAACGAATGGCTTTTATTTCGCTTTACCCGGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATa                                                     <  2:703149/151‑1 (MQ=255)                                                                                                 cTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                  >  1:1014163/1‑151 (MQ=255)                                                                                                  tCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                 >  2:953947/1‑151 (MQ=255)                                                                                                   cACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTCCGCTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCa                                                <  1:291412/151‑1 (MQ=255)                                                                                                    aCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                 >  1:1072689/1‑149 (MQ=255)                                                                                                     cTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAg                                               <  2:253879/150‑1 (MQ=255)                                                                                                     cTTGCAGTCGGCTTTCTCATTTTAAACAAATAACGTTTCCTTCGCTTTCCCTGGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAg                                               <  1:1138889/150‑1 (MQ=255)                                                                                                      ttGCGGTCGGCTTTCTCTTTTTAAACGAATGACTTTTCCTTCGTTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                 <  1:987408/147‑1 (MQ=255)                                                                                                       tGCAGTCGGTTTTCTCATTTAAAACGAATGTCTTTTTCTCCGTTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                             <  1:222729/150‑1 (MQ=255)                                                                                                       tGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                            >  1:348326/1‑151 (MQ=255)                                                                                                        gCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGc                                           >  1:128878/1‑151 (MQ=255)                                                                                                              ggttttCTCATTTAAAACGAATGACTTTTACTTCTCTTTCCCCGGGTTTAAAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                            <  1:592499/141‑1 (MQ=255)                                                                                                                cTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                   >  2:296989/1‑151 (MQ=255)                                                                                                                 tttCTCTTTTAAAACAAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCa                                     <  1:666997/148‑1 (MQ=255)                                                                                                                 tttCTCATTTTAAACGAATGCGTTTTACTTCGCTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                   <  2:730380/150‑1 (MQ=255)                                                                                                                 tttCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCa                                          <  2:911762/143‑1 (MQ=255)                                                                                                                 tttCTCATTTTAAACAAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCa                                     <  2:431721/148‑1 (MQ=255)                                                                                                                 tttCCCATTTAAAACTAATAACTTTTCCGTCCTTTTCCCCTTGTGTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCa                                     <  2:738787/148‑1 (MQ=255)                                                                                                                  ttCTCATTTTAAACGAATGACGTTTACTCCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                  <  2:686151/150‑1 (MQ=255)                                                                                                                       aTTTAAAACAAAGGACGTTTCTTTGGTTTTCCCGTGTTTTCCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCa                                     <  1:40410/142‑1 (MQ=255)                                                                                                                        tttaAAACGAATGACTTTTCCTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                             <  2:578378/149‑1 (MQ=255)                                                                                                                         tttAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                          >  2:548486/1‑151 (MQ=255)                                                                                                                            aaaCGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                         <  2:953805/149‑1 (MQ=255)                                                                                                                              aCGAATGACGTTTACTTCGCTTTACCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGc                      <  2:223411/150‑1 (MQ=255)                                                                                                                                gAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                    >  1:565794/1‑150 (MQ=255)                                                                                                                                 aaTGACTTTTCCTTCGCTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                    <  1:1252022/149‑1 (MQ=255)                                                                                                                                  aTGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                   >  1:623830/1‑149 (MQ=255)                                                                                                                                    gACGTTTACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGaa               <  2:312817/151‑1 (MQ=255)                                                                                                                                     atgTTTATTCCGTTTTCCCCT‑GTTTACAACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAg                 <  1:206489/146‑1 (MQ=255)                                                                                                                                        tttCCTGCGCTTACCCCTGGTTGGCAC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCCGCGGTTTGAAATTCTCGTCGAACTCATCCACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                <  2:243416/146‑1 (MQ=255)                                                                                                                                        tttACTTCGCTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtct            <  1:5284/150‑1 (MQ=255)                                                                                                                                        tttACTTCGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGaa               <  1:806357/147‑1 (MQ=255)                                                                                                                                             tccGCTTTCCCCTGTTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATctct          <  1:693060/145‑1 (MQ=255)                                                                                                                                               cGCTTTCCCCTGTTTTGAAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc     <  2:502057/150‑1 (MQ=255)                                                                                                                                               cGCTTTCCCCTGTTTGCCAT‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc     <  2:829205/150‑1 (MQ=255)                                                                                                                                               cGCTTTCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc     <  2:400134/150‑1 (MQ=255)                                                                                                                                                gCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc   >  2:889370/1‑151 (MQ=255)                                                                                                                                                  tttCCCCTGGTTTGCAC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc   <  1:656747/149‑1 (MQ=255)                                                                                                                                                  tttCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc     <  2:779972/147‑1 (MQ=255)                                                                                                                                                  tttCCCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt  <  1:945881/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCT  >  NC_000913/786707‑787001

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