BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,896,350	G→T	10.1%	A87E (GCA→GAA) ‡	adeP ←	adenine:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,896,350	0	G	T	10.1%	445.3 / 26.6	139	A87E (GCA→GAA) ‡	adeP	adenine:H(+) symporter
Reads supporting (aligned to +/- strand): ref base G (66/59); new base T (12/2); total (78/61)
Fisher's exact test for biased strand distribution p-value = 2.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

Only 100 of 142 total aligned reads displayed.

CATCCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGG  >  NC_000913/3896202‑3896500

catcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGc                                                                                                                                                       <  1:376831/150‑1 (MQ=255)

  tcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCaa                                                                                                                                                     >  1:944222/1‑150 (MQ=255)

     aGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAATAGCGAAGAACGCATTCAGGCCCATAGCGGGTGGCAGTGCAACTg                                                                                                                                                  <  2:93054/150‑1 (MQ=255)

        aGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCa                                                                                                                                               <  2:1043584/150‑1 (MQ=255)

          cGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAgg                                                                                                                                             <  1:291769/150‑1 (MQ=255)

           gAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTTGCACTGGCaccc                                                                                                                                            <  2:913119/150‑4 (MQ=255)

           gAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGt                                                                                                                                            <  2:1841215/150‑1 (MQ=255)

              cgcgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAg                                                                                                                                         >  1:2014746/1‑150 (MQ=255)

                cgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCa                                                                                                                                       <  2:446626/150‑1 (MQ=255)

                cgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCa                                                                                                                                       <  2:507932/150‑1 (MQ=255)

                  gAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGAACAACGACAAAAGCGAAGCCCGCAGTCAGGCCCAGAGCGCGTGCCAGGGCACCTGGCAGGTTAGCaaa                                                                                                                                     >  2:708430/1‑150 (MQ=255)

                    aaaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACa                                                                                                                                   >  2:197014/1‑150 (MQ=255)

                         gTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTccc                                                                                                                              >  1:1280696/1‑150 (MQ=255)

                         gTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTccc                                                                                                                              >  1:1679893/1‑150 (MQ=255)

                            aGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTNCNNCNACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcatc                                                                                                                           <  2:1206823/150‑1 (MQ=255)

                             gTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcatca                                                                                                                          <  1:514253/150‑1 (MQ=255)

                             gTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACAGGCAGGTTAGCAAACAGTCCcatca                                                                                                                          <  2:90252/150‑1 (MQ=255)

                               aaCAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCata                                                                                                                        <  2:1122796/150‑1 (MQ=255)

                                aCAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCataa                                                                                                                       >  1:595103/1‑150 (MQ=255)

                                 cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGCACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACCGTCCCATCataat                                                                                                                      >  2:1575967/1‑150 (MQ=255)

                                 cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCataat                                                                                                                      <  2:572352/150‑1 (MQ=255)

                                    cagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACt                                                                                                                   >  1:1422471/1‑150 (MQ=255)

                                        ccTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCg                                                                                                               <  2:790454/150‑1 (MQ=255)

                                             cgcgCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGc                                                                                                          <  1:1018697/150‑1 (MQ=255)

                                             cgcgCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGc                                                                                                          <  1:2088907/150‑1 (MQ=255)

                                                gcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCCGc                                                                                                       >  1:127987/1‑150 (MQ=255)

                                                  cccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGa                                                                                                     >  2:3291/1‑150 (MQ=255)

                                                  cccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCAACGCCTGTACAACGACAAAAGCGAACAACGCATTCAGGCCCAAATCGGTTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGa                                                                                                     <  2:2213182/150‑1 (MQ=255)

                                                    cAGAAGATTGCGCCCATCCCGACCTGCCACGGCAATCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTTCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGTGATc                                                                                                   <  2:2289615/150‑1 (MQ=255)

                                                     agaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGACCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCa                                                                                                  >  1:2130623/1‑150 (MQ=255)

                                                     agaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCa                                                                                                  <  2:1624187/150‑1 (MQ=255)

                                                        agaTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGAc                                                                                               >  2:621672/1‑150 (MQ=255)

                                                             gcgcCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTa                                                                                          >  1:1204623/1‑150 (MQ=255)

                                                              cgcCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAg                                                                                         >  2:2383614/1‑150 (MQ=255)

                                                              cgcCAATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTat                                                                                         >  2:1720509/1‑149 (MQ=255)

                                                                cccATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGCACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTg                                                                                       >  1:1985924/1‑150 (MQ=255)

                                                                cccATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTg                                                                                       >  2:1109235/1‑150 (MQ=255)

                                                                  cATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGAc                                                                                     <  1:506386/150‑1 (MQ=255)

                                                                    tCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGa                                                                                   >  1:332357/1‑150 (MQ=255)

                                                                     cccGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGaa                                                                                  >  2:341436/1‑150 (MQ=255)

                                                                      ccGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAg                                                                                 <  1:1848958/150‑1 (MQ=255)

                                                                         aCCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACg                                                                              <  2:841108/150‑1 (MQ=255)

                                                                          ccTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgg                                                                             >  2:2111143/1‑150 (MQ=255)

                                                                          ccTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgg                                                                             >  2:787640/1‑150 (MQ=255)

                                                                            tGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGgcg                                                                           >  1:1944206/1‑150 (MQ=255)

                                                                            tGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGgcg                                                                           >  1:949156/1‑150 (MQ=255)

                                                                             gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGgcgc                                                                          >  1:2302714/1‑150 (MQ=255)

                                                                             gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTG‑CGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCt                                                                         >  2:96776/1‑150 (MQ=255)

                                                                               cACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTg                                                                        >  1:1327397/1‑150 (MQ=255)

                                                                                aCGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTgg                                                                       <  1:2348297/150‑1 (MQ=255)

                                                                                aCGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTgg                                                                       <  2:866980/150‑1 (MQ=255)

                                                                                  ggCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTa                                                                     <  1:2486110/150‑1 (MQ=255)

                                                                                  ggCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTa                                                                     >  1:1900540/1‑150 (MQ=255)

                                                                                    cAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATc                                                                   >  2:2469725/1‑150 (MQ=255)

                                                                                       gCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCAt                                                                <  1:154842/150‑1 (MQ=255)

                                                                                        cccATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATg                                                               <  1:806780/150‑1 (MQ=255)

                                                                                             cGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGc                                                          <  1:1302260/150‑1 (MQ=255)

                                                                                                cTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGCATGGATACCAGcgcc                                                       <  1:1836051/150‑4 (MQ=255)

                                                                                                 tGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGCAGACGGCGCTGGTATCCATGCCAGCAACg                                                      <  2:1343549/150‑1 (MQ=255)

                                                                                                 tGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACg                                                      <  2:1422471/150‑1 (MQ=255)

                                                                                                  gTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGc                                                     >  2:888082/1‑150 (MQ=255)

                                                                                                    aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCTGCTATGAGACAGTTAGTGACGAAGACGGCGCTTGTATCCATGCCAGCAActcga                                                   <  2:779970/150‑5 (MQ=255)

                                                                                                    aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCa                                                   >  1:638369/1‑150 (MQ=255)

                                                                                                     cAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCaa                                                  >  1:2355024/1‑150 (MQ=255)

                                                                                                     cAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCaa                                                  >  2:1413934/1‑150 (MQ=255)

                                                                                                     cAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCaa                                                  >  2:1027274/1‑150 (MQ=255)

                                                                                                       aCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGa                                                >  2:267634/1‑150 (MQ=255)

                                                                                                          aCAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAAtt                                             <  1:724217/150‑1 (MQ=255)

                                                                                                             aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGc                                          >  1:2316981/1‑150 (MQ=255)

                                                                                                             aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGc                                          >  2:1208605/1‑150 (MQ=255)

                                                                                                              aaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCg                                         <  2:1439816/150‑1 (MQ=255)

                                                                                                              aaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCg                                         <  2:1208869/150‑1 (MQ=255)

                                                                                                               aGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCgg                                        <  2:1994845/150‑1 (MQ=255)

                                                                                                                gCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCggg                                       <  1:552474/150‑1 (MQ=255)

                                                                                                                     gaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCg                                         >  1:1990004/1‑143 (MQ=255)

                                                                                                                     gaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCg                                         <  2:1990004/143‑1 (MQ=255)

                                                                                                                      aaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCACACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCCGCAACGCCAAGAATTTGCGGGTTAACa                                 >  2:733606/1‑150 (MQ=255)

                                                                                                                        cGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaa                               >  1:1057967/1‑150 (MQ=255)

                                                                                                                         gCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaa                              <  2:1490859/150‑1 (MQ=255)

                                                                                                                            ttCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACg                           <  1:787640/150‑1 (MQ=255)

                                                                                                                              cAGGCCCATAGCGGGTGCCAGTGCAACTGGCGGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGAt                         <  2:2216746/150‑1 (MQ=255)

                                                                                                                                ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGt                       >  2:1523012/1‑150 (MQ=255)

                                                                                                                                ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGt                       >  2:1745792/1‑150 (MQ=255)

                                                                                                                                 gCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTa                      >  2:2219818/1‑150 (MQ=255)

                                                                                                                                   acATAGCGGGTTCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGACGAATGAAGCGATCATACAGGTAGTGACTAATACGGCGCTGGTATCCATGCCAGAAACGCCAAGAATTTGCGGGTTAAAAAAAACGATGTaaa                    <  2:1723458/149‑1 (MQ=255)

                                                                                                                                       ctaacctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                >  1:1957813/7‑150 (MQ=255)

                                                                                                                                       ctaacctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                >  2:2145970/7‑150 (MQ=255)

                                                                                                                                       ctaacctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                >  2:1349497/7‑150 (MQ=255)

                                                                                                                                           ggTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTc            >  1:459456/1‑150 (MQ=255)

                                                                                                                                           cctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGta            >  1:1520543/3‑149 (MQ=255)

                                                                                                                                           cctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTc            >  1:166141/3‑150 (MQ=255)

                                                                                                                                           cctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTc            >  1:73695/3‑150 (MQ=255)

                                                                                                                                           cctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTc            >  1:1237522/3‑150 (MQ=255)

                                                                                                                                           cctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATAAAAATACTGCCGAATGCAGCGATAAGACAGGTAGTGACGAAGACGGCGCTGGTATACCTGCCAGCAACGCCAACAAATTTCGGGTTAACAAAAACCATGTAAACCAACGTc            >  2:2285678/3‑150 (MQ=255)

                                                                                                                                             tGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAg          <  1:1109235/150‑1 (MQ=255)

                                                                                                                                               ccAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCTCTGGTATACATGCCAGCAACGCCAAGAATTTGCTGGTTAACAAAAACGATGTAAAACATCGTCAGGa        <  2:186180/150‑1 (MQ=255)

                                                                                                                                                 aGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGTCAGCAACGCCAAGAATTTGCGGGTT‑ACAAAAACGGTGTTAACCATGGTgagggagg     >  2:1811583/1‑150 (MQ=255)

                                                                                                                                                    gCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtg   >  1:1695982/1‑150 (MQ=255)

                                                                                                                                                    gCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtg   >  2:193477/1‑150 (MQ=255)

                                                                                                                                                    gCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTG‑CGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtgg  <  1:96776/150‑1 (MQ=255)

CATCCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGG  >  NC_000913/3896202‑3896500

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