BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,624,135	A→G	N100N (AAT→AAC)	b3485 ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,624,135	0	A	G	96.8%	475.2 / 2.1	154	N100N (AAT→AAC)	b3485	putative transporter
Reads supporting (aligned to +/- strand): ref base A (1/4); new base G (68/81); total (69/85)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 159 total aligned reads displayed.

AATGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATATTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGATA  >  NC_000913/3624046‑3624224

aaTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGttt                                                                                        <  1:1135265/93‑1 (MQ=255)

aaTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGttt                                                                                        >  1:1234155/1‑93 (MQ=255)

 aTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTa                                                                                       <  2:1327630/93‑1 (MQ=255)

 aTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTa                                                                                       <  1:1673903/93‑1 (MQ=255)

 aTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATATTTa                                                                                       <  1:3825596/93‑1 (MQ=255)

  tGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTAt                                                                                      >  2:2519089/1‑93 (MQ=255)

  tGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTAt                                                                                      <  1:1838233/93‑1 (MQ=255)

   gCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGCCCTCACGCTGAACATTAGGCGGAATGTTTATc                                                                                     >  1:1032899/1‑93 (MQ=255)

   gCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATc                                                                                     >  1:78922/1‑93 (MQ=255)

    ccTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATcg                                                                                    >  1:4598413/1‑93 (MQ=255)

    ccTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATcg                                                                                    <  1:97806/93‑1 (MQ=255)

     cTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATcgc                                                                                   >  2:148317/1‑93 (MQ=255)

     cTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATcgc                                                                                   >  1:996672/1‑93 (MQ=255)

       ggCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGa                                                                                 <  1:2837896/93‑1 (MQ=255)

        gCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGaa                                                                                <  1:100199/93‑1 (MQ=255)

        gCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGaa                                                                                <  2:2392901/93‑1 (MQ=255)

        gCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGaa                                                                                <  2:900177/93‑1 (MQ=255)

        gCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGaa                                                                                <  2:287660/93‑1 (MQ=255)

          tCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAgg                                                                              <  2:953038/93‑1 (MQ=255)

          tCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAgg                                                                              >  2:2779482/1‑93 (MQ=255)

          tCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAgg                                                                              <  1:3538863/93‑1 (MQ=255)

           cATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGt                                                                             <  2:523656/93‑1 (MQ=255)

            aTGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGCATGTTTATCGCGAAGGta                                                                            <  2:4643959/93‑1 (MQ=255)

            aTGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGta                                                                            <  2:619302/93‑1 (MQ=255)

             tGCGCGTGGCATCGATGTTCACCGGAATAGCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGtat                                                                           <  1:2081377/93‑1 (MQ=255)

             tGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGtat                                                                           <  1:2081378/93‑1 (MQ=255)

              gcgcGTGGCATCGACGTTCTCCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGtttc                                                                          <  2:4225962/93‑4 (MQ=255)

              gcgcGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGtata                                                                          <  2:3526639/93‑1 (MQ=255)

              gcgcGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGtata                                                                          <  1:3569364/93‑1 (MQ=255)

               cgcgTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAg                                                                         <  2:1251601/93‑1 (MQ=255)

                gcgTGGCATCGACGTCCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGc                                                                        <  2:3982233/93‑1 (MQ=255)

                  gTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCgc                                                                      >  1:1652654/1‑92 (MQ=255)

                  gTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCgc                                                                      >  1:315457/1‑92 (MQ=255)

                   tGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTc                                                                     <  1:2417334/93‑1 (MQ=255)

                   tGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGCATAGCGTc                                                                     >  2:711358/1‑93 (MQ=255)

                        tCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCg                                                                >  2:4070254/1‑93 (MQ=255)

                        tCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATAACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCg                                                                >  1:3284473/1‑93 (MQ=255)

                           aCGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTcc                                                             <  1:2885241/93‑1 (MQ=255)

                             gTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAg                                                           >  1:4038314/1‑93 (MQ=255)

                             gTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAg                                                           >  1:1454191/1‑93 (MQ=255)

                                cACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTcc                                                        <  2:1281958/93‑1 (MQ=255)

                                 aCCTGAATATCCGGCTGGCGTCCGGCGAGGACATCATGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCg                                                       <  1:420584/93‑1 (MQ=255)

                                  ccTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCgg                                                      >  1:697601/1‑93 (MQ=255)

                                   cTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGc                                                     >  2:3162476/1‑93 (MQ=255)

                                   cTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGc                                                     >  1:843907/1‑93 (MQ=255)

                                    tGAATATCCGGCTGGCGTCCGGCGAGGACATCCCGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCAGGca                                                    >  1:738044/1‑93 (MQ=255)

                                    tGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGca                                                    >  1:1791686/1‑93 (MQ=255)

                                    tGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGca                                                    >  2:4642383/1‑93 (MQ=255)

                                     gAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcat                                                   >  2:3404108/1‑93 (MQ=255)

                                      aaTATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGACGGTATAGCGTCCGGCGTCCAGTCCGGcatc                                                  >  1:2202369/1‑93 (MQ=255)

                                       atatCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatcc                                                 <  1:263999/93‑1 (MQ=255)

                                       atatCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatcc                                                 <  2:3133212/93‑1 (MQ=255)

                                        tatCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatcca                                                <  1:1686159/93‑1 (MQ=255)

                                        tatCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatcca                                                <  2:1709284/93‑1 (MQ=255)

                                         atCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatccat                                               >  1:3446817/1‑93 (MQ=255)

                                          tCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATATTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatccatc                                              >  2:2987791/1‑93 (MQ=255)

                                           ccGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCAtct                                             >  1:683758/1‑93 (MQ=255)

                                            cGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCAtctc                                            <  1:632566/93‑1 (MQ=255)

                                             ggCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCa                                           >  1:1190299/1‑93 (MQ=255)

                                               cTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATc                                         <  2:3785279/93‑1 (MQ=255)

                                               cTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATc                                         <  2:2271052/93‑1 (MQ=255)

                                                 ggCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATcgg                                       >  2:2148587/1‑93 (MQ=255)

                                                 ggCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATcgg                                       >  1:662199/1‑93 (MQ=255)

                                                   cGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATcggcg                                     >  2:627177/1‑93 (MQ=255)

                                                     gccGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGt                                   <  1:1705169/92‑1 (MQ=255)

                                                        ggCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGTGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGat                                >  2:3252780/1‑93 (MQ=255)

                                                        ggCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGat                                <  2:1452695/93‑1 (MQ=255)

                                                         gCGAGGACATCACGCTGAAAATTAGGCGGAATATTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatc                               <  1:2357650/93‑1 (MQ=255)

                                                            aGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatcatc                            <  2:3222744/93‑1 (MQ=255)

                                                            aGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatcatc                            <  2:2924540/93‑1 (MQ=255)

                                                             ggACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCAtct                           >  2:3326045/1‑93 (MQ=255)

                                                              gACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCAtctc                          <  2:845126/93‑1 (MQ=255)

                                                              gACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCAtctc                          <  2:4362133/93‑1 (MQ=255)

                                                              gACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCAtctc                          >  2:4246376/1‑93 (MQ=255)

                                                               aCATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTcc                         <  2:834085/93‑1 (MQ=255)

                                                               aCATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTcc                         <  2:1857142/93‑1 (MQ=255)

                                                                cATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCg                        <  2:4302643/93‑1 (MQ=255)

                                                                 aTCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCgg                       >  1:3174976/1‑93 (MQ=255)

                                                                   cACGCTTAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTg                     <  2:2258745/93‑1 (MQ=255)

                                                                   cACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTg                     >  1:3353696/1‑93 (MQ=255)

                                                                     cGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGc                   >  1:2792913/1‑93 (MQ=255)

                                                                       cTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCaa                 >  2:2542842/1‑93 (MQ=255)

                                                                       cTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCaa                 >  2:4539786/1‑93 (MQ=255)

                                                                        tGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCaaa                >  2:4688340/1‑93 (MQ=255)

                                                                        tGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCaaa                >  1:3108189/1‑93 (MQ=255)

                                                                          aaaaTTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCaaaaa              <  2:3030545/93‑1 (MQ=255)

                                                                           aaaTTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAAc             <  1:1072519/93‑1 (MQ=255)

                                                                             aTTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCa           <  1:487985/93‑1 (MQ=255)

                                                                              ttAGGCGGAATGTTTATCGCTAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCAc          <  2:2243480/93‑1 (MQ=255)

                                                                              ttAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCAc          >  2:2401525/1‑93 (MQ=255)

                                                                               tAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACg         <  2:2042514/93‑1 (MQ=255)

                                                                                 ggCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGa       >  2:2028405/1‑93 (MQ=255)

                                                                                 ggCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGa       <  1:1746279/93‑1 (MQ=255)

                                                                                 ggCGGAATGTTTATCGCAAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGa       <  1:4221960/93‑1 (MQ=255)

                                                                                  gCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGAc      >  2:232210/1‑93 (MQ=255)

                                                                                   cGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACg     <  2:897035/93‑1 (MQ=255)

                                                                                    ggAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGa    <  1:519442/93‑1 (MQ=255)

                                                                                    ggAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGa    >  1:550684/1‑93 (MQ=255)

                                                                                     gAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGAt   <  2:1747025/93‑1 (MQ=255)

                                                                                      aaTGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGATa  <  2:536196/93‑1 (MQ=255)

AATGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATATTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGATA  >  NC_000913/3624046‑3624224

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