BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,538,736	C→T	6.0%	G283D (GGC→GAC) ‡	cysM ←	cysteine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,538,736	0	C	T	6.0%	316.3 / 5.5	99	G283D (GGC→GAC) ‡	cysM	cysteine synthase B
Reads supporting (aligned to +/- strand): ref base C (50/43); new base T (2/4); total (52/47)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

TATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGA  >  NC_000913/2538589‑2538885

tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTTGAAAGGTAGCGATCGCCa                                                                                                                                                     >  1:1321231/1‑150 (MQ=255)

tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCa                                                                                                                                                     <  2:1052818/150‑1 (MQ=255)

tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCa                                                                                                                                                     <  2:1196446/150‑1 (MQ=255)

tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCa                                                                                                                                                     >  2:1329119/1‑150 (MQ=255)

tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCa                                                                                                                                                     >  1:213518/1‑150 (MQ=255)

tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCCCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCa                                                                                                                                                     >  2:335804/1‑150 (MQ=255)

 atTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGCTCGCCAc                                                                                                                                                    >  1:451632/1‑150 (MQ=255)

  ttGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACg                                                                                                                                                   <  2:1022578/150‑1 (MQ=255)

    gCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGAt                                                                                                                                                 >  1:734641/1‑150 (MQ=255)

    gCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGAt                                                                                                                                                 <  1:1329119/150‑1 (MQ=255)

    gCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCAGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGAt                                                                                                                                                 >  2:173219/1‑150 (MQ=255)

     cAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATc                                                                                                                                                >  2:1142524/1‑150 (MQ=255)

      aGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCg                                                                                                                                               <  2:979476/150‑1 (MQ=255)

       gAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGc                                                                                                                                              >  2:156464/1‑150 (MQ=255)

       gAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATACTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGc                                                                                                                                              <  1:822352/150‑1 (MQ=255)

         aaCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAg                                                                                                                                            >  2:1271347/1‑150 (MQ=255)

          aCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGTGGCGATCGctac                                                                                                                                           <  1:1361383/150‑4 (MQ=255)

          aCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGTGGCGATCGctac                                                                                                                                           <  1:1033514/150‑4 (MQ=255)

          aCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAga                                                                                                                                           <  1:1024314/150‑1 (MQ=255)

           ccTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGTGGCGATCGctacc                                                                                                                                          <  1:427258/150‑5 (MQ=255)

                tAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCg                                                                                                                                     <  2:929170/150‑1 (MQ=255)

                   gCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGcca                                                                                                                                  <  1:481718/150‑1 (MQ=255)

                     cTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccacc                                                                                                                                <  2:73391/150‑1 (MQ=255)

                       gATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGcccccac                                                                                                                              >  2:929351/1‑150 (MQ=255)

                       gATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccaccac                                                                                                                              >  1:555063/1‑150 (MQ=255)

                            gCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCTACGATCGCAGATGATCGCCACCACCACcg                                                                                                                         >  2:1026715/1‑150 (MQ=255)

                                  cgcATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAg                                                                                                                   <  1:1461274/150‑1 (MQ=255)

                                   gcATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAgg                                                                                                                  >  1:183859/1‑150 (MQ=255)

                                   gcATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAgg                                                                                                                  >  2:304765/1‑150 (MQ=255)

                                    cATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAggg                                                                                                                 >  2:234275/1‑150 (MQ=255)

                                    cATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAggg                                                                                                                 <  2:1532240/150‑1 (MQ=255)

                                         ggCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGc                                                                                                            >  1:105006/1‑150 (MQ=255)

                                            aGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCttt                                                                                                         <  1:366761/150‑1 (MQ=255)

                                              ttttGCGTTTGTCATCAGTCTCCGATGCTATTATTCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTg                                                                                                       >  2:1500/1‑150 (MQ=255)

                                               tttGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGc                                                                                                      <  1:1572782/150‑1 (MQ=255)

                                                         tcatcaGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGt                                                                                            <  2:1076166/150‑1 (MQ=255)

                                                          catcaGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTg                                                                                           >  2:555980/1‑150 (MQ=255)

                                                          catcaGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTg                                                                                           >  1:830363/1‑150 (MQ=255)

                                                              aGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTcc                                                                                       >  1:1523799/1‑150 (MQ=255)

                                                                tctcCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCgg                                                                                     <  2:818692/150‑1 (MQ=255)

                                                                tctcCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCgg                                                                                     <  2:285554/150‑1 (MQ=255)

                                                                tctcCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACACGCAGTGCTCCgg                                                                                     <  1:779264/150‑1 (MQ=255)

                                                                  tcCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCa                                                                                   <  2:1347568/150‑1 (MQ=255)

                                                                   ccGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCaa                                                                                  >  2:1277354/1‑150 (MQ=255)

                                                                    cGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACACCGGTAGAAAGGTAGCGATCGCCACCATCGCAGATGATCGCCCCCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTACGGCAAc                                                                                 >  1:1187890/1‑150 (MQ=255)

                                                                       tGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCa                                                                                   >  1:1621356/1‑145 (MQ=255)

                                                                       tGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCa                                                                                   <  2:1621356/145‑1 (MQ=255)

                                                                        gCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACcgcg                                                                             <  2:1367937/150‑1 (MQ=255)

                                                                          tATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGcgcc                                                                           >  1:109216/1‑150 (MQ=255)

                                                                           aTTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGcgccg                                                                          <  1:173219/150‑1 (MQ=255)

                                                                           aTTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGcgccg                                                                          >  2:1182389/1‑150 (MQ=255)

                                                                             tAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCGCGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCgccgcc                                                                        >  2:1208748/1‑150 (MQ=255)

                                                                             tAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCgccgcc                                                                        <  1:153767/150‑1 (MQ=255)

                                                                               aTCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCgg                                                                      <  1:313888/150‑1 (MQ=255)

                                                                                 ccTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAg                                                                    >  1:915346/1‑150 (MQ=255)

                                                                                 ccTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAg                                                                    >  1:1001602/1‑150 (MQ=255)

                                                                                  cTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGc                                                                   >  2:1126673/1‑150 (MQ=255)

                                                                                  cTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGc                                                                   <  1:1500/150‑1 (MQ=255)

                                                                                  cTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGc                                                                   >  2:574029/1‑150 (MQ=255)

                                                                                       aTCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACg                                                              >  2:1261383/1‑150 (MQ=255)

                                                                                        tCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGc                                                             >  1:1518808/1‑150 (MQ=255)

                                                                                          cccGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCCCCACCCCCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGccc                                                           >  1:754531/1‑149 (MQ=255)

                                                                                             gCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAg                                                        <  2:301176/150‑1 (MQ=255)

                                                                                             gCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAg                                                        <  1:1277354/150‑1 (MQ=255)

                                                                                               ccccTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGaa                                                      <  1:1142524/150‑1 (MQ=255)

                                                                                                 ccTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAata                                                    >  2:426150/1‑150 (MQ=255)

                                                                                                       tAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAttcctt                                              >  2:663312/1‑150 (MQ=255)

                                                                                                         aaaTGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGAGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAttccttcc                                            >  1:299557/1‑150 (MQ=255)

                                                                                                          aaTGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTccc                                           <  2:555063/150‑1 (MQ=255)

                                                                                                          aaTGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTccc                                           >  2:138184/1‑150 (MQ=255)

                                                                                                           aTGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCg                                          <  1:1317348/150‑1 (MQ=255)

                                                                                                           aTGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCg                                          <  2:105006/150‑1 (MQ=255)

                                                                                                            tGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGCGCTGACGCCACAGAATATTCCTTCCCGc                                         >  1:1554285/1‑150 (MQ=255)

                                                                                                            tGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGc                                         <  2:934956/150‑1 (MQ=255)

                                                                                                            tGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGc                                         <  2:1523799/150‑1 (MQ=255)

                                                                                                                ctTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCg                                     >  1:211669/1‑150 (MQ=255)

                                                                                                                  tCCCCAAACACCCCGGTAGAAAGGTTGCGATCGCCACTATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGcc                                   <  2:1187890/150‑1 (MQ=255)

                                                                                                                   ccccAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCa                                  >  1:823791/1‑150 (MQ=255)

                                                                                                                    cccAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAg                                 >  1:652753/1‑150 (MQ=255)

                                                                                                                            ccccGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCa                         <  2:734641/150‑1 (MQ=255)

                                                                                                                             cccGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCAt                        >  1:1011716/1‑150 (MQ=255)

                                                                                                                                ggTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGgt                     >  1:154293/1‑150 (MQ=255)

                                                                                                                                  tAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGgtgt                   <  2:1527147/150‑1 (MQ=255)

                                                                                                                                  tAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGgtgt                   >  1:216722/1‑150 (MQ=255)

                                                                                                                                      aaGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTc               <  1:550395/150‑1 (MQ=255)

                                                                                                                                       catctgCGATCGTGGCGATCGCAGATGATCGCCACCACCCCCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTcc              >  1:1390476/6‑150 (MQ=255)

                                                                                                                                        atctgCGATCGTGGCGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAttcc                                                <  1:815648/111‑1 (MQ=255)

                                                                                                                                        atctgCGATCGTGGCGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAttcc                                                >  2:815648/5‑115 (MQ=255)

                                                                                                                                         gTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGc            <  2:915346/150‑1 (MQ=255)

                                                                                                                                         gTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGc            >  2:1307817/1‑150 (MQ=255)

                                                                                                                                          tAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCa           >  2:141551/1‑150 (MQ=255)

                                                                                                                                          tAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCa           <  1:304765/150‑1 (MQ=255)

                                                                                                                                             cGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATcg        <  1:609157/150‑1 (MQ=255)

                                                                                                                                              gATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATcgc       <  1:156464/150‑1 (MQ=255)

                                                                                                                                              gATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATcgc       >  2:502747/1‑150 (MQ=255)

                                                                                                                                              gATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATcgc       >  1:1563254/1‑150 (MQ=255)

                                                                                                                                                 cGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCt    <  2:1001602/150‑1 (MQ=255)

                                                                                                                                                  gCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCtg   >  1:1038838/1‑150 (MQ=255)

                                                                                                                                                   ccACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCtga  <  2:183859/150‑1 (MQ=255)

                                                                                                                                                   ccACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCtga  <  1:234275/150‑1 (MQ=255)

TATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGA  >  NC_000913/2538589‑2538885

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: `‑`