Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,941	0	A	C	39.7%	121.6 / 18.2	73	E392D (GAA→GAC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (34/10);  new base C (0/29);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 3.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATA  >  NC_000913/3026795‑3027087                                                                                                                                                   |                                                                                                                                                    tACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑ccc                                                                                                                                                 >  2:1090462/1‑150 (MQ=255)  aaaaaGTATTGCAATTACAAGGCTATCGCCTCTCGCCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGCGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCAc                                                                                                                                               >  1:627213/3‑151 (MQ=255)  aCAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCAc                                                                                                                                               >  2:2059/1‑151 (MQ=255)            gCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCg                                                                                                                                     >  2:263352/1‑151 (MQ=255)             cAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGc                                                                                                                                    >  2:434530/1‑151 (MQ=255)              aaTGACAAAGCTATTGAATCTCGGCATATGAAGCGTTTTACCACGTCACGCTCGGCGGAGACACATCCCTGGGCCTTGACGCTTTGATTGGCAACTTTTTACCTGCAAAAGAGGCTGATTACGTCGTGATGGAA‑CCCACCGCCACTCCg                                                                                                                                     >  1:25846/1‑149 (MQ=255)                   cAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTAca                                                                                                                                >  1:1010069/1‑149 (MQ=255)                      ggCTATCGCCTCTCGGCATATGAAGCGTTTTCCCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCACCTTTTTACCGGGCAAAGAGGCTGATTTCGTGGTGATGGAC‑CCCACCGCCACTCCGCTAcagcag                                                                                                                            <  2:263597/150‑1 (MQ=255)                          aTCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTgcg                                                                                                                       >  2:1151734/1‑151 (MQ=255)                               ctctCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATg                                                                                                                  >  2:110223/1‑151 (MQ=255)                               ctctCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATg                                                                                                                  >  2:727210/1‑151 (MQ=255)                                  tcGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGAc                                                                                                                >  1:589235/1‑150 (MQ=255)                                   cGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACa                                                                                                               >  2:327993/1‑150 (MQ=255)                                    ggCATATGAGGTGTTTTCCCTGCCCCCGCTCGGCGGAGCGAAATCTCTGGCCTTTGCCGTTTTGTGTGGCACTTTTTTCCCGGGAAAAGGGGCTGATTGCGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGAc                                                                                                                <  2:559182/148‑1 (MQ=255)                                    ggCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAATTTTTTCCCGGGCAAAGAGGCTGTTTTCGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGa                                                                                                                 <  2:592208/147‑1 (MQ=255)                                    ggCATAAGAAGGTTTTTCCCTGGCCACCGTCGGCGGAGCGAAATCTCTGGGCCTGGACGATTTGATTGGCAACTTTTTCCCGGGAAAAGAGGCTGATTTCGTGGTGAGGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                             <  2:905649/151‑1 (MQ=255)                                    gacATATGACGCTTTTTCCCGGCCCCCGCTCGGCGGAGCGAAATCTCGGGGCCTTGCCGATTTGATTGGCAACTTTTTCCCGGGCAAAGAGGCTGTTTGCGTGGTGTTGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGa                                                                                                                 <  1:317636/145‑1 (MQ=255)                                     gCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAct                                                                                                            >  2:200230/1‑151 (MQ=255)                                       atatGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCCACTTTTTACCTGGAAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCAATCCGCTACAGCAGCTGCTCTATGACAAatct                                                                                                          >  1:305926/1‑151 (MQ=255)                                            aaGCTTTTTACCTGGCCACGCTGGGCGGAGCAAAATCTCGGGGCTTTGACGATTTGATTGGCAATTTTTTCCCTGGAAAAGGGGCTGTTTTGGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTc                                                                                                     <  2:981375/151‑1 (MQ=255)                                                 ttttCCCTGGCCACGCTCGGCGGAGCAAAATCGGTGGGCCTGGACGATTTTATTGCCAACTTTTGCCCGGGCAAAGAGGCTGATTTCGTGGTGAGGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                  <  1:905157/149‑1 (MQ=255)                                                 ttttACCTGCCCGCGCTGGGCGGGGCGAAATCTGGGGGCCTGGACGATTTGAGGGGCAACTTTTTACCGGGAAAAGGGGGGGTTTTGGGGGGAATGAAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTc                                                                                                     <  1:234024/146‑1 (MQ=255)                                                 ttttACCGGCCCCCGCTCGGCGGAGCGAATTCTCTGGGCCTTGACGTTTTGATTGGCACTTTTTTCCCGGGAAAAGAGGCTGATTTCGTGGTGAGGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCt                                                                                                    <  2:340241/147‑1 (MQ=255)                                                  tttACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt                                                                                               >  2:964751/1‑151 (MQ=255)                                                   taaCCTGGCCAGGCCGGGGGGGGGAAAACCTCTGGGCCTTGAAGTTTTGGTTGGCACTTTTTTCCCTGCAAAAAGGGTTGATTGGGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt                                                                                               <  2:355266/148‑1 (MQ=255)                                                      cccggCCCCGCCCGGCGAGGCAAAATCTCGGGGTCTTGACGATTTGTTGGGCATTTTTCTACCGGGAAAAGAGGCTGATTTGGGGGTGTTGAAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTc                                                                                              <  1:769465/145‑1 (MQ=255)                                                      ccTGGCCCCGCTCGGCGGAGCGAAATCTCGGGGCTTTGACGATTTGATTGGCAATTTTTTCCCGGGAAAAGAGGCTGTTTTCGTGGTGATGAAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTc                                                                                              <  2:589535/148‑1 (MQ=255)                                                      ccTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGAc                                                                                           >  2:107797/1‑151 (MQ=255)                                                      ccTGGCCACGCTCGGCGGAGCAAAATCTCTGGGCCTTGACGATTTGATTGGCACTTTTTTCCCTGGAAAAGGGGCTGATTTCGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCg                                                                                             <  2:282111/149‑1 (MQ=255)                                                         ggCCACGCTCGGCGGGGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACaa                                                                                         <  2:126598/150‑1 (MQ=255)                                                           ccccGCTGGGGGGAGCGAAATCTCTGGGCCTTGAGGATTTGTTTGGCAACTTTTTCCTGGGCAAGGGGGTTTATTTCGTGGTGATGGCC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAt                                                                                       <  1:1087794/147‑1 (MQ=255)                                                           ccACGCTCGGCGGGGCAAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGAAAGGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAt                                                                                       <  2:762703/150‑1 (MQ=255)                                                             aCGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGt                                                                                    >  1:1058887/1‑151 (MQ=255)                                                              cGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt                                                                                   >  1:604248/1‑151 (MQ=255)                                                                        aGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg                                                                         >  1:656777/1‑151 (MQ=255)                                                                         gaaaaaTCTCGGGCCCTTGACGATTTAATTGACAATTTTTCCCCGGAAAAAGAGGATGTTTTCGTGGTGAAGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt                                                                                   <  1:8487/137‑1 (MQ=255)                                                                            aaaTCTCTGGGCTTTGACGATTTATTTGGCATCTTTTAACCAGGGAAAGAGGCTGATTTCGGGGTGATGGAC‑CCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGAc                                                                       <  2:256924/149‑1 (MQ=255)                                                                            aaaTCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGt                                                                     >  1:889372/1‑151 (MQ=255)                                                                            aaaTCTCTGGCCCTTGACGATTTAATTGGCAACTTCTTCCCGGGAAAAGAGGTTG‑TTTTGTTGTTATTAAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGAc                                                                       <  1:219362/149‑1 (MQ=255)                                                                            aaaTCTCGGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGAc                                                                       <  2:49685/149‑1 (MQ=255)                                                                                    gggCTTTGACGATTTGATTGTCACCTTTTCCCCTGGAAAAGAGGTTGTTTTCGTGTTGATGCAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGAt                                                             <  2:926560/151‑1 (MQ=255)                                                                                    gggCTTTGACGATTTGATGGGCAGTTTTTTCCCTGGAAAAGGGGCTGATTTCGTGGTGAGGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGa                                                              <  2:117517/150‑1 (MQ=255)                                                                                       ccTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGAc                                                          >  2:16924/1‑151 (MQ=255)                                                                                          tGACGATTTGATGGGAAATTTTTGACCGGGCAAAGAGGCTTATTTCGTGGTGATGAAA‑CCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGAcc                                                         <  1:153752/149‑1 (MQ=255)                                                                                             cgttttGATTGGCAATTTTTTCCCGGGAAAAGGGGCTGATTTCGTGGTGATGAAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTc                                                     <  2:484611/147‑1 (MQ=255)                                                                                                tttgtttgACACTTTTTTACCTGGCAAGGAGGCTGTTTTGGGGGTGTTGAAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGAt                                                  <  2:24481/150‑1 (MQ=255)                                                                                                           aCTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACACATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCt                                       >  1:107255/1‑150 (MQ=255)                                                                                                             tttttCCCTGGAAAAGAGGTTGTTTTCGGGGTGAGGGAC‑CCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCa                                          <  1:443259/145‑1 (MQ=255)                                                                                                             tttttACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACg                                    >  1:498927/1‑151 (MQ=255)                                                                                                                ttcccTGGAAAAGAGGCTGATTTCGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGtt                                  <  2:195309/147‑1 (MQ=255)                                                                                                                ttACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTg                                 >  1:1143215/1‑151 (MQ=255)                                                                                                                     tGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgt                            >  1:920020/1‑151 (MQ=255)                                                                                                                     tGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgt                            >  2:618624/1‑151 (MQ=255)                                                                                                                     gggaaaggagGCTAATTGCGGGGGGATGAAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgt                            <  1:503977/144‑1 (MQ=255)                                                                                                                      ggCAAAGAGGCTGATTTCGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtc                           <  2:623475/151‑1 (MQ=255)                                                                                                                        cAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgt                         >  2:1114931/1‑151 (MQ=255)                                                                                                                        cAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgt                         >  1:469159/1‑151 (MQ=255)                                                                                                                         aaaGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgtc                        >  2:122143/1‑151 (MQ=255)                                                                                                                         aaaGAGGCTGATTTCGTGGTGATGAAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgt                         <  2:1073510/150‑1 (MQ=255)                                                                                                                           agagGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTg                      >  2:903427/1‑151 (MQ=255)                                                                                                                            gaggtttttttCTTGGTGATGAAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTg                      <  1:906800/142‑1 (MQ=255)                                                                                                                              ggttgtttgcggggtgATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgt                    <  1:312733/138‑1 (MQ=255)                                                                                                                                cTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtgt                  >  1:143618/1‑150 (MQ=255)                                                                                                                                cTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtgt                  >  2:914298/1‑150 (MQ=255)                                                                                                                                 tGATTTCGTGGTGATGGAA‑CCCACCGCAACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGAGATGATGACGTTGGGCGATGACTGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTAc                >  1:663335/1‑151 (MQ=255)                                                                                                                                  taTTTCTTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTa                 <  1:1151914/148‑1 (MQ=255)                                                                                                                                     ttCGTGGTGATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGaa             <  2:657093/150‑1 (MQ=255)                                                                                                                                        gtggtgATGGAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAAc            <  1:1082297/148‑1 (MQ=255)                                                                                                                                          ggtgATGAAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACg           <  2:498647/147‑1 (MQ=255)                                                                                                                                            tgATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACt      >  2:855139/1‑150 (MQ=255)                                                                                                                                                ggAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaata  <  2:185326/150‑1 (MQ=255)                                                                                                                                                ggAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaat   <  2:584653/149‑1 (MQ=255)                                                                                                                                                ggAC‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaa    <  2:586108/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAA‑CCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATA  >  NC_000913/3026795‑3027087

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