BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,493,567	T→C	F86L (TTT→CTT)	b1423 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,493,567	0	T	C	97.8%	422.8 / ‑0.8	136	F86L (TTT→CTT)	b1423	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/0); new base C (65/68); total (68/68)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 142 total aligned reads displayed.

AATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATTGAT  >  NC_000913/1493477‑1493656

aaTGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGctt                                                                                         <  1:3976254/93‑1 (MQ=255)

 aTGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTg                                                                                        >  1:2109549/1‑93 (MQ=255)

 aTGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGCCGGATGcttc                                                                                        >  1:2906507/1‑92 (MQ=255)

   gCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGcc                                                                                      <  1:34725/93‑1 (MQ=255)

      gACGCGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGCCTGCCACAGAGCTTGCTACTCTTCGGCGGATGCTTGctctt                                                                                   >  1:1161698/1‑89 (MQ=255)

      gACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCAtt                                                                                   >  1:4478994/1‑93 (MQ=255)

       aCGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTa                                                                                  >  1:657191/1‑93 (MQ=255)

        cGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTAt                                                                                 <  2:2888825/93‑1 (MQ=255)

        cGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTAt                                                                                 <  2:4278009/93‑1 (MQ=255)

         gAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATg                                                                                >  1:4271081/1‑93 (MQ=255)

          aGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATgg                                                                               <  1:4402626/93‑1 (MQ=255)

          aGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATCATgg                                                                               <  1:4697680/93‑1 (MQ=255)

           gCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATggg                                                                              >  2:3714314/1‑93 (MQ=255)

           gCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATggg                                                                              <  2:2200564/93‑1 (MQ=255)

            cTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATgggg                                                                             <  1:2388361/93‑1 (MQ=255)

              ggCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGAt                                                                           >  2:1457381/1‑93 (MQ=255)

               gcgcGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATg                                                                          <  2:4374253/93‑1 (MQ=255)

                 gcgACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTa                                                                        >  2:1198015/1‑93 (MQ=255)

                 gcgACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTa                                                                        <  2:1253098/93‑1 (MQ=255)

                  cgACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAc                                                                       >  1:2917024/1‑93 (MQ=255)

                   gACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAcc                                                                      >  1:1635852/1‑93 (MQ=255)

                    aCTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAccc                                                                     <  1:887598/93‑1 (MQ=255)

                    aCTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAccc                                                                     >  1:3621351/1‑93 (MQ=255)

                       gAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCggg                                                                  >  1:3404386/1‑92 (MQ=255)

                       gAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGt                                                                  <  1:4526738/93‑1 (MQ=255)

                       gAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGt                                                                  >  2:4646670/1‑93 (MQ=255)

                        aaTGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTg                                                                 >  2:3890993/1‑93 (MQ=255)

                           gTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGc                                                              <  1:2138630/93‑1 (MQ=255)

                           gTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGc                                                              <  2:4022785/93‑1 (MQ=255)

                               aaCGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTAtt                                                          <  1:2953444/93‑1 (MQ=255)

                               aaCGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTAtt                                                          <  1:3293498/93‑1 (MQ=255)

                                aCGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTa                                                         >  1:939983/1‑93 (MQ=255)

                                 cGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTAc                                                        <  2:3494889/93‑1 (MQ=255)

                                 cGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTAc                                                        <  2:2918227/93‑1 (MQ=255)

                                 cGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTAc                                                        <  1:272497/93‑1 (MQ=255)

                                  gTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGGACCCGGTGATCTATTACg                                                       >  2:2437205/1‑93 (MQ=255)

                                   tCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGa                                                      >  2:4190571/1‑93 (MQ=255)

                                      tGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTCCCCGGTGAGCTATTACGAtct                                                   <  1:554174/93‑1 (MQ=255)

                                      tGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtct                                                   <  2:2314665/93‑1 (MQ=255)

                                      tGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtct                                                   <  1:1128400/93‑1 (MQ=255)

                                        gggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctct                                                 <  1:2012052/93‑1 (MQ=255)

                                        gggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctct                                                 <  2:3337875/93‑1 (MQ=255)

                                          gCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCg                                               >  1:686811/1‑93 (MQ=255)

                                           cGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgc                                              <  1:2220448/93‑1 (MQ=255)

                                           cGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgc                                              <  2:3017072/93‑1 (MQ=255)

                                              ttCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgcagg                                           >  2:4553407/1‑93 (MQ=255)

                                               tCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgcaggc                                          <  2:3346577/93‑1 (MQ=255)

                                                 gcgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgccggcag                                        >  2:1113179/1‑93 (MQ=255)

                                                 gcgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgcaggcag                                        <  1:4671171/93‑1 (MQ=255)

                                                 gcgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgcaggcag                                        <  2:2564185/93‑1 (MQ=255)

                                                 gcgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgcaggcag                                        <  1:4104987/93‑1 (MQ=255)

                                                  cgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGcaggcagg                                       >  2:1870787/1‑93 (MQ=255)

                                                   gTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAggg                                      <  2:2050246/93‑1 (MQ=255)

                                                   gTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAggg                                      <  2:4139724/93‑1 (MQ=255)

                                                   gTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAggg                                      <  1:2181068/93‑1 (MQ=255)

                                                   gTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAggg                                      <  2:4129841/93‑1 (MQ=255)

                                                    ttGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAgggg                                     <  2:1302772/93‑1 (MQ=255)

                                                    gttggACTGCACAAGAGCTTGCTACTCTGCGGCGGATTCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAgggg                                     <  1:506246/90‑1 (MQ=255)

                                                     tGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAggggc                                    >  2:1417406/1‑92 (MQ=255)

                                                      gggACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTg                                   >  2:4185376/1‑93 (MQ=255)

                                                      gggACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTg                                   >  1:2091906/1‑93 (MQ=255)

                                                       ggACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGc                                  <  1:2330565/93‑1 (MQ=255)

                                                         aCTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCg                                >  2:1088617/1‑93 (MQ=255)

                                                         aCTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCg                                >  1:542999/1‑93 (MQ=255)

                                                         aCTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCg                                >  2:3423677/1‑93 (MQ=255)

                                                          cTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCgg                               >  1:3391242/1‑93 (MQ=255)

                                                            gCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTa                             >  1:1515799/1‑93 (MQ=255)

                                                            gCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTa                             >  2:273620/1‑93 (MQ=255)

                                                             cacaAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTAc                            >  2:57581/1‑93 (MQ=255)

                                                             cacaAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTAc                            <  2:1714160/93‑1 (MQ=255)

                                                              acaAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACa                           <  1:2514357/93‑1 (MQ=255)

                                                              acaAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACa                           <  1:4647121/93‑1 (MQ=255)

                                                                aaGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACAtt                         >  1:3079949/1‑93 (MQ=255)

                                                                 agagCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTc                        >  1:4443956/1‑93 (MQ=255)

                                                                  gagCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCg                       <  1:3598581/93‑1 (MQ=255)

                                                                  gagCTTGCTACTCTGCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCg                       <  2:3655278/93‑1 (MQ=255)

                                                                      ttGCTACTCTTCGGCGGATGCTTGCCATTATGGGGGTGTACCCGGTGAGCTATTACGATCTCTCGCACGCAGGGGTGCCGGTACATTCGACAg                   >  2:361320/1‑93 (MQ=255)

                                                                      ttGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAg                   <  2:249044/93‑1 (MQ=255)

                                                                      ttGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAg                   >  1:637428/1‑93 (MQ=255)

                                                                        gCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCa                 >  1:468918/1‑93 (MQ=255)

                                                                        gCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCa                 >  1:3077588/1‑93 (MQ=255)

                                                                        gCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCa                 <  2:3966081/93‑1 (MQ=255)

                                                                         cTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAt                >  1:3977156/1‑93 (MQ=255)

                                                                          tACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAtt               >  2:2325659/1‑93 (MQ=255)

                                                                          tACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAtt               <  1:876078/93‑1 (MQ=255)

                                                                          tACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAtt               >  2:2325660/1‑93 (MQ=255)

                                                                            ctctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGGGCCGGTACATTCGACAGCATTTc             >  1:3533232/1‑93 (MQ=255)

                                                                              ctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGGGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCgg           >  1:1993260/1‑93 (MQ=255)

                                                                               ttCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGc          >  1:4017448/1‑93 (MQ=255)

                                                                               ttCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGc          <  1:2454970/93‑1 (MQ=255)

                                                                                tCGGCGGATGCTTGTCATTATGGGTATGTACCCGGTGAGCTATTACGACCTCACGCAGAAGGGGGTGCCGGTACATACGACAGCATTTCggac         <  2:2606098/93‑3 (MQ=255)

                                                                                tCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGcc         <  2:4724156/93‑1 (MQ=255)

                                                                                 cggcggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGccc        >  1:629143/1‑93 (MQ=255)

                                                                                 cggcggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGccc        >  2:433501/1‑93 (MQ=255)

                                                                                  ggcggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCa       >  2:3639082/1‑93 (MQ=255)

                                                                                   gcggATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTTGGCCCAt      >  2:4037085/1‑93 (MQ=255)

                                                                                   gcggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCAt      <  2:2285709/93‑1 (MQ=255)

                                                                                    cggATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCAtt     >  1:3988538/1‑93 (MQ=255)

                                                                                      gATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtga   >  2:2243887/1‑93 (MQ=255)

                                                                                       aTGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtgat  >  1:1038886/1‑93 (MQ=255)

AATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATTGAT  >  NC_000913/1493477‑1493656

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: `‑`