BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,573,190	Δ1 bp	coding (1487/1542 nt)	yjiU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,189	0	G	.	99.5%	973.3 / ‑3.2	207	coding (1486/1542 nt)	yjiU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base . (107/99); minor base T (0/1); total (107/100)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

TTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAACCC  >  NC_000913/4573097‑4573277

ttatCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAg                                                                                          >  2:2531297/1‑93 (MQ=255)

ttatCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAg                                                                                          >  1:4269949/1‑93 (MQ=255)

  atCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAgag                                                                                        >  1:3215196/1‑91 (MQ=255)

  atCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAgag                                                                                        >  2:5484077/1‑91 (MQ=255)

   tCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAgagc                                                                                       <  2:2572022/93‑4 (MQ=255)

    cAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAgagca                                                                                      <  1:3803296/93‑5 (MQ=255)

     aaTACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCat                                                                                    <  2:1272008/93‑1 (MQ=255)

       tACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCatat                                                                                  >  2:3539335/1‑93 (MQ=255)

        aCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATg                                                                                 <  1:4923026/93‑1 (MQ=255)

          aTGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGAc                                                                               <  2:5277355/93‑1 (MQ=255)

           tGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAGCCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACt                                                                              <  2:3302712/93‑1 (MQ=255)

           tGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACt                                                                              >  2:3904098/1‑93 (MQ=255)

            ggCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACtt                                                                             >  2:4329184/1‑93 (MQ=255)

             gCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACttt                                                                            >  2:107931/1‑93 (MQ=255)

               tGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACttttt                                                                          >  2:5616967/1‑93 (MQ=255)

               tGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACttttt                                                                          >  2:4045640/1‑93 (MQ=255)

               tGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACttttt                                                                          >  1:4694876/1‑93 (MQ=255)

               tGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACttttt                                                                          >  1:4840858/1‑93 (MQ=255)

                gAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTc                                                                         >  1:1205463/1‑93 (MQ=255)

                       gCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATc                                                                  >  1:2665076/1‑93 (MQ=255)

                       gCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATc                                                                  <  1:1817413/93‑1 (MQ=255)

                        cATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCa                                                                 >  2:400960/1‑93 (MQ=255)

                         aTGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCAc                                                                >  2:2243656/1‑93 (MQ=255)

                          tGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACg                                                               <  1:1334494/93‑1 (MQ=255)

                          tGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACg                                                               <  2:3439127/93‑1 (MQ=255)

                            ggCCCAACTTACCGGACCGCTTGCAGTACTGAAGCCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGa                                                             <  1:1840569/93‑1 (MQ=255)

                            ggCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGa                                                             <  2:3667693/93‑1 (MQ=255)

                            ggCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGa                                                             >  2:1378083/1‑93 (MQ=255)

                             gCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGAt                                                            <  1:1018293/93‑1 (MQ=255)

                             gCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGAt                                                            <  2:4174559/93‑1 (MQ=255)

                               ccAACTTACCGGACCGCTTGCAGTACTGAAGCCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGc                                                          <  1:953951/93‑1 (MQ=255)

                                cAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGcc                                                         <  2:5729867/93‑1 (MQ=255)

                                cAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGcc                                                         >  1:5881443/1‑93 (MQ=255)

                                 aaCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCt                                                        <  2:1293897/93‑1 (MQ=255)

                                  aCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTg                                                       >  2:6194880/1‑93 (MQ=255)

                                  aCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTg                                                       >  1:5756037/1‑93 (MQ=255)

                                   cTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTgg                                                      <  1:1915364/93‑1 (MQ=255)

                                    ttACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGa                                                     >  1:6426505/1‑93 (MQ=255)

                                    ttACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGa                                                     >  1:188807/1‑93 (MQ=255)

                                    ttACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGa                                                     <  1:561685/93‑1 (MQ=255)

                                    ttACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGa                                                     >  1:528687/1‑93 (MQ=255)

                                      aCCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAtt                                                   <  1:5843721/93‑1 (MQ=255)

                                       ccGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAttt                                                  <  2:2634573/93‑1 (MQ=255)

                                        cGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTg                                                 >  2:2539917/1‑93 (MQ=255)

                                          gACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGcc                                               <  1:5719908/93‑1 (MQ=255)

                                          gACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGcc                                               <  2:6092176/93‑1 (MQ=255)

                                           gccGCTTGCAGTACTGAAACCCGATGTCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGccc                                              <  2:312061/92‑1 (MQ=255)

                                           aCCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGccc                                              >  1:623987/1‑93 (MQ=255)

                                           aCCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGccc                                              >  2:3771424/1‑93 (MQ=255)

                                           aCCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGccc                                              >  2:2935726/1‑93 (MQ=255)

                                           aCCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGccc                                              >  1:39778/1‑93 (MQ=255)

                                            ccGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCAAATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCg                                             >  1:2608390/1‑93 (MQ=255)

                                              gCTTGCAGCACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGt                                           >  2:1661382/1‑93 (MQ=255)

                                                ttGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTaa                                         >  2:5020988/1‑93 (MQ=255)

                                                ttGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTaa                                         >  2:1519726/1‑93 (MQ=255)

                                                ttGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTaa                                         >  1:3075588/1‑93 (MQ=255)

                                                  gCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAcc                                       <  1:2887295/93‑1 (MQ=255)

                                                    aGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCaa                                     <  2:5192589/93‑1 (MQ=255)

                                                    aGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCaa                                     <  1:2212927/93‑1 (MQ=255)

                                                       aCTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaa                                  <  1:4091740/93‑1 (MQ=255)

                                                        cTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaaa                                 >  1:3168554/1‑93 (MQ=255)

                                                         tGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaaaa                                >  2:1731439/1‑93 (MQ=255)

                                                              cccGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTg                           <  2:3285280/93‑1 (MQ=255)

                                                              cccGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTg                           >  1:2816226/1‑93 (MQ=255)

                                                               ccGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTgg                          >  2:3825004/1‑93 (MQ=255)

                                                               ccGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTgg                          >  1:3491067/1‑93 (MQ=255)

                                                               ccGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGCAACCAATAAAATCCTga                          >  1:3303614/1‑92 (MQ=255)

                                                                cGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGa                         <  2:6236319/93‑1 (MQ=255)

                                                                 gATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGAt                        <  2:3254839/93‑1 (MQ=255)

                                                                 gATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGAt                        >  2:4724979/1‑93 (MQ=255)

                                                                  aTGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATa                       >  2:4399165/1‑93 (MQ=255)

                                                                  aTGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATa                       >  2:4843522/1‑93 (MQ=255)

                                                                  aTGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATa                       <  1:6640379/93‑1 (MQ=255)

                                                                  aTGGCAACCGTTTTTATTTACCG‑AATGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATa                       <  2:6421817/93‑1 (MQ=255)

                                                                   tGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATAc                      <  2:4618064/93‑1 (MQ=255)

                                                                    ggCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACg                     <  2:278764/93‑1 (MQ=255)

                                                                       aaCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCt                  >  1:5796069/1‑93 (MQ=255)

                                                                        aCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCtt                 >  2:773593/1‑93 (MQ=255)

                                                                        aCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCtt                 >  2:4862284/1‑93 (MQ=255)

                                                                        aCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCtt                 <  2:4555780/93‑1 (MQ=255)

                                                                         ccGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCttt                <  1:5455984/93‑1 (MQ=255)

                                                                          cGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCtttt               >  2:5701584/1‑93 (MQ=255)

                                                                          cGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCtttt               >  1:1107261/1‑93 (MQ=255)

                                                                           gTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTa              >  1:502315/1‑93 (MQ=255)

                                                                             ttttATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTaaa            <  2:3514206/93‑1 (MQ=255)

                                                                                tatttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGc         >  2:1597602/1‑93 (MQ=255)

                                                                                tatttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGc         <  1:434009/93‑1 (MQ=255)

                                                                                tatttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGc         <  1:4123911/93‑1 (MQ=255)

                                                                                 atttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCt        <  2:1846577/93‑1 (MQ=255)

                                                                                  tttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCtt       >  2:4801438/1‑93 (MQ=255)

                                                                                  tttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCtt       >  1:2310447/1‑93 (MQ=255)

                                                                                  tgtaCCGAAA‑GAGCATATGACTTTTTCTGTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCtt       <  2:2502577/91‑1 (MQ=255)

                                                                                   ttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTa      >  2:6546598/1‑93 (MQ=255)

                                                                                    taCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTaa     <  1:2900570/93‑1 (MQ=255)

                                                                                     aCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAAc    <  2:4114079/93‑1 (MQ=255)

                                                                                     aCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAAc    <  2:936229/93‑1 (MQ=255)

                                                                                      ccGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAAcc   <  2:2145545/93‑1 (MQ=255)

                                                                                      ccGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAAcc   <  1:473402/93‑1 (MQ=255)

                                                                                       cGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAAccc  >  1:854660/1‑93 (MQ=255)

                                                                                       cGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAAccc  >  1:3449861/1‑93 (MQ=255)

TTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAACCC  >  NC_000913/4573097‑4573277

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