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.4%	1092.4 / NA	345	A52A (GCT→GCC)	b4161	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (164/179); total (166/179)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 377 total aligned reads displayed.

GACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388841‑4389021

gACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGGTTCCACATCAg                                                                                          >  2:326479/1‑93 (MQ=255)

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

    ccAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTg                                                                                      >  1:8491390/1‑93 (MQ=255)

    ccAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTg                                                                                      >  1:6469180/1‑93 (MQ=255)

      aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                    >  2:916551/1‑93 (MQ=255)

      aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                    >  2:6504867/1‑93 (MQ=255)

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

          tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGGATAcc                                                                                >  1:2339105/1‑93 (MQ=255)

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

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

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

           aCGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCa                                                                               >  1:9634888/1‑93 (MQ=255)

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

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

                  cGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGc                                                                        <  2:9063342/93‑1 (MQ=255)

                  cGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGc                                                                        <  1:912759/93‑1 (MQ=255)

                  cGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGc                                                                        <  1:3224793/93‑1 (MQ=255)

                  cGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGc                                                                        <  1:3554278/93‑1 (MQ=255)

                   gCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCt                                                                       <  2:4204153/93‑1 (MQ=255)

                   gCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAACAAAGCGGCt                                                                       >  2:1157496/1‑93 (MQ=255)

                    ccGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGGtg                                                                      >  1:3710933/1‑93 (MQ=255)

                      ggTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgcg                                                                    >  2:5884861/1‑91 (MQ=255)

                       gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCCAAGCGGCtgatg                                                                   >  2:3058202/1‑93 (MQ=255)

                       gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCCTGCATACCAAAGCGGctgctg                                                                   >  2:152506/1‑93 (MQ=255)

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

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

                           ccAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGa                                                               <  2:4029554/93‑1 (MQ=255)

                           ccAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGa                                                               <  2:8647158/93‑1 (MQ=255)

                            cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGAt                                                              >  2:6853430/1‑93 (MQ=255)

                               cGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAcc                                                           >  2:8873692/1‑93 (MQ=255)

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

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

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

                                 aaCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCGCATCGGCGTGCATACCAAAGCGGCTGATGACGATACCtt                                                         <  1:1942495/93‑1 (MQ=255)

                                 aaCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCTAAACGGCGGATGGCGATACCtt                                                         >  1:3483709/1‑93 (MQ=255)

                                 aaCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCtt                                                         <  2:4501052/93‑1 (MQ=255)

                                 aaCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCtt                                                         <  2:3662972/93‑1 (MQ=255)

                                  aCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCTGCGTGCATACCAAAACGGCTGATGGCGATACCTtc                                                        >  2:3481647/1‑93 (MQ=255)

                                   cGGGTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGGTGACGGTACGTtcc                                                       >  2:6685171/1‑92 (MQ=255)

                                    ggATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcat                                                      >  1:6630815/1‑93 (MQ=255)

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

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

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

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

                                       tGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                   >  1:9782853/1‑93 (MQ=255)

                                       tGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                   <  1:1450196/93‑1 (MQ=255)

                                        ggTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATCCCTTCATCAgg                                                  >  1:4766153/1‑93 (MQ=255)

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

                                            cGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCg                                              >  1:7175393/1‑93 (MQ=255)

                                            cGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGGCGATAACTTCATCAGGCTCg                                              >  1:4863503/1‑93 (MQ=255)

                                             ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGc                                             >  1:1236410/1‑93 (MQ=255)

                                              gCGAATATTGCAGCGGTTAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGcc                                            >  1:8611348/1‑93 (MQ=255)

                                              gCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGcc                                            >  1:9398562/1‑93 (MQ=255)

                                              gCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGcc                                            >  1:4680315/1‑93 (MQ=255)

                                                gAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCaa                                          >  1:785128/1‑93 (MQ=255)

                                                gAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCACGCTCGCCaa                                          >  2:8257415/1‑93 (MQ=255)

                                                   tatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACa                                       >  2:2195082/1‑93 (MQ=255)

                                                   tatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACa                                       >  2:6964338/1‑93 (MQ=255)

                                                    atTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAg                                      >  1:8915116/1‑93 (MQ=255)

                                                     ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                     >  1:2208704/1‑93 (MQ=255)

                                                     ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                     <  1:2168606/93‑1 (MQ=255)

                                                     ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAACTTCATCAGGCTCGCCAAACAGa                                     >  2:5094959/1‑93 (MQ=255)

                                                     tggCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                     <  1:7561067/91‑1 (MQ=255)

                                                      tGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAt                                    <  2:7665395/93‑1 (MQ=255)

                                                        cAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTg                                  >  1:5942269/1‑93 (MQ=255)

                                                         aGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgt                                 <  1:9981444/93‑1 (MQ=255)

                                                           cGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcg                               >  2:5325179/1‑93 (MQ=255)

                                                           cGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcg                               >  2:2337149/1‑93 (MQ=255)

                                                             gTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCCAACAGATTgtcgtc                             >  1:5434115/1‑93 (MQ=255)

                                                             gTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                             <  2:6362801/93‑1 (MQ=255)

                                                              tGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGtcgtcg                            <  1:112119/93‑1 (MQ=255)

                                                                 aCGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAATGTCGTCGTAg                         >  2:7317810/1‑93 (MQ=255)

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

                                                                   gTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTc                       <  2:6269713/93‑1 (MQ=255)

                                                                   ctCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTc                       <  2:4850706/92‑1 (MQ=255)

                                                                    tCGCCATCGGCGGATTCCATATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACCGATTGtcgtcgtcgtcg                      >  2:7881258/1‑93 (MQ=255)

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

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

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

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

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

                                                                       ccATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGc                   >  1:10041076/1‑93 (MQ=255)

                                                                       ccATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGGCGGGc                   >  1:1639636/1‑93 (MQ=255)

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

                                                                          tCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGTCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtc                <  2:2086729/93‑1 (MQ=255)

                                                                           cggcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtct               <  1:7446439/93‑1 (MQ=255)

                                                                            ggcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtctc              <  2:7178776/93‑1 (MQ=255)

                                                                             gcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTcc             <  1:4179274/93‑1 (MQ=255)

                                                                             gcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTcc             >  1:6879683/1‑93 (MQ=255)

                                                                              cggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATTCCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCt            <  1:1304144/93‑1 (MQ=255)

                                                                               ggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCACACAGATTGTCGTCGTAGCCGGGCTTCTCCtt           >  1:4476607/1‑93 (MQ=255)

                                                                               ggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTAGGGCTTCTCCtt           >  2:8010377/1‑93 (MQ=255)

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

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

                                                                                 aTTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAg         <  1:2531815/93‑1 (MQ=255)

                                                                                 aTTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAg         >  1:2782097/1‑93 (MQ=255)

                                                                                 aTTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAg         <  1:7266663/93‑1 (MQ=255)

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

                                                                                      acaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtt    >  1:2701211/1‑93 (MQ=255)

                                                                                      acaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtt    >  2:2543243/1‑93 (MQ=255)

GACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388841‑4389021

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: `‑`