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.7%	1064.1 / NA	336	E187E (GAA→GAG)	ytfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (170/165); minor base T (1/0); total (171/165)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

TTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGA  >  NC_000913/4440429‑4440608

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

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

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

    ttCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGGCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTAc                                                                                     >  1:2590390/1‑93 (MQ=255)

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

       gATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGc                                                                                  <  2:69974/93‑1 (MQ=255)

        aTAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCt                                                                                 >  2:2434454/1‑93 (MQ=255)

          aGCTAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCttt                                                                               <  1:2367740/93‑1 (MQ=255)

          aGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCttt                                                                               <  1:35704/93‑1 (MQ=255)

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

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

              aaTTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGc                                                                           >  1:1815097/1‑93 (MQ=255)

              aaTTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGc                                                                           >  2:7370136/1‑93 (MQ=255)

               aTTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCa                                                                          <  1:7555391/93‑1 (MQ=255)

                tttACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAt                                                                         >  1:7767275/1‑93 (MQ=255)

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

                  tACCAGAGCGCAGCTGGGCATTGCGCTCGTCCTGCATAAAGCCTTCTGCGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtgt                                                                       <  1:6576848/93‑1 (MQ=255)

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

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

                      aaaGAGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGAcc                                                                   >  2:4004668/1‑93 (MQ=255)

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

                        aGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCtt                                                                 <  1:2964226/93‑1 (MQ=255)

                        aGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCtt                                                                 <  2:3892409/93‑1 (MQ=255)

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

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

                           gcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGa                                                              >  1:631364/1‑93 (MQ=255)

                             aGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAg                                                            <  1:5960279/93‑1 (MQ=255)

                             aGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAg                                                            >  2:4450533/1‑93 (MQ=255)

                              gCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAgg                                                           >  1:2835830/1‑93 (MQ=255)

                               cTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGGCCTTTGAAGGa                                                          <  2:4568482/93‑1 (MQ=255)

                                tGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGAt                                                         >  1:1548461/1‑93 (MQ=255)

                                 gggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATc                                                        <  2:481730/93‑1 (MQ=255)

                                   gCATTGCGCTCGGCCTGCATAAAGCTTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCGTTGGGCATGTGACCGTTGAAGGATcac                                                      <  2:4988551/93‑1 (MQ=255)

                                   gCATTGCGCTCGGCCTGCATAAAGCCCTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcac                                                      >  1:1373785/1‑93 (MQ=255)

                                      ttGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACaaa                                                   >  1:3594320/1‑93 (MQ=255)

                                      ttGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACaaa                                                   >  2:1309716/1‑93 (MQ=255)

                                      ttGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACaaa                                                   >  2:1720448/1‑93 (MQ=255)

                                           cTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgc                                              >  2:6194496/1‑93 (MQ=255)

                                           cTCGGCCTGCATAAAGCCTTCNGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgc                                              <  1:416455/93‑1 (MQ=255)

                                            tCGGCCTGGATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgcg                                             >  2:1508287/1‑93 (MQ=255)

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

                                             cGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGa                                            <  2:1073408/93‑1 (MQ=255)

                                             cGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGa                                            <  2:5275957/93‑1 (MQ=255)

                                             cGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGa                                            <  1:7906198/93‑1 (MQ=255)

                                             cGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGa                                            <  1:5369071/93‑1 (MQ=255)

                                                ccTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGa                                         <  1:2504767/93‑1 (MQ=255)

                                                ccTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGa                                         >  2:7932224/1‑93 (MQ=255)

                                                 cTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGaa                                        >  2:4959968/1‑93 (MQ=255)

                                                 cTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGaa                                        <  1:3642710/93‑1 (MQ=255)

                                                 cTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGaa                                        <  2:2763637/93‑1 (MQ=255)

                                                  tGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGTGTTACCGCCTTGGGCATGTGACCTCTGAAGGATCACAAATCCGCGATGaaa                                       >  1:7001008/1‑92 (MQ=255)

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

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

                                                   gCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATa                                      <  1:1677475/93‑1 (MQ=255)

                                                    cATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAc                                     >  2:6259942/1‑93 (MQ=255)

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

                                                     aTAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATTAATAcc                                    <  2:7787562/93‑1 (MQ=255)

                                                     aTAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAcc                                    <  1:7731741/93‑1 (MQ=255)

                                                     aTAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAcc                                    >  1:6677842/1‑93 (MQ=255)

                                                     aTAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAcc                                    >  2:3008269/1‑93 (MQ=255)

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

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

                                                      tAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCt                                   <  1:7612756/93‑1 (MQ=255)

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

                                                       aaaGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTg                                  <  1:4343087/93‑1 (MQ=255)

                                                        aaGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGc                                 >  2:7640392/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                                                               cTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCt                          <  1:4066602/93‑1 (MQ=255)

                                                               cTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCt                          <  2:7491399/93‑1 (MQ=255)

                                                                tGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTg                         >  2:3117515/1‑93 (MQ=255)

                                                                tGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTg                         >  2:2043225/1‑93 (MQ=255)

                                                                  ggATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGt                       >  1:5448766/1‑93 (MQ=255)

                                                                  ggATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGt                       <  1:314917/93‑1 (MQ=255)

                                                                     tatTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGcc                    >  1:325186/1‑93 (MQ=255)

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

                                                                       ttGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGt                  <  1:7430595/93‑1 (MQ=255)

                                                                       ttGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGt                  <  1:4915571/93‑1 (MQ=255)

                                                                           ttATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTaa              >  1:5137504/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                 cAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAggg        <  2:876408/93‑1 (MQ=255)

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

                                                                                   gTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCg      >  1:886059/1‑93 (MQ=255)

                                                                                     ggCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGAt    <  1:572515/93‑1 (MQ=255)

                                                                                      gCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATg   <  2:5660385/93‑1 (MQ=255)

TTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGA  >  NC_000913/4440429‑4440608

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