BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,749,776	T→C	A693A (GCA→GCG)	b1673 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,749,776	0	T	C	99.5%	667.4 / ‑4.9	214	A693A (GCA→GCG)	b1673	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (111/102); total (112/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

Only 100 of 226 total aligned reads displayed.

GCCTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCTGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTCTTT  >  NC_000913/1749687‑1749865

gCCTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTgccgcc                                                                                        >  2:480196/1‑93 (MQ=255)

gCCTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTgccgcc                                                                                        >  2:3844482/1‑93 (MQ=255)

gCCTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTgccgcc                                                                                        >  2:7394726/1‑93 (MQ=255)

gCCTCGTAAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTTCCGcc                                                                                        <  1:1628307/93‑1 (MQ=255)

 ccTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCa                                                                                       <  2:5520212/93‑1 (MQ=255)

 ccTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCa                                                                                       <  1:2643999/93‑1 (MQ=255)

  cTCGTCAACTTCGCTGACTTTCTCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAg                                                                                      <  1:7097525/93‑1 (MQ=255)

  cTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAg                                                                                      >  2:6966738/1‑93 (MQ=255)

   tCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGa                                                                                     <  1:7784225/93‑1 (MQ=255)

    tgTCAACTTCGCTGACTTTCGGTAACGGCAATTCATACTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGAt                                                                                    <  1:5042085/92‑1 (MQ=255)

    cGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGAt                                                                                    <  2:4056478/93‑1 (MQ=255)

    cGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGAt                                                                                    <  2:7264838/93‑1 (MQ=255)

     gTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATc                                                                                   <  1:6391872/93‑1 (MQ=255)

     gTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATc                                                                                   >  2:2755332/1‑93 (MQ=255)

      tCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcg                                                                                  >  1:880111/1‑93 (MQ=255)

      tCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcg                                                                                  <  2:1997503/93‑1 (MQ=255)

      tCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcg                                                                                  <  1:754386/93‑1 (MQ=255)

      tCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcg                                                                                  <  2:621472/93‑1 (MQ=255)

       cAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcgg                                                                                 <  1:6598149/93‑1 (MQ=255)

       cAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcgg                                                                                 >  1:5354365/1‑93 (MQ=255)

        aaCTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcggc                                                                                >  1:6783992/1‑93 (MQ=255)

        aaCTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcggc                                                                                <  1:4154982/93‑1 (MQ=255)

           ttCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCggcggc                                                                             >  2:619347/1‑93 (MQ=255)

             cGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCaa                                                                           <  2:6432173/93‑1 (MQ=255)

              gCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAat                                                                          <  1:5126276/93‑1 (MQ=255)

              gCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCAGCAat                                                                          >  1:3402000/1‑93 (MQ=255)

               cTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAata                                                                         >  1:1546195/1‑93 (MQ=255)

               cTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAata                                                                         >  2:4528678/1‑93 (MQ=255)

                 gACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATc                                                                       <  1:1783661/93‑1 (MQ=255)

                 gACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATc                                                                       >  2:4314225/1‑93 (MQ=255)

                 gACTTTCCCTAACGGCAATTCATCCTGATGGTTCCTCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAattac                                                                       >  2:2914506/1‑90 (MQ=255)

                  aCTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCCATATCt                                                                      >  1:2169322/1‑93 (MQ=255)

                  aCTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCt                                                                      >  1:5624186/1‑93 (MQ=255)

                   ctgtCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCtt                                                                     <  2:1704896/90‑1 (MQ=255)

                     ttCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTcc                                                                   <  1:5230490/93‑1 (MQ=255)

                     ttCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTcc                                                                   <  1:7722325/93‑1 (MQ=255)

                      tCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCa                                                                  <  1:2112267/93‑1 (MQ=255)

                      tCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCa                                                                  <  1:3834802/93‑1 (MQ=255)

                      tCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCa                                                                  <  1:2684051/93‑1 (MQ=255)

                        gCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGc                                                                >  1:774782/1‑93 (MQ=255)

                         cTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTTCAccc                                                               >  1:7225364/1‑93 (MQ=255)

                          tAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGccc                                                              >  2:1398131/1‑93 (MQ=255)

                            aCGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCa                                                            >  2:2676316/1‑93 (MQ=255)

                            aCGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCa                                                            <  2:4328485/93‑1 (MQ=255)

                            aCGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCa                                                            >  1:2492724/1‑93 (MQ=255)

                              ggCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGa                                                          <  1:6657641/93‑1 (MQ=255)

                               gCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGAc                                                         >  1:6099635/1‑93 (MQ=255)

                                 aaTTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGc                                                       >  1:2223037/1‑93 (MQ=255)

                                 aaTTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGc                                                       >  1:3917470/1‑93 (MQ=255)

                                   ttCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTg                                                     >  1:472835/1‑93 (MQ=255)

                                   ttCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTg                                                     <  2:2403175/93‑1 (MQ=255)

                                   ttCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTg                                                     <  1:5808193/93‑1 (MQ=255)

                                     cATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCa                                                   >  1:929738/1‑93 (MQ=255)

                                      aTCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCaa                                                  >  1:1199459/1‑93 (MQ=255)

                                        ccTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCTGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATg                                                >  2:5487766/1‑93 (MQ=255)

                                         cTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGt                                               <  2:4741773/93‑1 (MQ=255)

                                          tGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGtt                                              >  1:7346222/1‑93 (MQ=255)

                                          tGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGtt                                              >  1:3582357/1‑93 (MQ=255)

                                           gATGTTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGttt                                             <  1:3355949/93‑1 (MQ=255)

                                           gATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCACTATCTTCCAGCCCCAGACGCTGCAATGttt                                             >  1:3684091/1‑93 (MQ=255)

                                            aTGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTc                                            <  1:744614/93‑1 (MQ=255)

                                             tGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGGCCCAGACGCTGCAATGTTTcg                                           >  2:3185236/1‑93 (MQ=255)

                                             tGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTcg                                           >  2:3581673/1‑93 (MQ=255)

                                              ggTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTcgg                                          >  2:2105614/1‑92 (MQ=255)

                                                 tCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCggg                                       >  2:2272238/1‑93 (MQ=255)

                                                 tCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCggg                                       >  1:5045810/1‑93 (MQ=255)

                                                  cATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGt                                      >  2:3536741/1‑93 (MQ=255)

                                                   aTCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGtt                                     >  1:3070054/1‑93 (MQ=255)

                                                   aTCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGtt                                     <  2:5745835/93‑1 (MQ=255)

                                                    tCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTg                                    <  1:1883410/93‑1 (MQ=255)

                                                      ttACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGa                                  >  1:1985587/1‑93 (MQ=255)

                                                       tccccATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGAc                                 <  1:765802/91‑1 (MQ=255)

                                                       tACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGAc                                 >  1:836336/1‑93 (MQ=255)

                                                        aCCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACa                                >  1:412022/1‑93 (MQ=255)

                                                           cATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAAcc                             >  1:1042585/1‑93 (MQ=255)

                                                           cATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAAcc                             >  1:4674173/1‑93 (MQ=255)

                                                           cATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAAcc                             >  1:6338760/1‑93 (MQ=255)

                                                           cATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAAcc                             >  2:6563123/1‑93 (MQ=255)

                                                             ttCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCaa                           >  1:850689/1‑93 (MQ=255)

                                                              tCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCaac                          >  1:1815462/1‑92 (MQ=255)

                                                              tCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAg                          <  2:4751265/93‑1 (MQ=255)

                                                              tCCTTACGCAGGCAGTAGATTATGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAg                          >  1:5610541/1‑93 (MQ=255)

                                                               ccTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGc                         <  1:2896428/93‑1 (MQ=255)

                                                                cTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCt                        <  2:3448094/93‑1 (MQ=255)

                                                                cTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCt                        <  2:803271/93‑1 (MQ=255)

                                                                cTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCCCGGGTTGGAGAACCAAGCt                        >  2:1635144/1‑93 (MQ=255)

                                                                 ttACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCtg                       >  2:4588163/1‑93 (MQ=255)

                                                                  tACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCtgt                      <  1:2713975/93‑1 (MQ=255)

                                                                   aCGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCtgtg                     >  1:551552/1‑93 (MQ=255)

                                                                    cGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTgg                    <  2:7852819/93‑1 (MQ=255)

                                                                      caggcagTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGt                  <  2:4277702/93‑1 (MQ=255)

                                                                       aggcagTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTc                 >  2:1320744/1‑93 (MQ=255)

                                                                        ggcagTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTcc                >  2:61086/1‑93 (MQ=255)

                                                                           agTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAg             <  1:4320554/93‑1 (MQ=255)

                                                                            gTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGc            >  1:7187349/1‑93 (MQ=255)

                                                                             tAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGcc           >  1:3041202/1‑93 (MQ=255)

                                                                              aGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGccc          <  2:6284259/93‑1 (MQ=255)

                                                                                aTTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCAt        <  2:2841713/93‑1 (MQ=255)

                                                                                    ggtgCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTCt    <  2:5994193/92‑1 (MQ=255)

                                                                                      tgccgcCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTCttt  >  2:4450443/1‑93 (MQ=255)

GCCTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCTGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTCTTT  >  NC_000913/1749687‑1749865

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