BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,533,355	G→A	A176A (GCC→GCT)	yddE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,533,355	0	G	A	99.6%	929.5 / ‑5.4	255	A176A (GCC→GCT)	yddE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (113/141); total (113/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

Only 100 of 262 total aligned reads displayed.

GTTTTTGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGGGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATC  >  NC_000913/1533265‑1533444

gTTTTTGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGAGAAAGAGc                                                                                         >  1:4944633/1‑93 (MQ=255)

 tttttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCg                                                                                        >  2:5031286/1‑93 (MQ=255)

  ttttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGt                                                                                       <  2:9893878/93‑1 (MQ=255)

    ttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCg                                                                                     <  1:450243/93‑1 (MQ=255)

      gCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGat                                                                                   <  1:4178727/93‑1 (MQ=255)

        cTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGatat                                                                                 <  1:8402287/93‑1 (MQ=255)

        cTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGatat                                                                                 >  1:2238131/1‑93 (MQ=255)

        cTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGatat                                                                                 <  1:1670971/93‑1 (MQ=255)

         tGGGCGGATCGGGAACGGGAAGAAACCATTCCAGCAAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATc                                                                                <  1:7037405/93‑1 (MQ=255)

          gggCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATcc                                                                               >  2:923109/1‑93 (MQ=255)

          gggCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGAGTCGATATcc                                                                               >  1:2533018/1‑93 (MQ=255)

           ggCGGATCTGGAACGGGAAGAACCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCa                                                                              <  1:266773/93‑1 (MQ=255)

           ggCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCa                                                                              <  2:2610267/93‑1 (MQ=255)

            gCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCAc                                                                             <  1:1513085/93‑1 (MQ=255)

               gATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTc                                                                          <  1:5133862/93‑1 (MQ=255)

               gATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTc                                                                          >  1:7774225/1‑93 (MQ=255)

               gATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTc                                                                          <  1:348890/93‑1 (MQ=255)

               gATATGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTc                                                                          <  1:9815383/93‑1 (MQ=255)

                aTCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTcc                                                                         <  2:7185665/93‑1 (MQ=255)

                 tCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCg                                                                        >  2:8081692/1‑93 (MQ=255)

                 tCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCg                                                                        >  2:6090175/1‑93 (MQ=255)

                   tGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGt                                                                      >  2:3227994/1‑93 (MQ=255)

                   tGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGt                                                                      >  2:2564877/1‑93 (MQ=255)

                     gAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGttt                                                                    <  1:7851206/93‑1 (MQ=255)

                     gAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGttt                                                                    <  2:6971287/93‑1 (MQ=255)

                      aaCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTc                                                                   >  1:5519354/1‑93 (MQ=255)

                      aaCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTc                                                                   <  2:7915581/93‑1 (MQ=255)

                       aCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCa                                                                  >  2:1848535/1‑93 (MQ=255)

                       aCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCa                                                                  <  1:9146345/93‑1 (MQ=255)

                       aCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCa                                                                  <  1:5926533/93‑1 (MQ=255)

                         gggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGt                                                                >  2:3965981/1‑93 (MQ=255)

                          ggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTg                                                               >  1:6584136/1‑93 (MQ=255)

                          ggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTg                                                               <  2:459337/93‑1 (MQ=255)

                          ggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTg                                                               <  1:11328/93‑1 (MQ=255)

                           gaagaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTgg                                                              <  2:6240611/93‑1 (MQ=255)

                              gaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGat                                                           <  1:1275977/93‑1 (MQ=255)

                              gaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGat                                                           >  2:4528150/1‑93 (MQ=255)

                              gaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGat                                                           >  2:8390530/1‑93 (MQ=255)

                              gaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGat                                                           <  1:3689365/93‑1 (MQ=255)

                               aacccATGGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatc                                                          <  2:1699411/90‑1 (MQ=255)

                                acccATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCTCATTAAGGTCAGGCGAAAGAGCGTCGATATCCGCTTCCGGTTTCAGTGGGatca                                                         <  2:677829/91‑1 (MQ=255)

                                aaCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatca                                                         <  1:5810323/93‑1 (MQ=255)

                                 aCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatcat                                                        <  2:5844010/93‑1 (MQ=255)

                                  ccATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatcatc                                                       >  1:9363141/1‑93 (MQ=255)

                                   cATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGAtcatca                                                      <  1:472497/93‑1 (MQ=255)

                                   cATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGAtcatca                                                      >  2:5935692/1‑93 (MQ=255)

                                    aTTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCAc                                                     >  1:2528538/1‑93 (MQ=255)

                                    aTTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCAc                                                     >  1:3451807/1‑93 (MQ=255)

                                      tGCAGCCAATCGTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGTGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACtt                                                   <  2:2092696/93‑1 (MQ=255)

                                       gCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACttt                                                  <  1:8513383/93‑1 (MQ=255)

                                        cAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACtttt                                                 <  1:3899566/93‑1 (MQ=255)

                                         aGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTg                                                <  2:6287862/93‑1 (MQ=255)

                                         aGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTg                                                >  1:6194936/1‑93 (MQ=255)

                                         aGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTg                                                <  2:6637280/93‑1 (MQ=255)

                                         aGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTg                                                >  2:6858205/1‑93 (MQ=255)

                                         aGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTg                                                >  1:8144080/1‑93 (MQ=255)

                                          gCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGa                                               >  1:3021064/1‑93 (MQ=255)

                                          gCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGa                                               <  1:6098764/93‑1 (MQ=255)

                                           ccAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAg                                              >  1:9534737/1‑93 (MQ=255)

                                            cAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAgt                                             >  2:364861/1‑93 (MQ=255)

                                             aaTCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAgtg                                            <  2:1413905/93‑1 (MQ=255)

                                             aaTCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAgtg                                            <  2:9740508/93‑1 (MQ=255)

                                              aTCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAgtgt                                           >  2:1988877/1‑93 (MQ=255)

                                               tCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTc                                          <  2:7486351/93‑1 (MQ=255)

                                                cTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTcc                                         <  2:9603420/93‑1 (MQ=255)

                                                cTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTcc                                         <  2:7044053/93‑1 (MQ=255)

                                                 tttttGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCt                                        <  2:1181150/93‑1 (MQ=255)

                                                   tttGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGTGATCATCACTTTTGAGTGTCCtgt                                      <  1:1971475/93‑1 (MQ=255)

                                                    ttGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCtgtt                                     <  2:3732692/93‑1 (MQ=255)

                                                    ttGCTGATAGCGGTCAGCGCATTAAGGTCAGGCCAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCtgtt                                     <  1:687688/93‑1 (MQ=255)

                                                      gCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCtgttgt                                   <  1:8859891/93‑1 (MQ=255)

                                                            agaggTCAGCGCATTAAGGTCAGGCGAAGGAGCGTCGATATCCCCTCCCGGTGTCAGTGGGATCATCCCTTTTGAGTGTCCTGTTGTTGCCaa                             <  2:8241272/90‑2 (MQ=255)

                                                                ggcAGCGCAATAAAGGCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGa                         <  1:3137367/91‑1 (MQ=255)

                                                                gTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGa                         >  1:9597086/1‑93 (MQ=255)

                                                                gTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGa                         <  2:87501/93‑1 (MQ=255)

                                                                 tCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGaa                        <  1:2152214/93‑1 (MQ=255)

                                                                  cAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAAt                       <  2:5494495/93‑1 (MQ=255)

                                                                   aGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAAtt                      <  2:1492669/93‑1 (MQ=255)

                                                                   aGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAAtt                      <  2:6803881/93‑1 (MQ=255)

                                                                    gcgcATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTg                     >  1:7275326/1‑93 (MQ=255)

                                                                    gcgcATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTg                     <  2:5597030/93‑1 (MQ=255)

                                                                        aatAAGGGCAGGAGAAAGAGCGTCCATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCaa                 <  1:7185293/91‑1 (MQ=255)

                                                                         ttAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAg                <  2:1319729/93‑1 (MQ=255)

                                                                         ttAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAg                <  2:2298943/93‑1 (MQ=255)

                                                                           aaGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGcc              >  2:646191/1‑93 (MQ=255)

                                                                           aaGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGcc              <  2:1100662/93‑1 (MQ=255)

                                                                           aaGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGcc              <  2:3108229/93‑1 (MQ=255)

                                                                             ggTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTg            <  2:506151/93‑1 (MQ=255)

                                                                             ggTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTg            <  2:7667855/93‑1 (MQ=255)

                                                                             ggTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTg            >  2:47658/1‑93 (MQ=255)

                                                                               tCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGc          >  2:9151759/1‑93 (MQ=255)

                                                                                 aGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAg        <  2:9695856/93‑1 (MQ=255)

                                                                                     gAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAata    <  1:8461689/93‑1 (MQ=255)

                                                                                      aaaGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAatat   >  1:2748760/1‑93 (MQ=255)

                                                                                      aaaGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAatat   >  2:5350094/1‑93 (MQ=255)

                                                                                      aaaGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAatat   <  1:8168298/93‑1 (MQ=255)

                                                                                      aaaGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAatat   >  1:6409040/1‑93 (MQ=255)

                                                                                       aaGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATc  >  2:6386775/1‑93 (MQ=255)

                                                                                       aaGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATc  >  1:5236290/1‑93 (MQ=255)

                                                                                       aaGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATc  <  2:3565166/93‑1 (MQ=255)

GTTTTTGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGGGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATC  >  NC_000913/1533265‑1533444

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