BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,758,261	T→C	E96G (GAA→GGA)	b1680 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,758,261	0	T	C	99.4%	510.3 / ‑4.8	165	E96G (GAA→GGA)	b1680	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (93/71); total (93/72)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 175 total aligned reads displayed.

TTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTTCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCATTTTCTCG  >  NC_000913/1758171‑1758350

ttAGCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGAtccctccc                                                                                         >  2:1875672/1‑93 (MQ=255)

  aGCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCcgt                                                                                       <  2:5079047/93‑1 (MQ=255)

   gCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCcgtc                                                                                      >  1:4093272/1‑93 (MQ=255)

   gCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCcgtc                                                                                      <  2:49397/93‑1 (MQ=255)

    cGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCcgtcg                                                                                     <  2:3691579/93‑1 (MQ=255)

      tggtgCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTg                                                                                   >  2:3489090/1‑93 (MQ=255)

      tggtgCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTg                                                                                   >  1:2283635/1‑93 (MQ=255)

       ggtgCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGc                                                                                  >  2:3118918/1‑93 (MQ=255)

       ggcgCTCCATCTGACTGATGATGATGTTATAGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGc                                                                                  >  2:974468/4‑93 (MQ=255)

        gtgCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGcc                                                                                 >  1:2073269/1‑93 (MQ=255)

         tgCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCcg                                                                                >  2:6166309/1‑93 (MQ=255)

          gCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCcgc                                                                               <  1:1687838/93‑1 (MQ=255)

          gCTCCATCTGACTGATGATGATGCTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCcgc                                                                               <  2:4328943/93‑1 (MQ=255)

            tCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCgg                                                                             >  2:7467254/1‑93 (MQ=255)

            tCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCgg                                                                             >  2:4459444/1‑93 (MQ=255)

            tCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCgg                                                                             >  2:7284251/1‑93 (MQ=255)

            tCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCgg                                                                             >  2:6055604/1‑93 (MQ=255)

              cATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGAc                                                                           >  2:12405/1‑93 (MQ=255)

              cATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGAc                                                                           >  1:854203/1‑93 (MQ=255)

               aTCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACg                                                                          >  2:6483578/1‑93 (MQ=255)

                tCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGa                                                                         >  2:2158301/1‑93 (MQ=255)

                 ctgactgaTGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGaa                                                                        <  1:2536780/93‑1 (MQ=255)

                  tgactgaTGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAac                                                                       >  1:5936852/1‑93 (MQ=255)

                  tgactgaTGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAac                                                                       >  1:506790/1‑93 (MQ=255)

                  tgactgaTGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAac                                                                       <  2:7770635/93‑1 (MQ=255)

                   gactgaTGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAaca                                                                      <  2:3128286/93‑1 (MQ=255)

                     ctgaTGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACTAGATTGATCCCTCCCGTCGTGCCGCGGACGAAcccc                                                                    <  2:5418724/93‑1 (MQ=255)

                     ctgaTGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACAcc                                                                    >  1:2804855/1‑93 (MQ=255)

                       gatgatgatGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCac                                                                  >  2:1732674/1‑92 (MQ=255)

                       gatgatgatGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAg                                                                  >  1:1844634/1‑93 (MQ=255)

                         tgatgatgTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGGCGAACACCACct                                                                >  1:1927365/1‑93 (MQ=255)

                         tgatgatgTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAgat                                                                >  2:4843055/1‑91 (MQ=255)

                          gatgatgTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGctc                                                               >  1:5432909/1‑93 (MQ=255)

                          gatgatgTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCCGACGAACACCCGctc                                                               >  1:5754823/1‑93 (MQ=255)

                             gatgTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTc                                                            <  2:6545756/93‑1 (MQ=255)

                             gatgTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTc                                                            <  2:2932408/93‑1 (MQ=255)

                             gatgTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTc                                                            <  1:6608288/93‑1 (MQ=255)

                             gatgTTACCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTc                                                            <  2:6128171/93‑1 (MQ=255)

                               tgTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTccg                                                          <  1:6390013/93‑1 (MQ=255)

                               tgTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTccg                                                          <  2:7666648/93‑1 (MQ=255)

                                gTTATCGCCGCCGCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTccgc                                                         <  2:6071245/93‑1 (MQ=255)

                                gTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTCTTGGCGACCAGCTTGATCCCTCCCGTCGTGCCGCGGACGAACCCCAGCTCTTccgc                                                         <  2:7038991/93‑1 (MQ=255)

                                gTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTccgc                                                         >  2:5051987/1‑93 (MQ=255)

                                 ttATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTccgcc                                                        <  1:126593/93‑1 (MQ=255)

                                  tATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTccgccg                                                       >  2:5178763/1‑93 (MQ=255)

                                  tATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTccgccg                                                       >  2:1072493/1‑93 (MQ=255)

                                   aTCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGa                                                      >  2:2248971/1‑93 (MQ=255)

                                   aTCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGa                                                      <  1:6500226/93‑1 (MQ=255)

                                   aTCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGa                                                      >  1:628141/1‑93 (MQ=255)

                                         cgcccgcACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGc                                                >  2:784945/1‑93 (MQ=255)

                                          gcccgcACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCa                                               >  1:4302965/1‑93 (MQ=255)

                                          gcccgcACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCa                                               <  1:5698134/93‑1 (MQ=255)

                                            ccgcACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCAtt                                             >  2:6730987/1‑93 (MQ=255)

                                              gcACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATtaa                                           >  2:267706/1‑93 (MQ=255)

                                               cACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATtaat                                          <  1:1352592/93‑1 (MQ=255)

                                               cACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATtaat                                          <  2:5454160/93‑1 (MQ=255)

                                                accttGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATtaata                                         >  1:5884408/4‑93 (MQ=255)

                                                aCGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATtaata                                         >  1:4528601/1‑93 (MQ=255)

                                                  gTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATaaa                                       <  2:5778033/93‑1 (MQ=255)

                                                   ttGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAAc                                      <  1:3172579/93‑1 (MQ=255)

                                                   ttGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAAc                                      >  2:6461636/1‑93 (MQ=255)

                                                   ttGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAAc                                      >  1:4712116/1‑93 (MQ=255)

                                                   ttGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAAc                                      <  2:5973971/93‑1 (MQ=255)

                                                    tGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACa                                     <  2:5558024/93‑1 (MQ=255)

                                                     gCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGCACACCAGCTCTTCCGCCGAACGGGCCTTAATAAagag                                    >  2:5742581/1‑93 (MQ=255)

                                                     gCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAg                                    <  2:1771655/93‑1 (MQ=255)

                                                      cTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGc                                   >  1:6145425/1‑93 (MQ=255)

                                                       tGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCNGCNNAACGGGCATTANTAAACANCg                                  >  1:562840/1‑93 (MQ=255)

                                                          tGCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATg                               <  1:38879/93‑1 (MQ=255)

                                                           gCCCCAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGc                              >  2:4232600/1‑93 (MQ=255)

                                                            ccccAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGcc                             >  2:4511519/1‑93 (MQ=255)

                                                             cccAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGccc                            >  1:442276/1‑93 (MQ=255)

                                                              ccAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCg                           >  1:5249520/1‑93 (MQ=255)

                                                              ccAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCg                           >  1:5673494/1‑93 (MQ=255)

                                                              ccAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCg                           >  1:2680371/1‑93 (MQ=255)

                                                              ccAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCg                           >  1:2870377/1‑93 (MQ=255)

                                                               cAGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGCACACCAGCTCTTCAGCCGAAAGGGCATTACTAAACAGCGATGccccc                          >  1:180709/1‑91 (MQ=255)

                                                                aGCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCt                         <  1:7722512/93‑1 (MQ=255)

                                                                 gCTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCtt                        >  2:5849996/1‑93 (MQ=255)

                                                                  cTATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTg                       <  1:3691392/93‑1 (MQ=255)

                                                                   tATTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTgc                      <  2:7018557/93‑1 (MQ=255)

                                                                    aTTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTgcg                     >  2:414715/1‑93 (MQ=255)

                                                                    aTTGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAAAACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGAGATGCCCGCGTgcg                     >  2:313860/1‑93 (MQ=255)

                                                                      tGGCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCa                   >  1:6965956/1‑93 (MQ=255)

                                                                       ggCGACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCAc                  >  2:5828006/1‑93 (MQ=255)

                                                                          gACCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGtt               >  1:2393456/1‑93 (MQ=255)

                                                                           aCCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTtc              >  1:3555951/1‑93 (MQ=255)

                                                                           aCCAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTtc              <  1:1580958/93‑1 (MQ=255)

                                                                             cggATTGACCCCACGCGTCTGGCCGCGGACGAACACCAGCTTTTCCGCCGAACGGGCATTAATAAACCGCGATGGCCGCTTGCGCACGTtctc            <  2:2274728/91‑1 (MQ=255)

                                                                             cAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTtctc            >  2:5225866/1‑93 (MQ=255)

                                                                             cAGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTtctc            <  1:5538658/93‑1 (MQ=255)

                                                                              aGATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTcc           >  1:4926197/1‑93 (MQ=255)

                                                                               gATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCa          <  2:1352805/93‑1 (MQ=255)

                                                                               gATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCa          <  1:2799342/93‑1 (MQ=255)

                                                                               gATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCa          >  1:1798113/1‑93 (MQ=255)

                                                                               gATTGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCa          <  1:1102381/93‑1 (MQ=255)

                                                                                  tGATCCCTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCAttt       <  1:5980046/93‑1 (MQ=255)

                                                                                     tccctcccGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCATTTtct    >  1:1038762/1‑93 (MQ=255)

                                                                                      cccTCCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCATTTtctc   <  2:2283634/93‑1 (MQ=255)

                                                                                       ccTTCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCATTTTCTCg  <  1:4486465/93‑1 (MQ=255)

TTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACGTTGCTGTTGCCCCAGCTATTGGCGACCAGATTGATCCCTTCCGTCGTGCCGCGGACGAACACCAGCTCTTCCGCCGAACGGGCATTAATAAACAGCGATGCCCGCTTGCGCACGTTCTCCATTTTCTCG  >  NC_000913/1758171‑1758350

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