Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,511,772	0	A	C	44.6%	92.9 / 10.8	57	V312G (GTG→GGG)	mntH	manganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base A (28/3);  new base C (0/25);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.18e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAAAGT  >  NC_000913/2511622‑2511922                                                                                                                                                        |                                                                                                                                                       agagCGATACCAAACCTTAACACTACCTGACTCAGAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGATAGC‑GAATGAAGCCCTGCCTcacca                                                                                                                                                        <  2:921749/151‑1 (MQ=255)  gagagaTACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCATAATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAAC‑GAATGAAGCCCTGCATcaccac                                                                                                                                                       >  2:272216/1‑151 (MQ=255)   agCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCaccacc                                                                                                                                                      >  1:601489/1‑151 (MQ=255)    gCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCt                                                                                                                                                     <  2:473223/151‑1 (MQ=255)     cGATACCAAAACTTAACAGTACCTAACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTg                                                                                                                                                    >  1:561143/1‑151 (MQ=255)      gATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGc                                                                                                                                                   >  2:218015/1‑151 (MQ=255)         aCCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGcccc                                                                                                                                                >  1:676753/1‑151 (MQ=255)         aCCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGcccc                                                                                                                                                >  1:45478/1‑151 (MQ=255)             aaaCTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGcc                                                                                                                                             >  1:888167/1‑150 (MQ=255)              aaCTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGcc                                                                                                                                             >  1:850959/1‑149 (MQ=255)               aCTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGcc                                                                                                                                             >  2:322083/1‑148 (MQ=255)                 ttAACAGTAACTGACTCATAACCAGAATCCGTGTCGGATCTAATCACATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGA‑GAATCAAGCCCTGCATCACCACCTGCCCCGCCAgtgt                                                                                                                                        >  2:615690/1‑151 (MQ=255)                         cccTGACTCATAACCAGAATCCGTGTCGGATCTAACCCCATCAGAATGACAATAAAGGACGGCAACATGGTGACTTTACGACGCCCCCACAGCGGGATATAAAAGC‑GAATGAAGCCCTGCATCCCCACCTGCCCCGCCAGTGTCCccacca                                                                                                                                <  1:65704/150‑1 (MQ=255)                           cTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGCGTCCCCaccacc                                                                                                                              >  2:511868/1‑151 (MQ=255)                               cTCATAACCAGAATCCGTCTCGGATCTTATCCCATCAGAATGATAATAAATGACGGCAACATGGAGACTGTACGACGCACCCACAGCGGAATATGAAAGC‑GAGTGAAGCCCTGCATACCCACCTGCCCCGCCAGTGTCCCCACCACCGtt                                                                                                                           >  2:934137/1‑150 (MQ=255)                                  aTAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAAAAAATGACGGCAACATGGTGACTGTCCGACGCCCCCACAGCGGGATATAAAAGC‑GAATGAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGaa                                                                                                                        <  1:741871/150‑1 (MQ=255)                                   tAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGa                                                                                                                         >  1:505902/1‑148 (MQ=255)                                        aGAATCCGTGTCGGATCTAATCCCATCAAAATGACAATAAATGACGGCAACATGGTGTCTGTACGACCCCCCCCCAGCGGGATATAAAGGC‑GAATGAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGc                                                                                                                  <  1:651512/150‑1 (MQ=255)                                           aTCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCg                                                                                                              >  1:229194/1‑151 (MQ=255)                                                 gtCGGATCAAACCCCACCAGAATGACAAAAAATGACGGCAACATGGTGACTGTACGCCGCCCCCACAGCGGGATATAAAAGC‑AAATGAACCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCgcagca                                                                                                         <  1:75295/150‑1 (MQ=255)                                                   cGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAAcc                                                                                                      >  1:486554/1‑151 (MQ=255)                                                   cGGATCAAATCCCCACAAAAAAAAAAAAAATGCGGACAACAGGGTCACTGTACGACGCCCCCACGGGGGGATATGAAACC‑AAATGACGCCCCGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAAcc                                                                                                      <  2:415276/151‑1 (MQ=255)                                                    ggATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAAcc                                                                                                      >  2:482485/1‑150 (MQ=255)                                                       tCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGc                                                                                                  >  1:131335/1‑151 (MQ=255)                                                          acccccATCAAAACCACAAAAAAGGACGGCAACATGGTGACGGGCGCGCGCCCCCACAGCGGGATATAAAAGC‑AAATGAACCCCAGCCCCCCCCCCCCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAgg                                                                                                   <  2:42671/144‑1 (MQ=255)                                                          aaccccATCAGAATGACAAAAAATGACGGCAACATGGTGACTGTACGACGCCCCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCACCACCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTa                                                                                               <  1:225929/148‑1 (MQ=255)                                                           accccacaaGAATGACAAAAAATGACGGCAACAGGGGGACTGTACGACGCCCCCACAGGGGGATAGGAAAGC‑GAAGGAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTaa                                                                                              <  1:125143/143‑1 (MQ=255)                                                            tCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAAt                                                                                             >  2:786149/1‑151 (MQ=255)                                                               caccAGAATGACAATAAATGGCGGCAACATGTTTCCTTTCCGCCCCCCCCCCAGCGGAATATGAAAGC‑GATTGAGCCCCTGCACCCCCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCa                                                                                          <  1:1479/148‑1 (MQ=255)                                                                   aGAATAACAATAAATGCCGGCAACGTGGTGACTGTCCGACGCCCCCACGGCGGGATATAAAAGC‑GAATGAAGCCTTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCa                                                                                          <  1:252776/147‑1 (MQ=255)                                                                   aGAAAGACAAAAAAAGCCGGCAACAGGGGGACTGGCCGCCGCCCCCACAGGGGGAGATGACAGA‑GAAGGACGCCCCGCCCCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATcc                                                                                           <  2:220999/146‑1 (MQ=255)                                                                     aaTGAAAAAAAATGCCGCCAACAGGGGAATTGTACGCCGCCCCCGCGGCGGAATATAAAGGC‑AAACAAGCCCCTCCCTCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGa                                                                                      <  1:932731/149‑1 (MQ=255)                                                                         aCAAAAAATGACGCCACCAGGGTGCCGTTACGCCCCACCCCCAGGGGGAAAT‑AAAGCGGAAGGACGCCCAACACCCCCCCCGCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTg                                                                                <  2:7244/151‑1 (MQ=255)                                                                            aaaaaTGAGGCAAAAATGGTGACGGTACGCGGCCCCCACGGTGGGATAGAAAAGC‑GAATGAAGCCATGCACCCCACCCCGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCg                                                                             <  2:236482/149‑1 (MQ=255)                                                                            aTAAATGCCGGCAACAGGGGTACTGTCCGACCCCCCCAGGGCGGGAAGTGAAGGC‑GAATGAACCCCGCCACCACCCCCGGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCg                                                                             <  2:592918/151‑1 (MQ=255)                                                                              aaaggccGGCAACATGTTGCCTGTACGACGCCCCCCCAGCGGGTTATGAACGC‑GAAGGAAGCCCTGCCCCCCCCCCCGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGc                                                                            <  2:814364/144‑1 (MQ=255)                                                                                    cGGCAACAGGGTGACTGTACGACCCCCCCCCAGCGGGATATGAAAGC‑GAAAGACCCCCACAACCACCCCCCCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCg                                                                        <  2:880044/148‑1 (MQ=255)                                                                                     ggCAACGTGGTAGCTGGAGGACGCCCCCACGGCGGTATATAAAAGC‑AAATGAAGCCCGGCACCACCCCCGCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGt                                                                       <  2:601334/148‑1 (MQ=255)                                                                                     ggCAACATGGTGCCTGTACGACGCCCCCACAGCGGGATATGAAAGA‑GAATGAACCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCg                                                                        <  1:782304/147‑1 (MQ=255)                                                                                     ggCAACATGGTGACTGTACGACGCCCCCACAGGGGGATATGAAAGC‑AAATGAACCCCTCCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTgg                                                                     <  1:5668/150‑1 (MQ=255)                                                                                     ggCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGc                                                                    >  2:756880/1‑151 (MQ=255)                                                                                                   tGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGtt                                                        >  2:50290/1‑149 (MQ=255)                                                                                                          cgCCCCCACAGCGGAATATGAAAGC‑GAAGGAACCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCg                                                  <  1:634760/148‑1 (MQ=255)                                                                                                          cccccccACAGCGGGATATAAAAGC‑GAATGAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGt                                                 <  1:163390/145‑1 (MQ=255)                                                                                                             aCCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTc                                                >  2:998037/1‑147 (MQ=255)                                                                                                              cccACAGCGGAATATAAAAGC‑AAATGAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGAt                                            <  2:945284/150‑1 (MQ=255)                                                                                                                    gCGGGATATGAAAGC‑AAAGGAAGCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtc                                     <  1:734051/151‑1 (MQ=255)                                                                                                                      gggATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtcat                                   >  1:912520/1‑151 (MQ=255)                                                                                                                                aGC‑GAAGGATTCCCTCCATCACCGCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGc                                        <  2:570858/136‑1 (MQ=255)                                                                                                                                  gcaaaggAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCgg                          <  1:650959/143‑1 (MQ=255)                                                                                                                                    gAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAac                      >  1:623378/1‑151 (MQ=255)                                                                                                                                     aagaaaGCCGTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCgg                          <  2:729775/142‑1 (MQ=255)                                                                                                                                      ataacccccTGCTCCACCCCCCGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGAt                             <  2:939394/136‑1 (MQ=255)                                                                                                                                       tGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACAcc                   >  2:146487/1‑151 (MQ=255)                                                                                                                                        gAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACAcct                  >  1:901983/1‑150 (MQ=255)                                                                                                                                         aggccccccccccccaccTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGTCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAg                 <  2:995994/138‑1 (MQ=255)                                                                                                                                          aGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGt                >  1:656080/1‑151 (MQ=255)                                                                                                                                          aGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGt                >  1:418306/1‑151 (MQ=255)                                                                                                                                               tGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGAc           >  2:235364/1‑151 (MQ=255)                                                                                                                                                ccccccccaccCGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTGGCCGCAGCGTGGCTTAACAGCGGGTGCAGCGTCAGATAAGCCTCATCAAGATCGTCAACACCAGTATGatc          <  2:243414/145‑3 (MQ=255)                                                                                                                                                ccaccaccccCCGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGGCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGa            <  2:287465/145‑1 (MQ=255)                                                                                                                                                 cATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGCCCGTTGCCGCAGCGTGGCTTAACAGCGGCTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCg         >  2:915593/1‑151 (MQ=255)                                                                                                                                                 cATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCg         >  1:182242/1‑151 (MQ=255)                                                                                                                                                      ccaccTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGaaa     >  2:453521/1‑150 (MQ=255)                                                                                                                                                        accTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAAAGt  >  1:172144/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       AGAGCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑GAATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAAAGT  >  NC_000913/2511622‑2511922

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