BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,202,768	A→G	D97G (GAT→GGT)	ymfL →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,202,768	0	A	G	98.8%	520.9 / NA	162	D97G (GAT→GGT)	ymfL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (82/78); total (83/79)
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 165 total aligned reads displayed.

TTTTACAGCGTGCTGGTGGCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGATGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAATCGAAG  >  NC_000913/1202678‑1202857

    aCAGCGTGCTGGTGGCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGa                                                                                     <  2:373085/93‑1 (MQ=255)

    aCAGCGTGCTGGTGGCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGa                                                                                     >  1:335641/1‑93 (MQ=255)

       gCGTGCTGGTGGCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACaa                                                                                  >  2:3371867/1‑93 (MQ=255)

         gTGCTGGTGGCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTTTGGACAACg                                                                                >  1:1990703/1‑93 (MQ=255)

          tGCTGGTGGCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGc                                                                               <  1:4040468/93‑1 (MQ=255)

          tGCTGGTGGCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGc                                                                               >  1:4057493/1‑93 (MQ=255)

           gCTGGTGGCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGcc                                                                              <  1:586589/93‑1 (MQ=255)

             tggtggCACTCACTTGATTGCTGATGCTGTGGCGGAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGa                                                                            <  1:4080237/93‑1 (MQ=255)

             tggtggCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGa                                                                            <  1:3865633/93‑1 (MQ=255)

             tggtggCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGa                                                                            <  2:1079104/93‑1 (MQ=255)

             tggtggCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGATGTGGACAACGCCGa                                                                            <  2:3788963/93‑1 (MQ=255)

             tggtggCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTGCCTGACGTCGAGGGTGTGGACAACGCCGa                                                                            <  1:1737233/93‑1 (MQ=255)

               gtggCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGaga                                                                          >  1:295916/1‑91 (MQ=255)

                tggCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGatat                                                                         <  2:1424974/93‑1 (MQ=255)

                 ggCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATAtt                                                                        <  1:511953/93‑1 (MQ=255)

                   cactcactTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTaa                                                                      <  1:814994/93‑1 (MQ=255)

                   cactcactTCATTGCTGATGCTGGGGCCCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTaa                                                                      <  1:814993/93‑1 (MQ=255)

                    actcactTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAAc                                                                     <  2:1599681/93‑1 (MQ=255)

                    actcactTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAAc                                                                     <  1:1439186/93‑1 (MQ=255)

                    actcactTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAAc                                                                     <  2:2491130/93‑1 (MQ=255)

                     ctcactTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAAcc                                                                    >  1:1747517/1‑93 (MQ=255)

                      tcactTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCa                                                                   >  1:1355720/1‑93 (MQ=255)

                      tcactTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCa                                                                   <  1:4457797/93‑1 (MQ=255)

                      tcactTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCa                                                                   <  1:4457796/93‑1 (MQ=255)

                       cactTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAg                                                                  <  1:4687490/93‑1 (MQ=255)

                         ctTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCg                                                                <  2:1709574/93‑1 (MQ=255)

                          ttCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGt                                                               >  2:4266387/1‑93 (MQ=255)

                          ttCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGt                                                               >  1:3809070/1‑93 (MQ=255)

                           tCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTc                                                              <  1:2994612/93‑1 (MQ=255)

                             aTTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTctg                                                            <  2:2929988/93‑1 (MQ=255)

                              ttGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTctgc                                                           >  1:1226348/1‑93 (MQ=255)

                                gCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTctgctg                                                         <  2:4597949/93‑1 (MQ=255)

                                 cTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTgg                                                        >  2:3464430/1‑93 (MQ=255)

                                  tGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGa                                                       >  1:3694551/1‑93 (MQ=255)

                                  tGATGCTGTGGCGCAGTCGGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGa                                                       <  2:2792782/93‑1 (MQ=255)

                                   gATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGaa                                                      >  2:4688564/1‑93 (MQ=255)

                                     tGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGt                                                    >  1:2921380/1‑93 (MQ=255)

                                      gCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTc                                                   <  2:2981383/93‑1 (MQ=255)

                                       cTGTGGCGCCGTCTGCCAATGGCGTCTTTGTGCCTCTTCCTGACGTCGAGGGTGTGGCCAACGCCGCTATTAACCAGCGTCTGCTGGAAGTCa                                                  <  1:3614523/93‑1 (MQ=255)

                                        tgtgGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCAt                                                 >  2:1442763/1‑93 (MQ=255)

                                        tgtgGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCAt                                                 <  1:1837645/93‑1 (MQ=255)

                                        tgtgGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCAt                                                 >  2:438668/1‑93 (MQ=255)

                                        tgtgGCGCAGTCTGCAAATGGCGTCTTTGTGTCGCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCAt                                                 <  1:4576756/93‑1 (MQ=255)

                                          tgGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTg                                               <  1:2685231/93‑1 (MQ=255)

                                           ggCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGCGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGa                                              >  1:4232218/1‑93 (MQ=255)

                                            gcgcAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGaa                                             >  1:2836026/1‑93 (MQ=255)

                                              gcAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACa                                           <  1:2143226/93‑1 (MQ=255)

                                                 gTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGAt                                        >  1:3221992/1‑93 (MQ=255)

                                                 gTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGAt                                        >  2:3313112/1‑93 (MQ=255)

                                                 gTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGAt                                        >  1:2898779/1‑93 (MQ=255)

                                                 gTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGAt                                        >  2:1904735/1‑93 (MQ=255)

                                                    tGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCgg                                     >  1:4335402/1‑93 (MQ=255)

                                                     gCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGc                                    >  2:1635296/1‑93 (MQ=255)

                                                      cAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCa                                   <  1:3399084/93‑1 (MQ=255)

                                                       aaaTGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAg                                  >  1:2760382/1‑93 (MQ=255)

                                                       aaaTGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAg                                  >  2:1560023/1‑93 (MQ=255)

                                                         aTGGCTTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGtt                                <  1:624329/93‑1 (MQ=255)

                                                          tGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGttt                               >  1:3688859/1‑92 (MQ=255)

                                                          tGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTa                               <  2:2092009/93‑1 (MQ=255)

                                                           ggCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTAt                              >  1:2604134/1‑93 (MQ=255)

                                                           ggCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTAt                              >  2:3081878/1‑93 (MQ=255)

                                                            gCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTAtt                             <  2:1997573/93‑1 (MQ=255)

                                                            gCGTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTAtt                             <  2:1724359/93‑1 (MQ=255)

                                                              gTCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCa                           >  2:1162530/1‑93 (MQ=255)

                                                               tCTTTGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCaa                          <  1:1241973/93‑1 (MQ=255)

                                                                 tttGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTTaaaa                        >  2:247448/1‑93 (MQ=255)

                                                                 tttGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCaaaa                        >  2:4250525/1‑93 (MQ=255)

                                                                 tttGTGTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCaaaa                        >  2:4174928/1‑93 (MQ=255)

                                                                   tgtgTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACa                      >  1:670728/1‑93 (MQ=255)

                                                                   tgtgTCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACa                      >  1:3977447/1‑93 (MQ=255)

                                                                   tgtgTCTCTTCCTGACGTCGAGGAAGTGGACAACTCCGGTATTAACCAGCGCCTGCTGGAAGTCATTAAAAAGATC‑GCGGTAATTCaaaaaaa                     >  2:3479172/1‑90 (MQ=255)

                                                                   tgtgCCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACa                      <  1:2454693/93‑1 (MQ=255)

                                                                    gtgtCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAg                     <  1:1673150/93‑1 (MQ=255)

                                                                     tgtCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGa                    <  2:1297086/93‑1 (MQ=255)

                                                                     tgtCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGa                    <  2:4314992/93‑1 (MQ=255)

                                                                     tgtCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGa                    <  2:2571359/93‑1 (MQ=255)

                                                                      gtCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGAt                   >  2:2630559/1‑93 (MQ=255)

                                                                      gtCTCTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGAt                   <  1:3092750/93‑1 (MQ=255)

                                                                       tctcTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGAtt                  >  2:3232607/1‑93 (MQ=255)

                                                                       tctcTTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGAtt                  <  2:2329154/93‑1 (MQ=255)

                                                                        ctctTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTc                 <  1:2090593/93‑1 (MQ=255)

                                                                         tctTCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCg                >  1:1561638/1‑93 (MQ=255)

                                                                            tCCTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTc             <  2:1398330/93‑1 (MQ=255)

                                                                             tcTGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCa            <  1:1498577/92‑1 (MQ=255)

                                                                               tGACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGc          >  2:1884222/1‑93 (MQ=255)

                                                                                gACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCa         >  1:431920/1‑93 (MQ=255)

                                                                                gACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCa         >  1:1807109/1‑93 (MQ=255)

                                                                                gACGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCa         >  2:4511881/1‑93 (MQ=255)

                                                                                 aCGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCaa        >  2:3568435/1‑93 (MQ=255)

                                                                                  cGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAAt       <  2:1885494/93‑1 (MQ=255)

                                                                                  cGTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAAt       <  1:1060524/93‑1 (MQ=255)

                                                                                   gTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAGACAGATTCGTTGAGCAATc      >  1:2304623/1‑93 (MQ=255)

                                                                                   gTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAATc      >  1:1123724/1‑93 (MQ=255)

                                                                                   gTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAATc      <  2:2557756/93‑1 (MQ=255)

                                                                                   gTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAATc      >  1:4738878/1‑93 (MQ=255)

                                                                                   gTCGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAAAAGATCGGCAGTTATTCAAAACAGATTGGTTTCGCAAtt      >  2:1733849/1‑92 (MQ=255)

                                                                                     cGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAATCGa    >  2:521443/1‑93 (MQ=255)

                                                                                     cGAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAATCGa    >  2:900013/1‑93 (MQ=255)

                                                                                      gAGGGTGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAATCGaa   <  2:2956451/93‑1 (MQ=255)

                                                                                      gAGGATGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAATCGaa   >  2:3028480/1‑93 (MQ=255)

TTTTACAGCGTGCTGGTGGCACTCACTTCATTGCTGATGCTGTGGCGCAGTCTGCAAATGGCGTCTTTGTGTCTCTTCCTGACGTCGAGGATGTGGACAACGCCGATATTAACCAGCGTCTGCTGGAAGTCATTGAACAGATCGGCAGTTATTCAAAACAGATTCGTTCAGCAATCGAAG  >  NC_000913/1202678‑1202857

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