BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,388,932	A→G	A52A (GCT→GCC)	b4161 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,388,932	0	A	G	99.5%	637.1 / NA	207	A52A (GCT→GCC)	b4161	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (82/123); total (84/123)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 218 total aligned reads displayed.

ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388842‑4389021

aCGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGc                                                                                         <  2:1159725/93‑1 (MQ=255)

 cGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCg                                                                                        >  1:4157472/1‑93 (MQ=255)

 cGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCg                                                                                        <  2:2861375/93‑1 (MQ=255)

  gCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGt                                                                                       <  2:238381/93‑1 (MQ=255)

  gCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGt                                                                                       <  1:255125/93‑1 (MQ=255)

  gCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGt                                                                                       >  2:776829/1‑93 (MQ=255)

  gCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGt                                                                                       <  1:2023878/93‑1 (MQ=255)

   ccAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTg                                                                                      <  2:4269833/93‑1 (MQ=255)

    cAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGc                                                                                     >  2:1880841/1‑93 (MQ=255)

    cAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGc                                                                                     >  2:3303026/1‑93 (MQ=255)

     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                    <  2:255936/93‑1 (MQ=255)

      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAt                                                                                   <  1:513875/93‑1 (MQ=255)

      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAt                                                                                   >  2:3158482/1‑93 (MQ=255)

      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACACCGGCGGGCaa                                                                                   >  1:4196494/1‑92 (MQ=255)

       aCTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAtc                                                                                  >  1:2801510/1‑92 (MQ=255)

       aCTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATa                                                                                  <  2:4243840/93‑1 (MQ=255)

        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAc                                                                                 <  1:3095279/93‑1 (MQ=255)

        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGGGCATAc                                                                                 >  2:2211229/1‑93 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGTTTCCACATCCGAGTGCATAcc                                                                                >  1:4361338/1‑93 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                                <  2:4116099/93‑1 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                                >  1:4277845/1‑93 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                                <  1:3780552/93‑1 (MQ=255)

           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaa                                                                              <  1:3021744/93‑1 (MQ=255)

           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaa                                                                              <  2:2731237/93‑1 (MQ=255)

           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaa                                                                              <  1:1298863/93‑1 (MQ=255)

            gcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaaa                                                                             <  1:2004383/93‑1 (MQ=255)

             cgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAg                                                                            >  2:812420/1‑93 (MQ=255)

             cgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAg                                                                            >  1:1623947/1‑93 (MQ=255)

              ggTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGc                                                                           <  2:2645090/93‑1 (MQ=255)

                  gCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCt                                                                       >  1:2529146/1‑93 (MQ=255)

                   ccGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtg                                                                      >  2:3893598/1‑93 (MQ=255)

                    cGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtga                                                                     <  2:3565413/93‑1 (MQ=255)

                    cGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtga                                                                     <  2:1647937/93‑1 (MQ=255)

                     ggTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgat                                                                    >  2:3984296/1‑93 (MQ=255)

                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   <  1:1465217/93‑1 (MQ=255)

                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   <  2:4469377/93‑1 (MQ=255)

                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   >  1:1924077/1‑93 (MQ=255)

                       ttACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatga                                                                  >  1:1040592/1‑93 (MQ=255)

                       ttACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatga                                                                  >  1:3658929/1‑93 (MQ=255)

                       ttACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatga                                                                  <  1:2665264/93‑1 (MQ=255)

                        tACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatgag                                                                 >  1:2623154/1‑92 (MQ=255)

                        tACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGAc                                                                 >  2:4544662/1‑93 (MQ=255)

                           cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGAt                                                              <  1:4533976/93‑1 (MQ=255)

                           cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGAt                                                              <  1:3457743/93‑1 (MQ=255)

                            aGCGAACGGGTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCCTCCGCGGATTTCACAACCGCGTGGCCACCCAAACCGCTGATGACGaaa                                                             >  2:138268/1‑91 (MQ=255)

                            aGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATa                                                             <  1:989274/93‑1 (MQ=255)

                              gcaaCGGATGGGACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAcc                                                           <  2:2669869/91‑1 (MQ=255)

                               gAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCt                                                          <  2:461845/93‑1 (MQ=255)

                               gAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCt                                                          >  1:3577155/1‑93 (MQ=255)

                                aaCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCtt                                                         >  1:4505940/1‑93 (MQ=255)

                                  cGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                       <  2:1263656/93‑1 (MQ=255)

                                  cGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                       <  2:870876/93‑1 (MQ=255)

                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatc                                                     >  1:3365132/1‑93 (MQ=255)

                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatc                                                     >  1:4575218/1‑93 (MQ=255)

                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatc                                                     <  2:3517182/93‑1 (MQ=255)

                                     aTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatca                                                    <  2:3766289/93‑1 (MQ=255)

                                     aTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatca                                                    <  1:1311587/93‑1 (MQ=255)

                                     aTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatca                                                    <  2:2171856/93‑1 (MQ=255)

                                        ggaCGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  1:961572/91‑1 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  2:1431133/93‑1 (MQ=255)

                                          aCGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTc                                               >  2:885895/1‑93 (MQ=255)

                                                 atatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAAc                                        >  2:3027524/1‑93 (MQ=255)

                                                  tatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACa                                       <  2:3019671/93‑1 (MQ=255)

                                                   atTGCAGCGGTGATCGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAg                                      <  2:4374639/93‑1 (MQ=255)

                                                    ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCCGCGTGCATACCAAAGCGGCTCATGACGATACCTTCATCAGGCTCGCCAAACAGa                                     >  1:765927/1‑93 (MQ=255)

                                                      gCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAtt                                   <  2:2331153/93‑1 (MQ=255)

                                                      gCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAtt                                   <  1:1847486/93‑1 (MQ=255)

                                                      gCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAtt                                   <  1:2132986/93‑1 (MQ=255)

                                                         gCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTTGCCAAAAAGATTgtc                                >  1:1640599/1‑93 (MQ=255)

                                                         gCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtc                                <  1:1133010/93‑1 (MQ=255)

                                                         gCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtc                                <  1:2196243/93‑1 (MQ=255)

                                                           ggTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgt                              <  2:1246282/93‑1 (MQ=255)

                                                               aaCGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTa                          >  1:1010945/1‑93 (MQ=255)

                                                                aCGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAg                         <  1:3877287/93‑1 (MQ=255)

                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        >  2:1926397/1‑93 (MQ=255)

                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        <  1:583948/93‑1 (MQ=255)

                                                                  gTCGCCATCGGCGGATTCCTCATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGGCGTAGTc                       >  2:4541573/1‑93 (MQ=255)

                                                                  gTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTc                       >  2:2832911/1‑93 (MQ=255)

                                                                   tCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCg                      >  1:4040878/1‑93 (MQ=255)

                                                                    cGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCgg                     <  1:1743869/93‑1 (MQ=255)

                                                                    cGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCgg                     >  1:3943960/1‑93 (MQ=255)

                                                                     gCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCggg                    >  2:177913/1‑93 (MQ=255)

                                                                       cATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCt                  <  1:3998222/93‑1 (MQ=255)

                                                                        aTCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCtt                 >  1:839388/1‑93 (MQ=255)

                                                                        aTCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCtt                 <  2:3597882/93‑1 (MQ=255)

                                                                        aTCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCtt                 >  2:1060676/1‑93 (MQ=255)

                                                                            gcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTAGTCGTAGTCGGGCTTCTcc             >  2:1972288/1‑93 (MQ=255)

                                                                              ggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCtt           >  2:712601/1‑93 (MQ=255)

                                                                               gATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGGTTCTCCTTa          >  2:2764555/1‑93 (MQ=255)

                                                                               gATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTa          <  1:975787/93‑1 (MQ=255)

                                                                               gATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTa          >  2:2495330/1‑93 (MQ=255)

                                                                               gATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTa          <  2:1696234/93‑1 (MQ=255)

                                                                                 ttCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGa        >  1:871212/1‑93 (MQ=255)

                                                                                  tCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGAc       <  1:3157930/93‑1 (MQ=255)

                                                                                  tCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGAc       >  1:2338412/1‑93 (MQ=255)

                                                                                  tCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGAc       <  1:1392943/93‑1 (MQ=255)

                                                                                  tCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGAc       >  2:3363468/1‑93 (MQ=255)

                                                                                    cacaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGt     <  2:2279089/93‑1 (MQ=255)

                                                                                       aTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtttt  <  1:4386267/93‑1 (MQ=255)

                                                                                       aTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtttt  >  2:3438657/1‑93 (MQ=255)

ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388842‑4389021

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