BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	719,820	A→G	L5L (TTA→TTG)	ybfK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	719,820	0	A	G	86.5%	670.6 / 90.9	275	L5L (TTA→TTG)	ybfK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (21/16); new base G (118/119); total (139/136)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 276 total aligned reads displayed.

GAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGAC  >  NC_000913/719728‑719910

  aTGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGcc                                                                                          >  1:2058508/1‑93 (MQ=255)

   tGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCa                                                                                         >  1:6527382/1‑93 (MQ=255)

   tGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCa                                                                                         >  1:4564373/1‑93 (MQ=255)

   tGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCa                                                                                         >  2:4957126/1‑93 (MQ=255)

   tGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCa                                                                                         >  1:7214166/1‑93 (MQ=255)

      gaAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTc                                                                                      >  1:7037295/1‑93 (MQ=255)

       aaaTCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCa                                                                                     <  1:5740987/93‑1 (MQ=255)

       aaaTCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCa                                                                                     <  2:3941461/93‑1 (MQ=255)

         atcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAgg                                                                                   <  2:5771428/93‑1 (MQ=255)

         atcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAgg                                                                                   <  2:5272358/93‑1 (MQ=255)

          tcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGt                                                                                  <  2:7096550/93‑1 (MQ=255)

          tcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGt                                                                                  <  2:2785451/93‑1 (MQ=255)

          tcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGt                                                                                  <  2:3883476/93‑1 (MQ=255)

          tcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGt                                                                                  <  2:6368414/93‑1 (MQ=255)

           cgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTg                                                                                 <  1:6868859/93‑1 (MQ=255)

             atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGc                                                                               >  2:1997475/1‑93 (MQ=255)

             atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGc                                                                               >  1:2552981/1‑93 (MQ=255)

             atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGc                                                                               <  1:4682054/93‑1 (MQ=255)

             atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGc                                                                               >  1:4624470/1‑93 (MQ=255)

             atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGc                                                                               <  2:4582354/93‑1 (MQ=255)

             atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGc                                                                               <  1:1385168/93‑1 (MQ=255)

              tcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCt                                                                              >  1:3864339/1‑93 (MQ=255)

              tcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCt                                                                              >  1:4254319/1‑93 (MQ=255)

              tcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCt                                                                              <  1:6075447/93‑1 (MQ=255)

              tcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCt                                                                              >  1:2686878/1‑93 (MQ=255)

              tcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCt                                                                              >  1:7274090/1‑93 (MQ=255)

                gTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGa                                                                            <  1:5915679/93‑1 (MQ=255)

                 tGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGat                                                                           >  2:6264172/1‑93 (MQ=255)

                  gCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGata                                                                          <  1:3745429/93‑1 (MQ=255)

                   ccGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGataa                                                                         >  2:6587438/1‑93 (MQ=255)

                    cGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGataat                                                                        <  2:2930148/93‑1 (MQ=255)

                       aaaaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACa                                                                     <  1:3977323/93‑1 (MQ=255)

                       aaaaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACa                                                                     >  1:3479190/1‑93 (MQ=255)

                        aaaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACaa                                                                    <  1:6173692/93‑1 (MQ=255)

                        aaaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACaa                                                                    <  2:2243123/93‑1 (MQ=255)

                         aaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAg                                                                   >  1:7483923/1‑93 (MQ=255)

                          aaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAgg                                                                  >  2:6819544/1‑93 (MQ=255)

                              cAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTa                                                              >  2:1864138/1‑93 (MQ=255)

                               aaaGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTaa                                                             >  1:5137009/1‑93 (MQ=255)

                               aaaGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTaa                                                             <  1:2856263/93‑1 (MQ=255)

                                aaGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAAc                                                            <  2:6961682/93‑1 (MQ=255)

                                aaGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAAc                                                            >  2:4806279/1‑93 (MQ=255)

                                    tACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACagag                                                        <  1:5515002/93‑1 (MQ=255)

                                     aCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAgaga                                                       >  1:6351023/1‑93 (MQ=255)

                                     aCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAgaga                                                       <  1:3150350/93‑1 (MQ=255)

                                       cccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAAt                                                     <  1:1488085/93‑1 (MQ=255)

                                       cccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAAt                                                     <  2:2686042/93‑1 (MQ=255)

                                             gATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCAc                                               >  2:2329130/1‑93 (MQ=255)

                                               taaataaaTCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTg                                             >  1:895065/1‑93 (MQ=255)

                                                 aataaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                                           >  2:3719118/1‑93 (MQ=255)

                                                 aataaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                                           >  2:798700/1‑93 (MQ=255)

                                                 aataaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                                           <  1:2117719/93‑1 (MQ=255)

                                                   taaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGa                                         <  1:740930/93‑1 (MQ=255)

                                                    aaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                        <  2:4201984/93‑1 (MQ=255)

                                                    aaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                        <  2:2496428/93‑1 (MQ=255)

                                                    aaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                        <  2:146234/93‑1 (MQ=255)

                                                      atCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATa                                      <  2:6004032/93‑1 (MQ=255)

                                                        cAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACa                                    <  1:6701744/93‑1 (MQ=255)

                                                        cAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACa                                    <  1:6580695/93‑1 (MQ=255)

                                                        cAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACa                                    <  1:5229892/93‑1 (MQ=255)

                                                        cAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACa                                    <  1:1010791/93‑1 (MQ=255)

                                                          ggTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAgg                                  <  1:2746211/93‑1 (MQ=255)

                                                           gTAGATTTCTTTTTTCAAAATTGAGAAACCCCTTCCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAgg                                  <  1:1204890/93‑1 (MQ=255)

                                                           gTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGc                                 >  1:5980436/1‑93 (MQ=255)

                                                              gATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCata                              <  2:6644215/93‑1 (MQ=255)

                                                              gATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCata                              <  1:4445284/93‑1 (MQ=255)

                                                              gATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCata                              >  1:5689879/1‑93 (MQ=255)

                                                              gATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCata                              <  2:7245010/93‑1 (MQ=255)

                                                              gATTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCata                              <  1:7127914/93‑1 (MQ=255)

                                                               aTTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCatat                             >  1:1084891/1‑93 (MQ=255)

                                                               aTTGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCatat                             >  2:3416745/1‑93 (MQ=255)

                                                                ttGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATg                            >  2:5381923/1‑93 (MQ=255)

                                                                ttGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATg                            >  1:6794423/1‑93 (MQ=255)

                                                                 tGCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGc                           >  2:4048500/1‑93 (MQ=255)

                                                                  gCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGcc                          <  1:4115559/93‑1 (MQ=255)

                                                                  gCTTTTTTCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGcc                          >  1:3995149/1‑93 (MQ=255)

                                                                       tttCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTaa                     <  2:31126/93‑1 (MQ=255)

                                                                       tttCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTaa                     <  1:249729/93‑1 (MQ=255)

                                                                        ttCAGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAAt                    <  2:167882/93‑1 (MQ=255)

                                                                           aGAATGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACa                 >  2:3658057/1‑93 (MQ=255)

                                                                             aaTGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATTAGCACTGCGGAATACAGGCATATGCCTGTAATACaga               <  1:1848066/93‑1 (MQ=255)

                                                                              aTGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagag              >  2:46712/1‑93 (MQ=255)

                                                                               tGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagaga             <  2:4377679/93‑1 (MQ=255)

                                                                               tGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagaga             >  2:7365019/1‑93 (MQ=255)

                                                                               tGCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagaga             >  2:1288448/1‑93 (MQ=255)

                                                                                gCAGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagagag            <  1:5454916/93‑1 (MQ=255)

                                                                                gCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagagag            >  2:7042049/1‑93 (MQ=255)

                                                                                  aGAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTc          <  1:2916241/93‑1 (MQ=255)

                                                                                  aGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTc          <  1:4742928/93‑1 (MQ=255)

                                                                                   gAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCa         >  1:6261409/1‑93 (MQ=255)

                                                                                   gAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCa         <  1:818767/93‑1 (MQ=255)

                                                                                   gAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCa         >  1:125969/1‑93 (MQ=255)

                                                                                   gAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCa         >  2:5740847/1‑93 (MQ=255)

                                                                                   gAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCa         <  2:5575155/93‑1 (MQ=255)

                                                                                   gAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCa         <  1:5395576/93‑1 (MQ=255)

                                                                                   gAAA‑CCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCa         <  1:3230712/93‑1 (MQ=255)

                                                                                   gAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCa         >  1:1648289/1‑93 (MQ=255)

                                                                                         ccTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAg    >  1:5809675/1‑93 (MQ=255)

                                                                                          cTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGa   >  2:2546428/1‑93 (MQ=255)

                                                                                          cTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGa   >  1:3922504/1‑93 (MQ=255)

GAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAA‑CCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGAC  >  NC_000913/719728‑719910

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