Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,941	0	A	C	41.2%	146.6 / 12.3	86	E392D (GAA→GAC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (41/9);  new base C (0/35);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 7.96e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAAT  >  NC_000913/3026793‑3027089                                                                                                                                                     |                                                                                                                                                      ccTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAAcc                                                                                                                                                     >  1:95259/1‑151 (MQ=255)  cTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAAccc                                                                                                                                                    >  2:777425/1‑151 (MQ=255)      aaaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGATCCC‑Acc                                                                                                                                                <  2:1026889/150‑1 (MQ=255)             tGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTc                                                                                                                                         >  1:258168/1‑150 (MQ=255)              gCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCg                                                                                                                                       >  2:592764/1‑151 (MQ=255)                aaTTACAGGCCTACCGCCTCTCGGCATATGAAGCTTTTTACCGGGCCCGGCCCGGCGGAGCAAAATCTCTGGGCCTGGCGGATTTGTTTGGCAACTTTTTCCCGGGCAAAGAGGCTGATTCGGTGGTGAGGAACCCC‑ACCGCCACTCCGc                                                                                                                                      <  1:144207/150‑1 (MQ=255)                aaTTACAAGGCTATCGCCTCGCGGAATATGAGGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGCCGATTTGATGGGCAACTTTTTCCCTGGCAAAGAGCCTTATTTCGTGTTGAGGGACCCC‑CCCGCCACTCCGCt                                                                                                                                     <  1:954102/151‑1 (MQ=255)                aaTTACAAGGCTATCCCCCCTCGGCATATAAGGCTTTTTCCCTGCCCCCGCCCGGCGGAGCAAAATCTCTGGGCCTTGACGATTTGATTGGCACCTTTTTCCCGGCAAAAGAGGCTGATTTCGTGGTGATGGACCCC‑CCCGCCACTcc                                                                                                                                        <  1:372367/148‑1 (MQ=255)                  ttACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTa                                                                                                                                    <  2:57129/150‑1 (MQ=255)                      aaGGCTATCGCCTCTCGGCATATGAGGCGTTTTCCCTGGCCACGCTCGGCGGAGCGAAATCTCTGGCCCTTGACGATTTGATTGGCACCTTTTTCCCTGCCAAAGAGGCTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTAcagc                                                                                                                                <  1:412324/150‑1 (MQ=255)                       aGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTAcagc                                                                                                                                >  1:253108/1‑149 (MQ=255)                        ggCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACagcagc                                                                                                                             >  2:1242679/1‑151 (MQ=255)                          cTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTg                                                                                                                           >  2:1244602/1‑151 (MQ=255)                           tATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTgc                                                                                                                          >  1:251113/1‑151 (MQ=255)                            aTCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTgcg                                                                                                                         >  1:1106470/1‑151 (MQ=255)                                    tcGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGAc                                                                                                                  >  1:495733/1‑150 (MQ=255)                                      ggCATATGAGGCGTTTTACCGGCCCACGCGCGGCGGACCAAAATCTCTGGGCTTTGACGTTTTGTTTGGCAATTTTTTCCCGGGCAAGGAGCTTGTTTTCGTGGTAATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                               <  1:1120912/151‑1 (MQ=255)                                      ggCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                               >  1:81529/1‑151 (MQ=255)                                        cATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActc                                                                                                             >  1:831971/1‑151 (MQ=255)                                          tatGAAGCTTTTTCCCGGGCCACGCTCGGCGGAGCGAAATCTCGGGGCCTTGACGATTTGATGGGAAACTTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActc                                                                                                             <  2:1293218/149‑1 (MQ=255)                                          tatGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActct                                                                                                            <  2:313835/150‑1 (MQ=255)                                           atGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGCCTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGt                                                                                                          >  1:88647/1‑151 (MQ=255)                                           atGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGt                                                                                                          >  1:1241976/1‑151 (MQ=255)                                              aaGCTTTTTCCCGGGCCCCGCGGGGCGGAGCGAAATCTCGGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGTTGTTTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTc                                                                                                       <  2:1196623/151‑1 (MQ=255)                                               aGCGTTTTACCTGGCCCCGCTCGGGGGGGGTAAATCTCTGGCCCTTGACTTTTTGATTGGAATTTTTTTACCTGCCAAAGAGTTTGTTTTGGTGTTGATGTACCCC‑TCCGCCACTCCGCTACAGCAGCTGCGCTATGACAActc                                                                                                             <  2:853400/144‑1 (MQ=255)                                                 cGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                    >  2:839353/1‑151 (MQ=255)                                                   ttttCCCAGGCCCCGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGTTTGGCAACTTTTTCCCGGACAAAGGGCCTGATTTCGTGGTGATGGACCCC‑CCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTa                                                                                                   <  2:965319/150‑1 (MQ=255)                                                   ttttACCTGGCCACGCTCGGCGGGGCGAAATCTCGGGCCCTTGACGATTTGATGGGCATCTTTTTCCCTGGCAAAGAGGCTGATTCCGTGGTGAGGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCGGTTTCt                                                                                                      <  1:241131/147‑1 (MQ=255)                                                   ttttACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                    >  2:797210/1‑149 (MQ=255)                                                   ttttACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGAAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTa                                                                                                   <  2:833259/150‑1 (MQ=255)                                                     ttcccTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCTTTGACGTTTTGATTGGCACTTTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt                                                                                                 <  2:1295245/147‑1 (MQ=255)                                                     ttACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTc                                                                                                >  1:739362/1‑151 (MQ=255)                                                      tACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCATCTTTTTCCCGGGAAAAGAGGCTGATTTCGTGGTGAGGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTc                                                                                                <  2:897480/150‑1 (MQ=255)                                                        ccTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGTTTTGATTGGCAACTTTTTCCCTGGCAAAGGGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGAc                                                                                             <  2:316219/151‑1 (MQ=255)                                                             ccccGCTCGGGGGGGCAAAATCTCGGGCCCTTGTCTTTTTGATTGCCAACTTTTTCCCGGGCAAAGAGTCTGATTTCGTGGTGATGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAtt                                                                                        <  1:711173/148‑1 (MQ=255)                                                             ccACGCTCGGCGGAGCAAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAt                                                                                         <  2:645263/150‑1 (MQ=255)                                                                 gCTCGGGGGAGCAAATTCTCTGGGCTTTGAGGATTTGATTGGCAACTTTTTACCTGGCAAAGAGTCTGTTTTCGTGTTGAGGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt                                                                                     <  2:736807/150‑1 (MQ=255)                                                                      gcggAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgat                                                                               <  2:960770/151‑1 (MQ=255)                                                                             gAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACg                                                                        >  1:963642/1‑151 (MQ=255)                                                                               aaTCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTAcag                                                                                                                                 >  1:781645/1‑92 (MQ=255)                                                                               aaTCTCTGGCCCTTGACGATTTGTTGGGCAACTTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGTTGACGtt                                                                      <  1:353628/151‑1 (MQ=255)                                                                                aTCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTg                                                                     >  2:551755/1‑151 (MQ=255)                                                                                  ctctGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTggg                                                                   >  2:110594/1‑151 (MQ=255)                                                                                      gggCTTTGACGTTTTGTTTGGCAACTTTTTCCCTGGAAAAGAGGCTGATTTCGGGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGc                                                                  <  1:336493/148‑1 (MQ=255)                                                                                      gggCTTTGACGATTTGATGGGCAATTTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCg                                                                 <  1:1087974/149‑1 (MQ=255)                                                                                      gggCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGAt                                                               <  2:444041/151‑1 (MQ=255)                                                                                      ggcccTTGCCGTTTGGATGGGAAATTTTTTCCCTGCAAAAGGGGCTGATTTCGTGGTGATGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGa                                                                <  1:969974/147‑1 (MQ=255)                                                                                        gCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATg                                                              >  2:583716/1‑150 (MQ=255)                                                                                           ttGGCGATTTGATGGGCAATTTTTTCCCGGGAAAAGAGGTTGTTTTGGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGAc                                                            <  1:1086194/149‑1 (MQ=255)                                                                                             gACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGtt                                                        >  2:582432/1‑151 (MQ=255)                                                                                              aCGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTc                                                       >  2:793095/1‑151 (MQ=255)                                                                                               cGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGGGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTc                                                       <  1:223178/150‑1 (MQ=255)                                                                                                  tttGATTGGCATCTTTTTCCCGGGAAAAGAGGTTGATTTCGGGGTGATGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGtt                                                        <  1:892326/146‑1 (MQ=255)                                                                                                  tttGATTGGCAACTTTTTCCCGGGCAAAGAGGCTGATTTCGTGGTGATGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGAt                                                    <  2:816559/150‑1 (MQ=255)                                                                                                        tGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGc                                             >  1:20928/1‑151 (MQ=255)                                                                                                           cAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCAcc                                          >  2:648724/1‑151 (MQ=255)                                                                                                              cTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTAc                                       >  1:1281568/1‑151 (MQ=255)                                                                                                               tttttACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACg                                      >  1:355685/1‑151 (MQ=255)                                                                                                                   tACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTg                                   >  1:810363/1‑150 (MQ=255)                                                                                                                   tACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTg                                   >  1:647350/1‑150 (MQ=255)                                                                                                                   tACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTg                                   >  1:574009/1‑150 (MQ=255)                                                                                                                       gggCAAAGAGGCTTATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgt                              <  2:770870/150‑1 (MQ=255)                                                                                                                        ggaaaaGAGGTTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATgg                               <  1:1066943/146‑1 (MQ=255)                                                                                                                        ggCAAAGAGGTTGATTTCGGGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtc                             <  1:371962/151‑1 (MQ=255)                                                                                                                        ggCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtc                             >  2:169339/1‑151 (MQ=255)                                                                                                                              gatgcttattgggtggggATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCGCCTACGTTGATGgtcg                            <  1:500810/134‑1 (MQ=255)                                                                                                                              gagGCTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTgg                       <  2:81535/151‑1 (MQ=255)                                                                                                                              gagGCTGATTTCGTGGTGATGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTg                        <  1:734337/150‑1 (MQ=255)                                                                                                                                ggcggttttCGTGGTGAGGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgt                      <  1:158062/144‑1 (MQ=255)                                                                                                                                ggCTGTTTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtg                     <  1:920909/151‑1 (MQ=255)                                                                                                                                   tGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTAc                  >  1:700563/1‑151 (MQ=255)                                                                                                                                    gATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACg                 >  1:1155558/1‑151 (MQ=255)                                                                                                                                     aTTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACg                 >  1:763057/1‑150 (MQ=255)                                                                                                                                          gtggtgATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACg             <  1:336423/149‑1 (MQ=255)                                                                                                                                          gtggtgATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGc            >  2:126142/1‑150 (MQ=255)                                                                                                                                           tggtgATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCa           <  1:1284956/150‑1 (MQ=255)                                                                                                                                            ggtgATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACg             <  1:645962/147‑1 (MQ=255)                                                                                                                                            ggtgATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAAc         >  2:240490/1‑151 (MQ=255)                                                                                                                                            ggtgATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAAc         >  2:726110/1‑151 (MQ=255)                                                                                                                                              tgatgaACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCaa          <  1:239870/148‑1 (MQ=255)                                                                                                                                               gATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACt        <  2:1232999/149‑1 (MQ=255)                                                                                                                                               gATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACta       >  1:567851/1‑150 (MQ=255)                                                                                                                                                aggaaaCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaat     <  1:773998/147‑1 (MQ=255)                                                                                                                                                aTGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaa      <  1:927203/150‑1 (MQ=255)                                                                                                                                                 tGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaat     <  1:116713/150‑1 (MQ=255)                                                                                                                                                  ggAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaata    >  1:1111320/1‑150 (MQ=255)                                                                                                                                                  ggAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaata    >  2:961071/1‑150 (MQ=255)                                                                                                                                                    aaCCCAACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCGCCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaataat  <  1:30074/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAAT  >  NC_000913/3026793‑3027089

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