BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,361,189	G→C	5.3%	Q76E (CAG→GAG)	rhmD ←	L‑rhamnonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,361,189	0	G	C	5.3%	350.4 / 7.4	108	Q76E (CAG→GAG)	rhmD	L‑rhamnonate dehydratase
Reads supporting (aligned to +/- strand): ref base G (49/52); new base C (3/3); total (52/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 110 total aligned reads displayed.

GAGCCGGAGTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCAC  >  NC_000913/2361041‑2361338

gagCCGGAGTAATACAGGGTGGCACTGAGCATATGATCGTGGATCATTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTgg                                                                                                                                                      <  2:1444391/150‑1 (MQ=255)

 agCCGGAGTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGc                                                                                                                                                     <  1:1068654/150‑1 (MQ=255)

 agCCGGAGTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTATCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGc                                                                                                                                                     >  2:1962093/1‑150 (MQ=255)

   ccGGAGTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCg                                                                                                                                                   <  2:1409102/150‑1 (MQ=255)

       aGTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGtttt                                                                                                                                               <  2:191670/150‑1 (MQ=255)

       aGTAATACAGGGTGGCACGGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCCAGACCGGATTCGCCGtttc                                                                                                                                               <  2:325612/150‑2 (MQ=255)

        gTAATACAGGGTGGCACTGAGCATTTGCTCGGGGAGCAGTTTGATATCACTGCCACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTc                                                                                                                                              <  1:727689/150‑1 (MQ=255)

        gTAATACAGGGTGGCACTGAGCATTTGATCGTGGCTCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTATCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCCAGACCGGATTCGCCGtttcg                                                                                                                                              <  1:181482/150‑3 (MQ=255)

        gTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTc                                                                                                                                              <  2:1425948/150‑1 (MQ=255)

        gTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGACGTTTTc                                                                                                                                              >  1:588960/1‑150 (MQ=255)

          aaTACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCCAGACCGGATTCGCCGtttcgac                                                                                                                                            <  1:1691704/150‑5 (MQ=255)

          aaTACAGGGTGGCACGGAGCCTTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTATTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTg                                                                                                                                            <  1:1065798/150‑1 (MQ=255)

           aTACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACTCATTTACCCTCAATGAAACGGTTTAGATGTTTTTCCACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGc                                                                                                                                           >  2:120112/1‑150 (MQ=255)

           aTACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGc                                                                                                                                           >  2:1090571/1‑150 (MQ=255)

            tACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCt                                                                                                                                          <  2:1487808/150‑1 (MQ=255)

              ccgggTGGCACTGAGCATTTGATAGTGGATCAGTTTGTTATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTc                                                                                                                                        <  1:946950/148‑1 (MQ=255)

              cAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTTCCCTCTTTGAAACGGTTAAGATGTTTTTCTACAATAACACAGCACATTTCACCTGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTc                                                                                                                                        >  1:285187/1‑150 (MQ=255)

              cAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTc                                                                                                                                        >  1:298317/1‑150 (MQ=255)

                 gtgggCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGAc                                                                                                                                     <  1:484301/147‑1 (MQ=255)

                    ggCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAGTAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTc                                                                                                                                  >  1:1071886/1‑150 (MQ=255)

                     gCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTca                                                                                                                                 <  1:624270/150‑1 (MQ=255)

                     gCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCGGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTca                                                                                                                                 <  1:1292261/150‑1 (MQ=255)

                           aGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTGTTCTGCTTCGACTTCAAcaacc                                                                                                                           >  2:411256/1‑150 (MQ=255)

                             cATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCAACTTCAAcaaccaa                                                                                                                         >  1:137704/1‑150 (MQ=255)

                              aTTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCAATTCACCGGCTCTCGAAACGGCGAATCAGGTCTGGCCGTTTTCTGCTTCGACTTCAAcaaccaac                                                                                                                        >  1:987697/1‑150 (MQ=255)

                                 tGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTGAACAACCAACGTg                                                                                                                     >  2:33147/1‑150 (MQ=255)

                                 tGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTg                                                                                                                     <  2:1404203/150‑1 (MQ=255)

                                 tGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTg                                                                                                                     <  1:674361/150‑1 (MQ=255)

                                   aTCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGcc                                                                                                                   <  1:1709486/150‑1 (MQ=255)

                                    tCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCa                                                                                                                  <  1:1620694/150‑1 (MQ=255)

                                     cGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCTGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTACGACTTCAACAACCAACGTGCCaa                                                                                                                 >  2:1569246/1‑150 (MQ=255)

                                      gTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCaat                                                                                                                >  1:1365118/1‑149 (MQ=255)

                                             aGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTa                                                                                                         <  1:1825696/150‑1 (MQ=255)

                                              gTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTaa                                                                                                        >  1:1982556/1‑150 (MQ=255)

                                                ttGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGTGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATg                                                                                                      <  1:192173/150‑1 (MQ=255)

                                                ttGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATg                                                                                                      <  2:1194544/150‑1 (MQ=255)

                                                ttGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATg                                                                                                      <  1:704187/150‑1 (MQ=255)

                                                 tGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGc                                                                                                     <  2:1706004/150‑1 (MQ=255)

                                                 tGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTATCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTATCTGGTTCGACTTCAACACCCAACGTGCCAAGAACGTTAATGc                                                                                                     >  2:964768/1‑150 (MQ=255)

                                                   atatCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCCACGTGCCAAGAACGTTAATGCCa                                                                                                   >  1:1331521/1‑150 (MQ=255)

                                                    tatCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCaa                                                                                                  <  2:1538997/150‑1 (MQ=255)

                                                    tatCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCACATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCaa                                                                                                  >  1:1482940/1‑150 (MQ=255)

                                                      tCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGTCAAAt                                                                                                >  2:1616845/1‑150 (MQ=255)

                                                        aCTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGa                                                                                              >  2:1100936/1‑150 (MQ=255)

                                                        aCTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGa                                                                                              >  2:1452304/1‑150 (MQ=255)

                                                         cTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACTGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGAc                                                                                             >  2:1250408/1‑150 (MQ=255)

                                                         cTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCACCGTGCCAAGAACGGTAATGCCAAATgcc                                                                                             >  2:463249/1‑148 (MQ=255)

                                                         cTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGAc                                                                                             >  2:1498541/1‑150 (MQ=255)

                                                          tGACACATTTACCCTCAATGAAACGGTTATGATGTGTTTCGACAATCAAACATGACATTTCACCGGCTTTCGAAACGTAGAATCGGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCATCGTGCCAAGAACGTTAATGCCAAATGACt                                                                                            <  1:129033/150‑1 (MQ=255)

                                                          tGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACt                                                                                            <  2:297900/150‑1 (MQ=255)

                                                           gACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTg                                                                                           <  2:834641/150‑1 (MQ=255)

                                                                aTTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTg                                                                                      <  1:1945803/150‑1 (MQ=255)

                                                                  ttACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACATCCTACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTTAc                                                                                    <  2:83277/150‑1 (MQ=255)

                                                                  ttACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGAc                                                                                    <  1:1946135/150‑1 (MQ=255)

                                                                  ttACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGAc                                                                                    <  1:1090571/150‑1 (MQ=255)

                                                                   tACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACt                                                                                   <  2:1355907/150‑1 (MQ=255)

                                                                   tACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACt                                                                                   <  2:298317/150‑1 (MQ=255)

                                                                   tACCCTCAATGAAACGGTTAAGATGTTTTTCGAAAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACt                                                                                   <  1:856799/150‑1 (MQ=255)

                                                                        tCAATGAAACGGTTAAGATGGTTTTCGACAATAACACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCg                                                                              <  2:376881/150‑1 (MQ=255)

                                                                           aTGAAACGGTTAATATGTTTTTCGACAATAAAACAGCCCATTTCACCGTCTGTCGAAACGTCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATTACTGGCGTGACTGCTCGTaa                                                                           <  1:20182/150‑1 (MQ=255)

                                                                           aTGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGAGTGACTGCTCGTaa                                                                           >  1:1974140/1‑150 (MQ=255)

                                                                            tGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAAt                                                                          >  1:1513067/1‑150 (MQ=255)

                                                                            tGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAAt                                                                          >  2:1251239/1‑150 (MQ=255)

                                                                             gAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGAGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGGTAATGCCAAATGACTGGCGTCACTGCTCGGAATc                                                                         >  1:954779/1‑150 (MQ=255)

                                                                               aaCGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATcgc                                                                       <  2:928500/150‑1 (MQ=255)

                                                                                 cGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCgg                                                                     <  1:850951/150‑1 (MQ=255)

                                                                                  ggTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGt                                                                    >  2:1455787/1‑150 (MQ=255)

                                                                                  ggTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGt                                                                    >  1:330622/1‑150 (MQ=255)

                                                                                  ggTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGt                                                                    >  1:954876/1‑150 (MQ=255)

                                                                                   gTTAAGATGTTTTTCGACAATAAAACAGCCCATTACACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTa                                                                   <  2:422176/150‑1 (MQ=255)

                                                                                             ttttCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATc                                                         >  2:1638406/1‑150 (MQ=255)

                                                                                                tCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCggg                                                      <  1:1452304/150‑1 (MQ=255)

                                                                                                 cGACAATAAAACAGCCCATTTCCCCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCACCGTGCCAAGAACGTTAATGCCAAATGCCTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCgggg                                                     >  2:784580/1‑150 (MQ=255)

                                                                                                 cGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGAAATCCCGGCATTTACTAATCgggg                                                     >  1:59843/1‑150 (MQ=255)

                                                                                                 cGACAATAAAACAGCCCATTTCACCGCCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCgggg                                                     <  1:1801875/150‑1 (MQ=255)

                                                                                                  gACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGCCTGCTCGTAATCGCGGTATTTACTCATCggggg                                                    >  1:1551951/1‑149 (MQ=255)

                                                                                                   aCAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGCACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTg                                                   <  2:1070144/150‑1 (MQ=255)

                                                                                                    cAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCCACAACCAACGTGCCAAGAACGTTAAAGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTgg                                                  >  2:339917/1‑150 (MQ=255)

                                                                                                    cAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTgg                                                  >  1:491714/1‑150 (MQ=255)

                                                                                                      aTAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCa                                                <  1:938017/150‑1 (MQ=255)

                                                                                                            cAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAAAGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGa                                          >  2:980769/1‑150 (MQ=255)

                                                                                                              gCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATc                                        >  2:222306/1‑150 (MQ=255)

                                                                                                               cccATTTCACCGGCTGTCGAAATGGAGAATCCGGTATTGCCGTTTTCTGCTTCGACTTCAACAACCAACTAGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCg                                       <  2:285187/150‑1 (MQ=255)

                                                                                                               cccATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCAACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCg                                       <  2:137704/150‑1 (MQ=255)

                                                                                                                  aTTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCt                                    <  2:693260/150‑1 (MQ=255)

                                                                                                                     tCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATc                                 >  1:1644432/1‑150 (MQ=255)

                                                                                                                        ccGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAg                              >  1:1430757/1‑150 (MQ=255)

                                                                                                                         cGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACCACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGAGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCATTATGATCGTCTATCAAGt                             >  1:1908112/1‑150 (MQ=255)

                                                                                                                          ggCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTg                            >  1:1686068/1‑150 (MQ=255)

                                                                                                                              gTCGAAACGGCGAATCCGGGCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTc                        <  2:801272/150‑1 (MQ=255)

                                                                                                                               agaaaaCGGCCAGCCCGGATTCGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGGCTATCCAGTGATTCg                       >  2:437487/4‑150 (MQ=255)

                                                                                                                                 gAAACGGCGAATCCGGTCTCGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAATTGACTGGCCGGACTGCTCGTAATCGCGGTATTTACTCAACGGGGTGGCAATAAGATCGTCTATCCAGTGATTCGcc                     >  2:518447/1‑150 (MQ=255)

                                                                                                                                  aaaCGGCCAGACCGGATTCGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGGTAATGCCAAATGACAGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGAGCGTCTATCCAGTGATTCGccc                    >  2:1681686/1‑150 (MQ=255)

                                                                                                                                    aCGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGccccc                  >  1:1847591/1‑150 (MQ=255)

                                                                                                                                        cGAATCCGGTCTGGCGGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGt              <  1:92339/150‑1 (MQ=255)

                                                                                                                                        cGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGt              <  2:1071886/150‑1 (MQ=255)

                                                                                                                                          aaTCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCg            >  1:1474698/1‑150 (MQ=255)

                                                                                                                                             ccGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACATTCAAAGAACGTTAATGACAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTAATCGGGGTGGCAATATGATCTTCTATCCAGTGATTCGCCCCCTGGTCGTGa         >  1:95804/1‑150 (MQ=255)

                                                                                                                                              cGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGAt        <  2:1616039/150‑1 (MQ=255)

                                                                                                                                                    ggCCGTTTTCTCCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCCTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTCATTCGCCCCCTGGTCGAGATAAacac  >  2:1483399/1‑150 (MQ=255)

GAGCCGGAGTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCAC  >  NC_000913/2361041‑2361338

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