BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,720,435	A→G	L821L (TTA→TTG)	yfiQ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,720,435	0	A	G	98.8%	514.9 / ‑4.6	163	L821L (TTA→TTG)	yfiQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (84/77); total (86/77)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 168 total aligned reads displayed.

GAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAA  >  NC_000913/2720346‑2720524

  gaTCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCtt                                                                                      >  2:1034948/1‑93 (MQ=255)

  gaTCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCtt                                                                                      >  2:3811659/1‑93 (MQ=255)

   aTCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTa                                                                                     >  2:1338357/1‑93 (MQ=255)

    tcctccGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAg                                                                                    <  2:1228407/93‑1 (MQ=255)

    tCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAg                                                                                    <  1:3194383/93‑1 (MQ=255)

     ccTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAgg                                                                                   <  1:4076680/93‑1 (MQ=255)

     ccTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAgg                                                                                   >  2:2886798/1‑93 (MQ=255)

      cTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCACGGATCTCAAAGGGGTTGGCTGAgga                                                                                  >  2:2650862/1‑92 (MQ=255)

        cGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGGTGGGCTTAGGTCg                                                                                >  2:174917/1‑93 (MQ=255)

         ggCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGa                                                                               >  1:2191791/1‑93 (MQ=255)

         ggCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGa                                                                               <  2:646706/93‑1 (MQ=255)

          gCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGAc                                                                              <  1:9112/93‑1 (MQ=255)

          gCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGAc                                                                              >  1:527787/1‑93 (MQ=255)

          gCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGAc                                                                              >  1:3446716/1‑93 (MQ=255)

           cgacACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTAGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGAct                                                                             <  2:2761265/90‑2 (MQ=255)

            gTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGc                                                                            >  2:2218192/1‑93 (MQ=255)

            gTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGc                                                                            <  1:2402805/93‑1 (MQ=255)

            gTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGc                                                                            <  2:1787614/93‑1 (MQ=255)

              cACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCtt                                                                          <  2:3401891/93‑1 (MQ=255)

              cACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCtt                                                                          <  1:1735823/93‑1 (MQ=255)

               aCGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTa                                                                         >  1:3939471/1‑93 (MQ=255)

                 gcgTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAAt                                                                       <  1:3835725/93‑1 (MQ=255)

                  cgTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATg                                                                      >  1:449362/1‑93 (MQ=255)

                   gTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTAGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATgg                                                                     <  2:3021650/93‑1 (MQ=255)

                    tGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGa                                                                    >  2:4439105/1‑93 (MQ=255)

                     gCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGaa                                                                   <  2:664165/93‑1 (MQ=255)

                        aTTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTCGGTCGACGCTTAATGGAAAAg                                                                >  1:1676341/1‑93 (MQ=255)

                          ttCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGtt                                                              >  1:3254315/1‑93 (MQ=255)

                          ttCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGtt                                                              >  2:3982479/1‑93 (MQ=255)

                           tCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTg                                                             >  2:3451214/1‑93 (MQ=255)

                           tCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTg                                                             >  1:1805242/1‑93 (MQ=255)

                             cGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAt                                                           <  2:2137973/93‑1 (MQ=255)

                             cGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAt                                                           <  1:2159701/93‑1 (MQ=255)

                              gATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAtt                                                          >  1:1271872/1‑93 (MQ=255)

                              gATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAtt                                                          <  2:2936620/93‑1 (MQ=255)

                              gATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAtt                                                          >  2:3673438/1‑93 (MQ=255)

                               aTCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTa                                                         >  1:375590/1‑93 (MQ=255)

                                 ccTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTAcc                                                       <  1:2220422/93‑1 (MQ=255)

                                  cTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCt                                                      >  1:2450493/1‑93 (MQ=255)

                                  cTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCt                                                      >  1:1501652/1‑93 (MQ=255)

                                   tGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCta                                                     <  1:1893415/93‑1 (MQ=255)

                                   tGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCta                                                     >  2:683469/1‑93 (MQ=255)

                                   tGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCta                                                     >  2:2986713/1‑93 (MQ=255)

                                    gATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCtat                                                    >  2:676979/1‑93 (MQ=255)

                                     aTAACATCGATGGCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTGCCtata                                                   <  1:4034181/93‑1 (MQ=255)

                                     aTAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCtata                                                   <  1:4671052/93‑1 (MQ=255)

                                     aTAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCtata                                                   >  1:2813332/1‑93 (MQ=255)

                                     aTAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCtata                                                   >  1:381500/1‑93 (MQ=255)

                                     aTAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCtata                                                   <  2:41298/93‑1 (MQ=255)

                                      tgaCAACGAGGCGGAATTGTCGGTGCTGGTTCGCTCGGATCGCAAAGGGTTGGGCTTAGGTCGACAATTAATGGAAAAGTTGATTACCTATAc                                                  <  1:2697092/91‑1 (MQ=255)

                                      tAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATAc                                                  <  1:2697093/93‑1 (MQ=255)

                                         cATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATAcgcg                                               >  2:985819/1‑93 (MQ=255)

                                         cATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATAcgcg                                               >  2:482840/1‑93 (MQ=255)

                                         cATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATAcgcg                                               >  2:482842/1‑93 (MQ=255)

                                           tCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCgag                                             >  2:1816482/1‑93 (MQ=255)

                                           tCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGGTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTCCCTATACGCgag                                             >  1:160479/1‑93 (MQ=255)

                                            cGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCgaga                                            <  1:793568/93‑1 (MQ=255)

                                            cGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCgaga                                            >  2:424656/1‑93 (MQ=255)

                                            cGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCgaga                                            >  1:2322325/1‑93 (MQ=255)

                                             gATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGAt                                           <  2:2381126/93‑1 (MQ=255)

                                              aTGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATc                                          <  2:338286/93‑1 (MQ=255)

                                               tGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCa                                         <  2:1444616/93‑1 (MQ=255)

                                               tGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCa                                         <  2:4079816/93‑1 (MQ=255)

                                                gCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCAc                                        >  2:3783021/1‑93 (MQ=255)

                                                 ccGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACg                                       <  1:4245519/93‑1 (MQ=255)

                                                  cGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACgg                                      <  1:1967699/93‑1 (MQ=255)

                                                   gAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGa                                     >  2:1060032/1‑93 (MQ=255)

                                                   gAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGa                                     <  1:2928310/93‑1 (MQ=255)

                                                   gAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGa                                     >  1:3692028/1‑93 (MQ=255)

                                                     aTTTGCTGTACTGGTTCGGCCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACt                                   <  1:683938/93‑1 (MQ=255)

                                                     aTTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACt                                   >  1:146442/1‑93 (MQ=255)

                                                      tttGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTa                                  >  1:2775197/1‑93 (MQ=255)

                                                      tttGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTa                                  <  1:367073/93‑1 (MQ=255)

                                                         gCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACaa                               <  1:738522/93‑1 (MQ=255)

                                                         gCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACaa                               <  1:1146487/93‑1 (MQ=255)

                                                           tGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACg                             >  2:2566419/1‑93 (MQ=255)

                                                            gTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGt                            <  2:1914483/93‑1 (MQ=255)

                                                            gTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGt                            >  2:4531267/1‑93 (MQ=255)

                                                              aCTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGTGAGATCACGGACTACAACGTCt                          <  2:2992242/93‑1 (MQ=255)

                                                               cTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTg                         >  2:3414238/1‑93 (MQ=255)

                                                                tGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGa                        >  2:2355736/1‑93 (MQ=255)

                                                                 ggTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGaa                       >  2:3460302/1‑93 (MQ=255)

                                                                     cGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGt                   <  2:3155597/93‑1 (MQ=255)

                                                                     cGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGt                   <  2:4562449/93‑1 (MQ=255)

                                                                     cGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGt                   <  1:83997/93‑1 (MQ=255)

                                                                      gCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTa                  >  2:2894866/1‑93 (MQ=255)

                                                                      gCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTa                  >  1:438666/1‑93 (MQ=255)

                                                                      gCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCCTAATGGCAAAGCTGATTACCTATACGCGAGATCACGGACTACAACGCCTGACTggga                  >  2:2894865/1‑91 (MQ=255)

                                                                        gcGGAACTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTGCCTATACGCGAGATCACGGACTACAACGTCTGAATGGTAtt                <  2:2085558/92‑1 (MQ=255)

                                                                          ggATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAc              <  2:1532931/93‑1 (MQ=255)

                                                                            aTCTCAAAGAGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATGACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGa            <  1:2146165/93‑1 (MQ=255)

                                                                             tctcAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTCCAACGTCTGAATGGTATTACGag           >  1:1090200/1‑92 (MQ=255)

                                                                             tctcAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGAt           >  2:1070897/1‑93 (MQ=255)

                                                                              ctcAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATg          >  1:215559/1‑93 (MQ=255)

                                                                               tcAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGc         >  1:3088144/1‑93 (MQ=255)

                                                                                   aGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCaaa     >  2:1238516/1‑93 (MQ=255)

                                                                                    gggTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAAc    >  1:2891302/1‑93 (MQ=255)

                                                                                    gggTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAAc    <  1:262280/93‑1 (MQ=255)

                                                                                      gTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACaa  <  1:2319309/93‑1 (MQ=255)

                                                                                      gTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACaa  <  2:1192023/93‑1 (MQ=255)

GAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAA  >  NC_000913/2720346‑2720524

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