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

CCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACT  >  NC_000913/3026793‑3027083                                                                                                                                                      |                                                                                                                                               ccTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAAcc                                                                                                                                              >  4:341114/1‑151 (MQ=255)    aCAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCa                                                                                                                                            <  2:557692/150‑1 (MQ=255)      aaaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCg                                                                                                                                         <  1:638499/151‑1 (MQ=255)       aaGAATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGc                                                                                                                                        >  2:2230/1‑151 (MQ=255)        aGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGcc                                                                                                                                       <  1:247244/151‑1 (MQ=255)        aGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGcc                                                                                                                                       >  1:95098/1‑151 (MQ=255)           aTTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCAc                                                                                                                                     >  2:316482/1‑150 (MQ=255)               cAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGc                                                                                                                                >  4:239114/1‑151 (MQ=255)                aaTTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCt                                                                                                                               >  2:159309/1‑151 (MQ=255)                      aaGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcagc                                                                                                                          >  1:533705/1‑150 (MQ=255)                       aGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcagc                                                                                                                          >  3:424537/1‑149 (MQ=255)                         gCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCt                                                                                                                      >  2:385971/1‑151 (MQ=255)                          cTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTg                                                                                                                     >  1:65530/1‑151 (MQ=255)                           tATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTgc                                                                                                                    >  3:24935/1‑151 (MQ=255)                            aTCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTgcg                                                                                                                   >  1:957057/1‑151 (MQ=255)                              cGCCTCTCGGCATATGAAGCGTTTTCCCTGGCCACGCTCGGCGGAGCAAAATCTCTGGGCTTT‑GACGTTTTGTTTGGCATTTTTTTCCCTGGCAAGGAGGCTGATTTCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTgc                                                                                                                    <  4:81104/148‑1 (MQ=255)                                ccTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTGGGCGGAGCGAAATCTCTGGGCTTT‑GACGTTTTGTTTGGCAACTTTTTACCTGCCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTAt                                                                                                               <  3:239817/151‑1 (MQ=255)                                    tcGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGAc                                                                                                            >  2:123228/1‑150 (MQ=255)                                     cGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACa                                                                                                           >  2:1086444/1‑150 (MQ=255)                                     cGGAATATGAAGCGTTTTCCCTGGCCCCGCGCGGCGGAGCGAAATCTCGGGGCCTT‑GACGTTTGGATTGGCAATTTTTTACCGGGAAAAGAGGTTGTTTTCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACa                                                                                                           <  1:981750/150‑1 (MQ=255)                                       gCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAct                                                                                                        >  3:227159/1‑151 (MQ=255)                                        cATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActc                                                                                                       >  2:1053154/1‑151 (MQ=255)                                          tatGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTg                                                                                                     >  1:845928/1‑151 (MQ=255)                                           ataaaGCTTTTTCCCTGCCCCCGCCCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGCCACTTTTTTACCTGGAAAAGAGGCTAATTTCGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAct                                                                                                        <  1:281693/144‑1 (MQ=255)                                           atGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGt                                                                                                    >  2:139280/1‑151 (MQ=255)                                              aggttttttCCCGGCCCACGCTGGGGGGAGCAAAATCTCGGGGCCTG‑GACTTTTTGTTGGGCACTTTTTTCCCGGGAAAAGAGGTTGATTTTGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGttt                                                                                                  <  2:384390/145‑1 (MQ=255)                                              aaGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTc                                                                                                 >  2:34539/1‑151 (MQ=255)                                               aGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCt                                                                                                >  1:998438/1‑151 (MQ=255)                                                gCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCtt                                                                                               >  4:136885/1‑151 (MQ=255)                                                   ttttCCCGGGCCGCGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGTTTTGATTGCCAACTTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                              <  2:46477/149‑1 (MQ=255)                                                   ttttCCCGGGCCCCGCTCGGGGGAGCAAAATCTCTGGGCCTT‑GACGTTTTGATTGGCATCTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTa                                                                                             <  2:1039494/150‑1 (MQ=255)                                                    tttACTTGGCCCCGCCCGGGGAACCAAAATCTCTGGGCATT‑GACATTTTAATTGGAATTTTTTTCCCGGGAAAAAGGGCTATTTTTTGGGAAATGACACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAg                                                                                            <  2:352989/150‑1 (MQ=255)                                                    ttacccTGGCCACGCTGGGCGGAGCAAATTCTCGGGGCCTT‑GACGTTTTGATTGGCAATTTTTTCCCTGCAAAAGAGGCGGTTTTCGGGGTGATGGAACCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCtt                                                                                               <  2:1086274/143‑1 (MQ=255)                                                     ttACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTc                                                                                          >  2:882846/1‑151 (MQ=255)                                                        ccggcccccGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCATCTTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCg                                                                                         <  1:546861/141‑1 (MQ=255)                                                        ccTGGCCACGCTGGGCGGAGCGAAATCTCTGGGCTTT‑GACGATTTTATTGCCAATTTTTTACCGGGCAAAGGGGATGATTTGGGGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCg                                                                                         <  1:861781/149‑1 (MQ=255)                                                        ccTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGAc                                                                                       >  4:416941/1‑151 (MQ=255)                                                          gggCCACTCGCGGCGGAGGAAAACCCCGGGCCCTG‑GACGATTTGATTGTCATTTTTTTCCCGGGAAAAGGGGTTGTTTCCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCg                                                                                         <  1:557129/146‑1 (MQ=255)                                                           ggCCACGCTGGGGGGACGAAAATGTTGGGGCTTT‑GACGATTTGATGGAAAATTTTTGACCTGGAAAAGAGGTTGATTCCGTGGTAGTGAACCCCGCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                                                                                      <  1:823386/149‑1 (MQ=255)                                                             ccACCCTGGGCGGGGCGAAATCTCGGGGCCTT‑GACGTTTTGATTGCCACTTTTTCCCCTGGCAAAGGGGCTGTTTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                                                                                      <  2:261461/147‑1 (MQ=255)                                                                cccTCGGGGGCGCGAAACCTCGGGCCCTT‑GACGATTTGATTGCCATCTTTTTACCGGGAAAAGAGGCTTTTTTCGTGGTGATGCACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt                                                                               <  2:170112/149‑1 (MQ=255)                                                                cGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGt                                                                                <  4:113305/150‑1 (MQ=255)                                                                 gCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc                                                                              >  2:980537/1‑151 (MQ=255)                                                                     gggggAGCGAATTCTCGGGGCTTT‑GAAGATTTGTTGGGAAATTTTTTCCCGGGAAAAGAGGTTGTTTTGGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTg                                                                                 <  2:844236/141‑1 (MQ=255)                                                                     ggcggcGCGAATTCTCGGGCCGTT‑GAGGGTTTGATTGCCACATTTTTACCGGGCAAAAGGGCTTATTCCGTGGTGGTGGACCCCGCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                                                                                      <  1:961111/139‑1 (MQ=255)                                                                     ggcggAGCAAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGAAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCg                                                                             <  4:227106/148‑1 (MQ=255)                                                                     ggcggAGCAAAATCTCGGGGCCTT‑GACGTTTTGATTGGCACCTTTTTCCCTGGAAAAGAGGCTTTTTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtga                                                                          <  2:257413/151‑1 (MQ=255)                                                                        ggAGCGAAATCTCGGGGCCTT‑GACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGTTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatg                                                                        <  3:274783/150‑1 (MQ=255)                                                                            caaaaTCTGTGGGCTTT‑GACGATTTGATGGGCAATTTTTTCCCTGGAAAAGAGGCTTATTTCGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatga                                                                    <  2:482442/148‑1 (MQ=255)                                                                              aaaTCTCTGGGCTTT‑GAGTATTTGATGGGCAACTTTTTCCCTGGAAAGGAGGGTTTTTTGGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGt                                                                 <  2:609774/151‑1 (MQ=255)                                                                                       ggCCTTAGACGTTTTGATGGGCACTTTTTTCCCTGGAAAAGGGGCTGTTTTCGGGGTGAGGCACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTgg                                                              <  2:1128657/146‑1 (MQ=255)                                                                                           tt‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCg                                                    >  1:254942/1‑151 (MQ=255)                                                                                           tt‑GACGATTTGATTGGAAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGAcc                                                     <  4:215596/150‑1 (MQ=255)                                                                                            t‑GACGATTTGATTGGCAACTTTTTACCTGGAAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCAACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGTTGATGACGTTGGGCGATGACCGt                                                   >  1:802198/1‑151 (MQ=255)                                                                                               aCGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTc                                                 >  2:421713/1‑151 (MQ=255)                                                                                                   tttGATTGGCAACTTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGa                                               <  4:308544/149‑1 (MQ=255)                                                                                                     tGATTGCCACCTTTTTCCCTGCCAATGGGGCTAATTGCTTGTGGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCt                                            <  2:56597/150‑1 (MQ=255)                                                                                                      gATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTAc                                          >  1:453475/1‑151 (MQ=255)                                                                                                          gccatttttttCCCTGGAAAAGAGGTTATTTTCGTGGTGATAGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTAc                                          <  3:209922/141‑1 (MQ=255)                                                                                                            catttttttCCCTGGAAAAGGGGCTGTTTTCGTGGTGATGCACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGc                                       <  2:206241/144‑1 (MQ=255)                                                                                                             aattttttACCTGGAAAAGAGGCTGTTTTGGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCAcc                                    <  1:626558/147‑1 (MQ=255)                                                                                                             aattttttACCTGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCAc                                     <  1:1149969/146‑1 (MQ=255)                                                                                                                tttttACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACg                                >  1:1207969/1‑151 (MQ=255)                                                                                                                      ccTGGAAAAGAGGCTGTTTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCATCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATg                          <  3:51483/151‑1 (MQ=255)                                                                                                                         ggCAAAGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATgg                         <  3:281059/149‑1 (MQ=255)                                                                                                                         ggCAAAGAGGCTGATTTCGTGTTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgt                        <  4:306883/150‑1 (MQ=255)                                                                                                                            aaggagGTTGTTTTCGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgt                     <  3:410201/147‑1 (MQ=255)                                                                                                                            aaaGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcg                      <  1:883021/149‑1 (MQ=255)                                                                                                                             aaGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCt                   >  4:266763/1‑151 (MQ=255)                                                                                                                                 ggCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTg                  <  1:688416/148‑1 (MQ=255)                                                                                                                                 ggCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtg               <  4:272393/151‑1 (MQ=255)                                                                                                                                   cTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtgt              >  1:988528/1‑150 (MQ=255)                                                                                                                                    tgttttCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgt                <  2:963698/144‑1 (MQ=255)                                                                                                                                    tGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTAc            >  2:178898/1‑151 (MQ=255)                                                                                                                                    tGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTAc            >  1:383665/1‑151 (MQ=255)                                                                                                                                      aTTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGTTGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGCGTACg           >  4:258581/1‑150 (MQ=255)                                                                                                                                       tttCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACg           <  1:1130483/149‑1 (MQ=255)                                                                                                                                       tttCGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtgt              <  2:845758/146‑1 (MQ=255)                                                                                                                                         tCGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAAc        <  1:612228/150‑1 (MQ=255)                                                                                                                                          ggtggtgATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAAc        <  1:510504/148‑1 (MQ=255)                                                                                                                                           gtggtgATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGc      <  2:722182/150‑1 (MQ=255)                                                                                                                                           ggggtgAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCa     <  2:953685/149‑1 (MQ=255)                                                                                                                                              gtgATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACt  <  1:159332/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                               CCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTT‑GACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACT  >  NC_000913/3026793‑3027083

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