BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	981,566	C→T	R433C (CGC→TGC)	ycbB →	putative amidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	981,566	0	C	T	89.1%	838.3 / 57.7	275	R433C (CGC→TGC)	ycbB	putative amidase
Reads supporting (aligned to +/- strand): ref base C (13/17); new base T (117/128); total (130/145)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 275 total aligned reads displayed.

GGTCGCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCAGAGAAGC  >  NC_000913/981476‑981657

ggtcgCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGc                                                                                           <  1:1125702/93‑1 (MQ=255)

ggtcgCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGc                                                                                           <  2:594851/93‑1 (MQ=255)

  tcgCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCaa                                                                                         <  1:3534115/93‑1 (MQ=255)

    gCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCACCTCTGGCATGCAAAg                                                                                       >  1:3411492/1‑93 (MQ=255)

     ccccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGa                                                                                      <  2:3359850/93‑1 (MQ=255)

     ccccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGa                                                                                      >  2:2196924/1‑93 (MQ=255)

     ccccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGa                                                                                      <  2:4407529/93‑1 (MQ=255)

      cccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGat                                                                                     >  2:258040/1‑93 (MQ=255)

       ccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGata                                                                                    >  2:5346792/1‑93 (MQ=255)

       ccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGata                                                                                    >  2:3486051/1‑93 (MQ=255)

       ccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGata                                                                                    >  1:1987080/1‑93 (MQ=255)

       ccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGata                                                                                    >  1:5838791/1‑93 (MQ=255)

        cGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGatat                                                                                   >  1:4941467/1‑93 (MQ=255)

        cGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGatat                                                                                   >  2:2344628/1‑93 (MQ=255)

         gATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATAtt                                                                                  <  1:3772569/93‑1 (MQ=255)

         gATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATAtt                                                                                  <  2:2429939/93‑1 (MQ=255)

         gATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGATAtt                                                                                  >  2:7371370/1‑93 (MQ=255)

          aTCGCAAACCGCCGATGATGAGCAGTCCCCTTAACAAAGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTc                                                                                 <  1:1240166/93‑1 (MQ=255)

          aTCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTc                                                                                 >  2:1867533/1‑93 (MQ=255)

          aTCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTc                                                                                 >  1:4669518/1‑93 (MQ=255)

           tCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCt                                                                                >  2:1365203/1‑93 (MQ=255)

           tCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCt                                                                                <  1:716171/93‑1 (MQ=255)

           tCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCt                                                                                <  2:91514/93‑1 (MQ=255)

            cGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTg                                                                               <  2:5592657/93‑1 (MQ=255)

             gCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGc                                                                              >  2:7117879/1‑93 (MQ=255)

             gCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGc                                                                              >  2:6828871/1‑93 (MQ=255)

              cAAAACGCCGATGAGGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGcc                                                                             <  1:907082/93‑1 (MQ=255)

               aaaaCGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGcnn                                                                            >  1:1789122/1‑91 (MQ=255)

               aaaaCGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCa                                                                            <  1:2932672/93‑1 (MQ=255)

                  aCGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCaaaa                                                                         <  2:2949909/93‑1 (MQ=255)

                  aCGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCaaaa                                                                         <  2:4153553/93‑1 (MQ=255)

                   cGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAg                                                                        <  2:2051738/93‑1 (MQ=255)

                   cGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAATAg                                                                        <  2:4735437/93‑1 (MQ=255)

                       gatgatGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTgcg                                                                    <  1:1187421/93‑1 (MQ=255)

                       gatgatGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTgcg                                                                    >  2:214198/1‑93 (MQ=255)

                         tgatgaGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCa                                                                  >  1:5933477/1‑93 (MQ=255)

                         tgatgaGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCa                                                                  >  1:4118321/1‑93 (MQ=255)

                           atgaGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAAc                                                                <  2:1429871/93‑1 (MQ=255)

                           atgaGCAGTGCCCTAAACAGCGTAGTGTTAATCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAAc                                                                <  1:1527991/93‑1 (MQ=255)

                            tgaGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACg                                                               >  2:4997464/1‑93 (MQ=255)

                             gaGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGa                                                              <  1:7227452/93‑1 (MQ=255)

                             gaGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGa                                                              >  2:5460365/1‑93 (MQ=255)

                              aGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGAt                                                             >  2:6643986/1‑93 (MQ=255)

                               gCAGTGCCCTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATc                                                            <  1:3355777/93‑1 (MQ=255)

                                cAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATcc                                                           <  2:2212093/93‑1 (MQ=255)

                                  gTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCgg                                                         <  1:4726723/93‑1 (MQ=255)

                                    gCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGa                                                       <  1:2104302/93‑1 (MQ=255)

                                     cccTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGat                                                      >  2:1058402/1‑93 (MQ=255)

                                      ccTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGata                                                     <  1:4020832/93‑1 (MQ=255)

                                      ccTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGata                                                     >  1:3488804/1‑93 (MQ=255)

                                      ccTAAACAACTTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGata                                                     <  1:6023045/93‑1 (MQ=255)

                                       cTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGatat                                                    >  1:2452271/1‑93 (MQ=255)

                                        ttAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATAtc                                                   >  2:3625542/1‑93 (MQ=255)

                                        ttAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATAtc                                                   >  2:3103395/1‑93 (MQ=255)

                                          aacaacGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATAtctc                                                 <  1:3604447/93‑1 (MQ=255)

                                           acaacGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCg                                                >  1:3621111/1‑93 (MQ=255)

                                               cGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAg                                            >  2:6515579/1‑93 (MQ=255)

                                                gTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGc                                           >  2:4925487/1‑93 (MQ=255)

                                                  aGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCa                                         <  1:3117686/93‑1 (MQ=255)

                                                   gTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGACAt                                        >  2:4285284/1‑93 (MQ=255)

                                                    tGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATg                                       >  1:7329450/1‑93 (MQ=255)

                                                      gTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGc                                     <  1:6348821/93‑1 (MQ=255)

                                                        aaaCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCta                                   >  1:2141576/1‑93 (MQ=255)

                                                         aaCCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCtat                                  <  1:1615618/93‑1 (MQ=255)

                                                          aCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCtata                                 >  2:3916307/1‑93 (MQ=255)

                                                          aCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCtata                                 <  1:6631735/93‑1 (MQ=255)

                                                           cccGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATAc                                >  2:1698779/1‑93 (MQ=255)

                                                            ccgccgTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACg                               <  1:6519472/93‑1 (MQ=255)

                                                            ccgccgTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACg                               >  1:3262769/1‑93 (MQ=255)

                                                              gccgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGt                             <  1:5304407/93‑1 (MQ=255)

                                                              gccgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGt                             <  1:7000925/93‑1 (MQ=255)

                                                              gccgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGt                             <  2:326048/93‑1 (MQ=255)

                                                                cgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGa                           <  2:6076513/93‑1 (MQ=255)

                                                                cgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGa                           <  2:6286136/93‑1 (MQ=255)

                                                                cgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGa                           <  2:1025841/93‑1 (MQ=255)

                                                                 gTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGAt                          <  1:1013521/93‑1 (MQ=255)

                                                                 gTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGAt                          >  2:706267/1‑93 (MQ=255)

                                                                  tGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATg                         >  1:5436171/1‑93 (MQ=255)

                                                                  tGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATg                         >  2:5243901/1‑93 (MQ=255)

                                                                     aaCGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATgcgc                      >  2:5766702/1‑93 (MQ=255)

                                                                     aaCGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATgcgc                      <  1:5466149/93‑1 (MQ=255)

                                                                     aaCGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATgcgc                      <  2:3213019/93‑1 (MQ=255)

                                                                      aCGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGcgcg                     >  2:3715474/1‑93 (MQ=255)

                                                                      aCGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGcgcg                     <  1:851763/93‑1 (MQ=255)

                                                                       cGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCgg                    <  2:1566712/93‑1 (MQ=255)

                                                                        gTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGc                   >  2:6588290/1‑93 (MQ=255)

                                                                        gTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGc                   >  1:708756/1‑93 (MQ=255)

                                                                         tACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCt                  >  1:4451588/1‑93 (MQ=255)

                                                                           ccTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTgg                >  2:4673498/1‑93 (MQ=255)

                                                                             tCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGaa              <  2:2570141/93‑1 (MQ=255)

                                                                             tCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGaa              >  1:389940/1‑93 (MQ=255)

                                                                              ccAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAAc             <  1:6170183/93‑1 (MQ=255)

                                                                              ccAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAAc             >  2:1215493/1‑93 (MQ=255)

                                                                                aaCTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAAcag           >  1:1392515/1‑93 (MQ=255)

                                                                                   tctGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAAcagcag        >  1:2800526/1‑93 (MQ=255)

                                                                                    ctGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCaga       <  2:5181720/93‑1 (MQ=255)

                                                                                     tGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCagag      >  1:554810/1‑93 (MQ=255)

                                                                                      ggCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCAgaga     <  2:6766673/93‑1 (MQ=255)

                                                                                      ggCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCAgaga     <  1:692872/93‑1 (MQ=255)

                                                                                       gCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCAGAGaa    <  1:4228266/93‑1 (MQ=255)

GGTCGCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCAGAGAAGC  >  NC_000913/981476‑981657

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