BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,919,279	A→G	5.1%	V590V (GTA→GTG)	yebT →	MCE domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,919,279	0	A	G	5.1%	795.4 / 20.8	235	V590V (GTA→GTG)	yebT	MCE domain protein
Reads supporting (aligned to +/- strand): ref base A (130/92); new base G (8/4); total (138/97)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 235 total aligned reads displayed.

GTGAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTC  >  NC_000913/1919131‑1919427

gTGAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTAc                                                                                                                                                     >  2:2032343/1‑150 (MQ=255)

gTGAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTAc                                                                                                                                                     >  2:2003281/1‑150 (MQ=255)

  gAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAg                                                                                                                                                   <  2:104627/150‑1 (MQ=255)

   aaGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAgg                                                                                                                                                  >  1:3036030/1‑150 (MQ=255)

   aaGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAgg                                                                                                                                                  >  2:2861610/1‑150 (MQ=255)

   aaGTGATTACCGAGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAgg                                                                                                                                                  >  1:1350063/1‑150 (MQ=255)

      tGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCAt                                                                                                                                               <  2:196441/150‑1 (MQ=255)

       gATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATc                                                                                                                                              >  1:3190464/1‑150 (MQ=255)

       gATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATc                                                                                                                                              <  1:203330/150‑1 (MQ=255)

       gATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATc                                                                                                                                              >  1:4183725/1‑150 (MQ=255)

        aTTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATcc                                                                                                                                             >  2:4186582/1‑150 (MQ=255)

         ttACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATccc                                                                                                                                            >  1:168054/1‑150 (MQ=255)

           aCCGAGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCg                                                                                                                                          <  2:2993089/150‑1 (MQ=255)

            ccGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGc                                                                                                                                         >  1:1792113/1‑150 (MQ=255)

                gCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCtctc                                                                                                                                     <  1:89001/150‑1 (MQ=255)

                 cGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTcc                                                                                                                                    >  2:3368080/1‑150 (MQ=255)

                    ccGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCaga                                                                                                                                 <  1:3818930/150‑1 (MQ=255)

                    ccGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCaga                                                                                                                                 <  2:759317/150‑1 (MQ=255)

                     cgcgAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCagag                                                                                                                                >  2:573880/1‑150 (MQ=255)

                     cgcgAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCagag                                                                                                                                >  2:2198315/1‑150 (MQ=255)

                      gcgAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGTCCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGc                                                                                                                               <  1:3536446/150‑1 (MQ=255)

                      gcgAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGc                                                                                                                               >  1:256568/1‑150 (MQ=255)

                      gcgAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGc                                                                                                                               <  1:4005587/150‑1 (MQ=255)

                      gcgAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCACGCTCTCCAGAGc                                                                                                                               <  1:2662796/150‑1 (MQ=255)

                            tAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTaaa                                                                                                                         <  1:3597369/150‑1 (MQ=255)

                              aCGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAgg                                                                                                                       >  1:2475147/1‑150 (MQ=255)

                                gcgTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGa                                                                                                                     >  1:3251654/1‑150 (MQ=255)

                                    ttGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCa                                                                                                                 >  1:1399208/1‑150 (MQ=255)

                                        tatCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAg                                                                                                             >  2:2261423/1‑150 (MQ=255)

                                         atCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGc                                                                                                            <  1:697504/150‑1 (MQ=255)

                                         atCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGc                                                                                                            >  2:2848476/1‑150 (MQ=255)

                                            gATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTc                                                                                                         <  2:1277961/150‑1 (MQ=255)

                                            gATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTc                                                                                                         >  2:190554/1‑150 (MQ=255)

                                            gATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTc                                                                                                         <  1:1538486/150‑1 (MQ=255)

                                            gATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTc                                                                                                         <  2:2401048/150‑1 (MQ=255)

                                                tGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCTGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGTCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACa                                                                                                     >  1:54391/1‑150 (MQ=255)

                                                 gCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACac                                                                                                    >  1:3969155/1‑149 (MQ=255)

                                                 gCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACaa                                                                                                    >  1:3700256/1‑150 (MQ=255)

                                                  cATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAAc                                                                                                   >  2:3112363/1‑150 (MQ=255)

                                                      ttAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCa                                                                                               <  2:223739/150‑1 (MQ=255)

                                                          gCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCgg                                                                                           >  2:2687360/1‑150 (MQ=255)

                                                          gCCGGACTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCCGCGCGAAAGTACAGCTCACTGGTAGTCGTCTGACCGTACAGGCATCACCGCTCTCCAGAGCATTAAAGGGAGCCATAAGCTTCTACAACCTCAGCgc                                                                                           >  2:158959/1‑149 (MQ=255)

                                                           ccGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGt                                                                                          >  2:3522375/1‑150 (MQ=255)

                                                           ccGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGt                                                                                          <  2:2534127/150‑1 (MQ=255)

                                                            cGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTg                                                                                         >  2:382211/1‑150 (MQ=255)

                                                              gAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGcc                                                                                       >  1:3738024/1‑150 (MQ=255)

                                                                  aTCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAgcg                                                                                   <  2:2753657/150‑1 (MQ=255)

                                                                   tCGCAACCTTCTGACCAGCAAGAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAgcgc                                                                                  <  2:3322626/150‑1 (MQ=255)

                                                                     gCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCa                                                                                <  1:1746587/150‑1 (MQ=255)

                                                                      cAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAg                                                                               >  1:4082622/1‑150 (MQ=255)

                                                                      cAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAg                                                                               >  2:3239182/1‑150 (MQ=255)

                                                                      cAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGACCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAg                                                                               >  1:1135922/1‑150 (MQ=255)

                                                                       aaCCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGt                                                                              <  2:3567209/150‑1 (MQ=255)

                                                                       aaCCTTCGGACCAGCAAGAGCGTGTTCTGGGCAGCAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCAGCCCCGCTCTCCAGAGCATTACAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGt                                                                              <  2:1529562/150‑1 (MQ=255)

                                                                        aCCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTc                                                                             >  1:589528/1‑150 (MQ=255)

                                                                          cTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGTGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAg                                                                           >  2:1884964/1‑150 (MQ=255)

                                                                           ttCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGc                                                                          <  2:2778320/150‑1 (MQ=255)

                                                                           ttCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGc                                                                          >  1:1464448/1‑150 (MQ=255)

                                                                           ttCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGc                                                                          <  1:1001421/150‑1 (MQ=255)

                                                                                 ccAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAg                                                                    >  1:3859434/1‑150 (MQ=255)

                                                                                    gCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCg                                                                 >  1:2922874/1‑150 (MQ=255)

                                                                                             tgtTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTa                                                        <  2:1693774/150‑1 (MQ=255)

                                                                                               ttCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTAtt                                                      >  1:4170474/1‑150 (MQ=255)

                                                                                                  tGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTg                                                   >  1:1916550/1‑150 (MQ=255)

                                                                                                   gggCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGt                                                  <  1:2198315/150‑1 (MQ=255)

                                                                                                   gggCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGt                                                  <  1:2100497/150‑1 (MQ=255)

                                                                                                    ggCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTa                                                 <  1:3487099/150‑1 (MQ=255)

                                                                                                      cAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATg                                               >  1:2533337/1‑150 (MQ=255)

                                                                                                       aGAAGGCGGGGCGAAAGTTCAGCTGACTGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGc                                              <  2:1806360/150‑1 (MQ=255)

                                                                                                       aGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGc                                              <  1:3257589/150‑1 (MQ=255)

                                                                                                       aGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGc                                              >  1:4223085/1‑150 (MQ=255)

                                                                                                        gAAGGCGGGGCGAAAGTTCAGCTGACTGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCt                                             >  2:2939466/1‑150 (MQ=255)

                                                                                                         aaGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCtt                                            >  2:4175090/1‑150 (MQ=255)

                                                                                                            gCGGGGCGAAAGTTCAGCTGACTGGTAGTGGTCTGACCGGACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACACCCTCAGCGGTGCCCGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCg                                         <  2:4205514/150‑1 (MQ=255)

                                                                                                            gCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCAGTAAAGGGAGCCATTAGCTTCGACACCCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCg                                         >  1:2139552/1‑150 (MQ=255)

                                                                                                               gggCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGaaa                                      <  2:318736/150‑1 (MQ=255)

                                                                                                                 gCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACa                                    >  1:733905/1‑150 (MQ=255)

                                                                                                                  cGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAg                                   >  1:2668763/1‑150 (MQ=255)

                                                                                                                   gAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATTAGCGTATTCTGTATGCTTCCGAAACAGc                                  <  1:2801750/150‑1 (MQ=255)

                                                                                                                    aaaGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCg                                 <  1:1446193/150‑1 (MQ=255)

                                                                                                                      aGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGc                               <  1:2774768/150‑1 (MQ=255)

                                                                                                                        ttaaGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGccc                             >  1:391721/4‑150 (MQ=255)

                                                                                                                             cTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCg                        >  2:1673517/1‑150 (MQ=255)

                                                                                                                              tGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCgg                       >  1:1420839/1‑150 (MQ=255)

                                                                                                                              tGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCgg                       >  2:3208475/1‑150 (MQ=255)

                                                                                                                              tGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCgg                       >  1:1220195/1‑150 (MQ=255)

                                                                                                                                    gTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCg                 >  1:3065163/1‑150 (MQ=255)

                                                                                                                                     tAGTGGTCTGACCGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTTTGCTTCCGAAACAGCGGCCCGTGCGGTTGGCgg                >  1:3300290/1‑150 (MQ=255)

                                                                                                                                      aGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCggg               <  1:4035033/150‑1 (MQ=255)

                                                                                                                                        tGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCa             >  1:1736614/1‑150 (MQ=255)

                                                                                                                                          gTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGa           <  1:3984016/150‑1 (MQ=255)

                                                                                                                                           tCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGAt          >  2:1445630/1‑150 (MQ=255)

                                                                                                                                            cTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGAtt         >  2:2442582/1‑150 (MQ=255)

                                                                                                                                             tGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTa        >  1:326369/1‑150 (MQ=255)

                                                                                                                                             tGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTa        <  2:3924461/150‑1 (MQ=255)

                                                                                                                                                ccGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGc     >  1:2079466/1‑150 (MQ=255)

                                                                                                                                                ccGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGc     <  1:2228257/150‑1 (MQ=255)

                                                                                                                                                ccGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGc     >  2:1593410/1‑150 (MQ=255)

                                                                                                                                                 cGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCt    >  1:1877809/1‑150 (MQ=255)

                                                                                                                                                 cGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCt    >  2:2333881/1‑150 (MQ=255)

GTGAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTC  >  NC_000913/1919131‑1919427

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