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	98.8%	522.9 / ‑3.0	166	E187E (GAA→GAG)	ytfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (69/95); total (69/97)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 170 total aligned reads displayed.

TATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGA  >  NC_000913/4440430‑4440608

 agttCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGtt                                                                                       <  1:865193/91‑1 (MQ=255)

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

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

        tAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCtt                                                                                >  2:3931968/1‑93 (MQ=255)

        tAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCtt                                                                                >  1:2294466/1‑93 (MQ=255)

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

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

          gCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTg                                                                              >  2:3135274/1‑93 (MQ=255)

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

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

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

            gAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTggg                                                                            <  2:1316254/93‑1 (MQ=255)

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

             aaTGTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGc                                                                           <  1:3453934/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

                       aGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGGGGCGAGCGTTACCGCTTTGGGCATGTGACCtt                                                                 <  2:4061890/93‑1 (MQ=255)

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

                           cAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGaa                                                             <  1:3913876/93‑1 (MQ=255)

                           cAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGaa                                                             <  2:1172507/93‑1 (MQ=255)

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

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

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

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

                              cTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGa                                                          <  2:309524/93‑1 (MQ=255)

                              cTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGa                                                          <  1:4541646/93‑1 (MQ=255)

                              cTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGa                                                          <  2:4409466/93‑1 (MQ=255)

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

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

                               tGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGAt                                                         <  2:1500443/93‑1 (MQ=255)

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

                               tGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGAt                                                         <  1:4516960/93‑1 (MQ=255)

                                gggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATc                                                        >  2:420943/1‑93 (MQ=255)

                                gggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATc                                                        >  2:286002/1‑93 (MQ=255)

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

                                    aTTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACaa                                                    >  2:1103253/1‑93 (MQ=255)

                                     atGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACaaa                                                   <  1:1098448/92‑1 (MQ=255)

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

                                          cTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgc                                              <  1:1820899/93‑1 (MQ=255)

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

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

                                              gCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATg                                          >  2:2191889/1‑93 (MQ=255)

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

                                                 tGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGCTGAAt                                       >  2:439112/1‑93 (MQ=255)

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

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

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

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

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

                                                     tAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGTCATGTGACCTCTGAAGGATCACAAATCCGCGATGAATACCt                                   >  1:444034/1‑93 (MQ=255)

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

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

                                                     tAAAGCCTTCTGGGATATTGATTATAACAGTAGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCt                                   >  2:2437849/1‑93 (MQ=255)

                                                        aGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCa                                <  1:3405723/93‑1 (MQ=255)

                                                        aGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCa                                <  1:861715/93‑1 (MQ=255)

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

                                                          ccTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCaaa                              <  2:444722/93‑1 (MQ=255)

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

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

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

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

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

                                                                gggATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTgg                        <  1:4594817/93‑1 (MQ=255)

                                                                gggATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTgg                        <  2:742765/93‑1 (MQ=255)

                                                                 ggATATTGATTATAACAGTGGCGAGCGTTCCCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGt                       <  1:2658632/93‑1 (MQ=255)

                                                                  gATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTg                      <  2:2613035/93‑1 (MQ=255)

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

                                                                    tatTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGcc                    <  2:607247/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                           tataACAGTGGCGAGCGTTACCGCTTTGGGCATATGACCTTTGAAGGATTACAAATCCGCGATGGATACCTGCAAAATCTGGTGCCGTTTaaa             <  1:1173473/93‑1 (MQ=255)

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

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

                                                                              agcAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAagag          <  1:143837/91‑1 (MQ=255)

                                                                              aaCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAagag          >  1:1958680/1‑93 (MQ=255)

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

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

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

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

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

                                                                                 aGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGc       >  1:1372469/1‑93 (MQ=255)

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

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

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

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

                                                                                   tGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGa     >  2:591851/1‑93 (MQ=255)

                                                                                   gggCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGa     <  2:3580044/92‑1 (MQ=255)

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

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

                                                                                      cGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGa  >  1:204388/1‑93 (MQ=255)

TATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGA  >  NC_000913/4440430‑4440608

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