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.3%	482.3 / ‑5.1	137	A176A (GCC→GCT)	yddE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (75/61); total (76/61)
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.56e-01

Only 100 of 138 total aligned reads displayed.

GTTTTTGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGGGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATC  >  NC_000913/1533265‑1533444

gTTTTTGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGc                                                                                         >  1:3288507/1‑93 (MQ=255)

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

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

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

  ttttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGt                                                                                       >  2:1199031/1‑93 (MQ=255)

   tttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTc                                                                                      <  1:2081784/93‑1 (MQ=255)

    ttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGGGCATTAAGGCCAGGCGAAAGAGGGTCg                                                                                     >  1:575396/1‑93 (MQ=255)

     tGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGa                                                                                    <  2:1168122/93‑1 (MQ=255)

     tGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGa                                                                                    <  2:4637343/93‑1 (MQ=255)

       ccTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGata                                                                                  <  1:1716172/93‑1 (MQ=255)

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

         tGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATc                                                                                >  2:71071/1‑93 (MQ=255)

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

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

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

              ggATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACtt                                                                           >  2:2039625/1‑93 (MQ=255)

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

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

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

                    ggAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGtt                                                                     <  2:2944890/93‑1 (MQ=255)

                    ggAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGtt                                                                     >  2:4296402/1‑93 (MQ=255)

                     gAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGttt                                                                    >  1:2283116/1‑93 (MQ=255)

                      aaCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGCAAGAGCGTCGATATCCACTTCCGGTTTc                                                                   >  2:2321199/1‑93 (MQ=255)

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

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

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

                        cGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAg                                                                 <  2:3190361/93‑1 (MQ=255)

                         gggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGt                                                                <  1:2640536/93‑1 (MQ=255)

                         gggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATACCCACTTCCGGTTTCAGt                                                                <  1:4640132/93‑1 (MQ=255)

                          ggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCCCTTCCGGTTTCAGTg                                                               >  2:1815677/1‑93 (MQ=255)

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

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

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

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

                           gaagaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTgg                                                              >  2:3779123/1‑93 (MQ=255)

                            aagaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTggg                                                             <  1:2025869/93‑1 (MQ=255)

                              taaaCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGat                                                           >  2:4694503/2‑93 (MQ=255)

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

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

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

                              gaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGAGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGat                                                           >  2:576174/1‑93 (MQ=255)

                               aaaCTATTGCAGGCAATGTTGTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGAGATCCACTTCCGGTTTCAGTGGGatc                                                          <  1:4465872/93‑1 (MQ=255)

                               aaaCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatc                                                          <  2:2800357/93‑1 (MQ=255)

                                 aCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatcat                                                        >  2:1789724/1‑93 (MQ=255)

                                  ccATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCTCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatcatc                                                       <  1:3272303/93‑1 (MQ=255)

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

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

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

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

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

                                    aTTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCAc                                                     <  1:4112229/93‑1 (MQ=255)

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

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

                                        cAGCCAATCTTGTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTCCCGGTTTCAGTGGGATCATCACtttt                                                 <  2:582499/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                               tCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTc                                          >  1:2675468/1‑93 (MQ=255)

                                                 tttttGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCt                                        >  2:4697472/1‑93 (MQ=255)

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

                                                   tttGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCtgt                                      >  2:3020967/1‑93 (MQ=255)

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

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

                                                     tGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCtgttg                                    >  1:2996312/1‑93 (MQ=255)

                                                        tGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGttgttg                                 <  2:1148299/93‑1 (MQ=255)

                                                        tGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGttgttg                                 >  2:2635959/1‑93 (MQ=255)

                                                        tGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGttgttg                                 >  1:1537034/1‑93 (MQ=255)

                                                         gATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGc                                >  2:3717630/1‑93 (MQ=255)

                                                         gATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGc                                >  2:4123229/1‑93 (MQ=255)

                                                         gATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGc                                >  1:974754/1‑93 (MQ=255)

                                                         gATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGc                                <  1:729477/93‑1 (MQ=255)

                                                             gCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCAcc                            >  1:2776925/1‑93 (MQ=255)

                                                             gCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCAcc                            <  2:2011870/93‑1 (MQ=255)

                                                             gCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCAcc                            >  1:1642773/1‑93 (MQ=255)

                                                             gCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCAcc                            <  1:1802409/93‑1 (MQ=255)

                                                               ggTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTg                          <  1:3676905/93‑1 (MQ=255)

                                                               ggTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTg                          <  2:4068332/93‑1 (MQ=255)

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

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

                                                                      gcATTAAGGTCAGGCGAAAGGGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGc                   >  2:2468971/1‑93 (MQ=255)

                                                                       cATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCa                  <  2:2219839/93‑1 (MQ=255)

                                                                        aTTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCaa                 <  1:2702848/93‑1 (MQ=255)

                                                                        aTTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCaa                 <  1:2209075/93‑1 (MQ=255)

                                                                          tAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGc               >  2:3346852/1‑93 (MQ=255)

                                                                            aGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCt             <  1:3667401/93‑1 (MQ=255)

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

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

                                                                                  ggCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGa       >  2:4488824/1‑93 (MQ=255)

                                                                                  ggCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGa       >  1:2842632/1‑93 (MQ=255)

                                                                                    cGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCCAGCCTGGCAGAat     >  2:635938/1‑93 (MQ=255)

                                                                                    cGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAat     <  1:2709045/93‑1 (MQ=255)

                                                                                    cGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAat     >  2:635937/1‑93 (MQ=255)

                                                                                     gAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAata    >  1:1996502/1‑93 (MQ=255)

                                                                                     gAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAata    >  2:2236208/1‑93 (MQ=255)

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

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

GTTTTTGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGGGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATC  >  NC_000913/1533265‑1533444

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