BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,440,519	A→G	E187E (GAA→GAG)	ytfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,440,519	0	A	G	99.0%	599.1 / ‑3.2	193	E187E (GAA→GAG)	ytfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (103/88); total (104/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 197 total aligned reads displayed.

TTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGA  >  NC_000913/4440429‑4440608

ttATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCg                                                                                         >  1:1819566/1‑93 (MQ=255)

 tATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGt                                                                                        <  2:3578802/93‑1 (MQ=255)

   tttCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTa                                                                                      <  1:2567702/93‑1 (MQ=255)

   tttCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTa                                                                                      <  2:3547285/93‑1 (MQ=255)

   tttCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTa                                                                                      >  1:1582836/1‑93 (MQ=255)

   tttCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTa                                                                                      >  2:3275132/1‑93 (MQ=255)

    ttCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTAc                                                                                     <  1:982867/93‑1 (MQ=255)

     tCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTAcc                                                                                    <  2:1585935/93‑1 (MQ=255)

         tAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGCTACCGCtt                                                                                <  2:4581800/93‑1 (MQ=255)

          aGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCttt                                                                               >  1:4152635/1‑93 (MQ=255)

           gCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTg                                                                              <  2:739577/93‑1 (MQ=255)

            cGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTgg                                                                             <  1:226819/93‑1 (MQ=255)

             gAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTggg                                                                            >  1:158918/1‑93 (MQ=255)

             gAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTggg                                                                            >  1:1629355/1‑93 (MQ=255)

             gAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTggg                                                                            >  1:3841727/1‑93 (MQ=255)

               aTTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCa                                                                          >  1:2790073/1‑93 (MQ=255)

                tttACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAt                                                                         <  1:3747313/93‑1 (MQ=255)

                 ttACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtg                                                                        <  2:4180143/93‑1 (MQ=255)

                  tACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtgt                                                                       <  1:4071885/93‑1 (MQ=255)

                  tACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtgt                                                                       <  2:3432304/93‑1 (MQ=255)

                  tACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtgt                                                                       <  1:1112110/93‑1 (MQ=255)

                  tACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtgt                                                                       <  1:2852364/93‑1 (MQ=255)

                   aCCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtgtg                                                                      >  1:835515/1‑93 (MQ=255)

                     cAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGAc                                                                    >  1:1973718/1‑93 (MQ=255)

                      aaaGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGAcc                                                                   >  2:1198482/1‑93 (MQ=255)

                      aaaGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGAcc                                                                   >  1:3006551/1‑93 (MQ=255)

                       aaGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCt                                                                  >  2:2700728/1‑93 (MQ=255)

                         gcgcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCttt                                                                >  1:45557/1‑93 (MQ=255)

                         gcgcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCttt                                                                >  1:2316231/1‑93 (MQ=255)

                         gcgcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCttt                                                                <  1:3488911/93‑1 (MQ=255)

                           gcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGa                                                              <  1:3080685/93‑1 (MQ=255)

                                 gggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGGGTTACCGCTTTGGGCATGAGACCTTTGAAGGAta                                                        >  2:401420/1‑92 (MQ=255)

                                   gCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcac                                                      >  1:1413198/1‑93 (MQ=255)

                                   gCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcac                                                      >  2:1367265/1‑93 (MQ=255)

                                   gCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcac                                                      <  2:2405649/93‑1 (MQ=255)

                                      ttGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACaaa                                                   <  2:3130051/93‑1 (MQ=255)

                                         cgcTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATcc                                                <  1:4522496/93‑1 (MQ=255)

                                          gcTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcg                                               >  2:2132288/1‑93 (MQ=255)

                                            tCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgcg                                             >  2:547819/1‑93 (MQ=255)

                                              ggCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGAt                                           <  2:2764195/93‑1 (MQ=255)

                                               gCCTGTATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATg                                          <  1:573582/93‑1 (MQ=255)

                                                  tGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAAt                                       <  2:2934542/93‑1 (MQ=255)

                                                  tGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAAt                                       <  1:2926353/93‑1 (MQ=255)

                                                  tGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAAt                                       >  1:2685548/1‑93 (MQ=255)

                                                   gCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATa                                      >  2:4642948/1‑93 (MQ=255)

                                                   gCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATa                                      >  2:4508582/1‑93 (MQ=255)

                                                    cATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAc                                     <  2:1665175/93‑1 (MQ=255)

                                                    cATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAc                                     <  1:2922135/93‑1 (MQ=255)

                                                      tAAAGCCTTCTGGGATATTGATTATAACGGTGGCGAGTGTTACTGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATGCCt                                   <  2:82888/93‑1 (MQ=255)

                                                      tAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCt                                   >  2:2599536/1‑93 (MQ=255)

                                                      tAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCt                                   >  2:2025646/1‑93 (MQ=255)

                                                       aaaGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCtc                                  >  1:124142/1‑92 (MQ=255)

                                                         aGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCa                                >  2:3231317/1‑93 (MQ=255)

                                                         aGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCa                                >  1:4018587/1‑93 (MQ=255)

                                                         aGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATAAATACCTGCa                                <  2:1255137/93‑1 (MQ=255)

                                                          gCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCaa                               >  2:4125421/1‑93 (MQ=255)

                                                           ccTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCaaa                              >  1:2165005/1‑93 (MQ=255)

                                                            cTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCaaaa                             >  2:1088647/1‑93 (MQ=255)

                                                            cTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCaaaa                             >  1:4682306/1‑93 (MQ=255)

                                                             ttCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAAt                            >  2:4356563/1‑93 (MQ=255)

                                                             ttCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAAt                            >  1:3447259/1‑93 (MQ=255)

                                                             ttCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAAt                            <  1:3322013/93‑1 (MQ=255)

                                                              tCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATc                           >  1:4186411/1‑93 (MQ=255)

                                                              tCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATc                           <  1:2472123/93‑1 (MQ=255)

                                                                tGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTg                         <  2:4312640/93‑1 (MQ=255)

                                                                 gggATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTgg                        >  1:3710579/1‑93 (MQ=255)

                                                                  ggATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTggg                       >  2:1928075/1‑92 (MQ=255)

                                                                   gATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTg                      >  1:1337694/1‑93 (MQ=255)

                                                                    atatTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGc                     <  1:1942082/93‑1 (MQ=255)

                                                                    atatTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGc                     <  2:1717239/93‑1 (MQ=255)

                                                                     tatTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGAGcc                    >  1:2005578/1‑93 (MQ=255)

                                                                      atTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCg                   >  1:3382671/1‑93 (MQ=255)

                                                                      atTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCg                   >  2:3699426/1‑93 (MQ=255)

                                                                       ttGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGt                  >  1:3195981/1‑93 (MQ=255)

                                                                       ttGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGt                  >  1:695410/1‑93 (MQ=255)

                                                                         gATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGttt                >  1:78612/1‑93 (MQ=255)

                                                                         gATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGttt                >  2:1334313/1‑93 (MQ=255)

                                                                          aTTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTa               >  1:1121233/1‑93 (MQ=255)

                                                                          aTTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTa               >  1:2231690/1‑93 (MQ=255)

                                                                            tataACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTaaa             <  1:1912616/93‑1 (MQ=255)

                                                                            tataACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTaaa             >  2:2846638/1‑93 (MQ=255)

                                                                             ataACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAag            <  2:2260123/93‑1 (MQ=255)

                                                                             ataACAGTGGCGAGCGTTACCGCTTTCGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAag            <  2:3792391/93‑1 (MQ=255)

                                                                              taACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAaga           >  2:2116427/1‑93 (MQ=255)

                                                                              taACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAaga           <  1:1353616/93‑1 (MQ=255)

                                                                              taACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAaga           >  2:1408900/1‑93 (MQ=255)

                                                                              taACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAaga           >  2:558231/1‑93 (MQ=255)

                                                                              taACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAaga           >  1:3948083/1‑93 (MQ=255)

                                                                              taACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAaga           >  1:3508314/1‑93 (MQ=255)

                                                                               aaCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAagag          <  2:4294440/93‑1 (MQ=255)

                                                                               aaCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAagag          <  2:636100/93‑1 (MQ=255)

                                                                               aaCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAagag          <  1:2735476/93‑1 (MQ=255)

                                                                               aaCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAagag          <  2:272513/93‑1 (MQ=255)

                                                                                aCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAgg         >  2:4455628/1‑93 (MQ=255)

                                                                                aCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAgg         >  1:2133027/1‑93 (MQ=255)

                                                                                aCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAgg         <  2:1639374/93‑1 (MQ=255)

                                                                                 cAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAggg        >  2:3738502/1‑93 (MQ=255)

                                                                                  aGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGc       <  2:179630/93‑1 (MQ=255)

                                                                                    tGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGa     <  1:4247690/93‑1 (MQ=255)

                                                                                       cGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGa  <  1:628226/93‑1 (MQ=255)

TTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGA  >  NC_000913/4440429‑4440608

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