BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,025,603	T→C	G154G (GGT→GGC)	ygfQ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,025,603	0	T	C	95.1%	408.0 / 6.7	142	G154G (GGT→GGC)	ygfQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/4); new base C (47/88); total (50/92)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 169 total aligned reads displayed.

CGTCATCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑CAGTGGCGCTCGGTGCG  >  NC_000913/3025513‑3025692

cGTCATCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAggtatcggcat                                                                                          <  2:4382067/93‑1 (MQ=255)

 gTCATCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAggtatcggcatc                                                                                         >  2:444904/1‑93 (MQ=255)

 gTCATCTTCACAGCCATTTCCGTAACCGGTGTGCGTATCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAggtatcggcatc                                                                                         <  1:1199811/93‑1 (MQ=255)

  tCATCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGgtatcggcatcg                                                                                        >  1:907695/1‑93 (MQ=255)

  tCATCTTCACCGCCATTTCCATAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGgtatcggcatcg                                                                                        >  2:4571747/1‑93 (MQ=255)

   cATCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGtatcggcatcgg                                                                                       <  1:2350338/93‑1 (MQ=255)

    aTCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCggg                                                                                      <  1:1070858/93‑1 (MQ=255)

    aTCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCggg                                                                                      <  1:995015/93‑1 (MQ=255)

     tCTTCGCCGCCATTTCCGTAGCCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGc                                                                                     <  2:3047644/93‑1 (MQ=255)

     tCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGTGTATCGCTCACGGTACAGGTATCGGCATCGGGc                                                                                     <  2:1966517/93‑1 (MQ=255)

     tCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGTGTATCGCTCACGGTACAGGTATCGGCATCGGGc                                                                                     <  2:1966518/93‑1 (MQ=255)

     tCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGc                                                                                     >  2:2570715/1‑93 (MQ=255)

     tCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGc                                                                                     >  2:3981648/1‑93 (MQ=255)

       ttCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTg                                                                                   >  1:1695456/1‑93 (MQ=255)

        tCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGt                                                                                  >  1:3646225/1‑93 (MQ=255)

        tCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGt                                                                                  >  2:2833039/1‑93 (MQ=255)

          aCCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGttt                                                                                >  2:1702215/1‑93 (MQ=255)

           ccGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTc                                                                               >  2:2984464/1‑93 (MQ=255)

           ccGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTc                                                                               <  1:3649286/93‑1 (MQ=255)

             gCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTctg                                                                             <  2:3603172/93‑1 (MQ=255)

               cATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTctgct                                                                           >  2:383876/1‑93 (MQ=255)

                aTTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGTATCGGGCTGTTTctgctg                                                                          >  1:2154169/1‑93 (MQ=255)

                aTTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTctgctg                                                                          >  1:1735928/1‑93 (MQ=255)

                aTTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTCctccgc                                                                          >  1:2461236/1‑89 (MQ=255)

                 tttCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTctgctgc                                                                         <  1:110285/93‑1 (MQ=255)

                 tttCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTctgctgc                                                                         <  1:1739038/93‑1 (MQ=255)

                  ttCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTctgctgct                                                                        <  1:3473436/93‑1 (MQ=255)

                    ccGTAACCGGTGTGCGTACCTGGCTCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGa                                                                      <  2:2590554/93‑1 (MQ=255)

                    ccGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGa                                                                      <  2:2157969/93‑1 (MQ=255)

                    ccGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGa                                                                      <  1:1964129/93‑1 (MQ=255)

                     cGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGAt                                                                     <  2:2981273/93‑1 (MQ=255)

                      gTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGAtt                                                                    <  1:483936/93‑1 (MQ=255)

                      gTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGAtt                                                                    >  1:3106699/1‑93 (MQ=255)

                         aCCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATtgct                                                                 <  2:3183519/93‑1 (MQ=255)

                         aCCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATtgct                                                                 >  2:1840917/1‑93 (MQ=255)

                           cGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATtgctgc                                                               <  1:2689398/93‑1 (MQ=255)

                            ggTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTgctgct                                                              >  1:1843718/1‑93 (MQ=255)

                            ggTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTgctgct                                                              >  2:3907165/1‑93 (MQ=255)

                             gtgtGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTa                                                             >  1:4737006/1‑93 (MQ=255)

                              tgtgCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTaa                                                            <  1:1540318/93‑1 (MQ=255)

                              ggtgCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTaa                                                            <  1:2354978/92‑1 (MQ=255)

                               gtgCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAAc                                                           <  1:589265/93‑1 (MQ=255)

                                tgCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACg                                                          <  1:2802836/93‑1 (MQ=255)

                                 gCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACgg                                                         <  1:3161581/93‑1 (MQ=255)

                                  cGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGTATCGGCCTGTTTCTGCTGCTGATTGCTGCTAACGgt                                                        <  1:4116077/93‑1 (MQ=255)

                                  cGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGgt                                                        <  1:3881874/93‑1 (MQ=255)

                                   gTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGgtg                                                       >  1:749021/1‑93 (MQ=255)

                                    tACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGgtgt                                                      <  1:4237801/93‑1 (MQ=255)

                                    tACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGgtgt                                                      <  2:3895654/93‑1 (MQ=255)

                                        tGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGt                                                  <  2:3710740/93‑1 (MQ=255)

                                        tGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGt                                                  <  2:4114509/93‑1 (MQ=255)

                                        tGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGt                                                  <  2:3316078/93‑1 (MQ=255)

                                         ggATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTa                                                 >  2:3995034/1‑93 (MQ=255)

                                          gATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTAt                                                <  1:1040157/93‑1 (MQ=255)

                                          gATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTAt                                                <  2:2745906/93‑1 (MQ=255)

                                          gATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTAt                                                <  2:4502738/93‑1 (MQ=255)

                                           aTCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATg                                               >  2:616987/1‑93 (MQ=255)

                                           aTCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATg                                               >  2:4038652/1‑93 (MQ=255)

                                           aTCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATg                                               <  1:4134056/93‑1 (MQ=255)

                                             cTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGt                                             <  1:578364/93‑1 (MQ=255)

                                               tACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTa                                           <  2:1014984/93‑1 (MQ=255)

                                                aCGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTAt                                          <  2:4246598/93‑1 (MQ=255)

                                                   tAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCaa                                       >  2:2054108/1‑93 (MQ=255)

                                                    aaTTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCaca                                      >  1:1667807/1‑91 (MQ=255)

                                                    aaTTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCaaa                                      >  1:244671/1‑93 (MQ=255)

                                                      tttGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCaaaaa                                    <  1:1492758/93‑1 (MQ=255)

                                                       ttGCCGATGGGTATCGCTCGCGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAacac                                   <  2:4212959/93‑1 (MQ=255)

                                                       ttGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAAc                                   <  1:3267504/93‑1 (MQ=255)

                                                        tGCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAAcc                                  <  2:3374913/93‑1 (MQ=255)

                                                         gCCGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAAccc                                 <  2:821550/93‑1 (MQ=255)

                                                           cGATGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGa                               >  2:526509/1‑93 (MQ=255)

                                                             aTGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGAtt                             >  2:1303403/1‑93 (MQ=255)

                                                             aTGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGAtt                             <  2:4635536/93‑1 (MQ=255)

                                                              tGGGTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTg                            >  2:3418069/1‑93 (MQ=255)

                                                               gggTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGa                           <  1:2368694/93‑1 (MQ=255)

                                                                ggTATCGCTCACGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGaa                          >  2:4104002/1‑93 (MQ=255)

                                                                ggTATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGaa                          <  2:1639852/93‑1 (MQ=255)

                                                                  tATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAgg                        <  2:1699761/93‑1 (MQ=255)

                                                                  tATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAgg                        <  1:2420806/93‑1 (MQ=255)

                                                                  tATCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAAACCGATTGAAgg                        >  2:2486012/1‑93 (MQ=255)

                                                                    tCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCt                      >  1:2973257/1‑93 (MQ=255)

                                                                    tCGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCt                      <  1:3388752/93‑1 (MQ=255)

                                                                     cGCTCACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCtt                     <  2:3878555/93‑1 (MQ=255)

                                                                         cACGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑cca                <  1:1178276/93‑3 (MQ=255)

                                                                         cACGGTACAGGGATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGTTATGGTTATCAAAAACCCGATTGAAGGCTTG‑cca                <  1:3319350/93‑3 (MQ=255)

                                                                          aCGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑ccag               >  1:270651/1‑90 (MQ=255)

                                                                           cGGTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑ccagt              <  2:1059524/93‑5 (MQ=255)

                                                                             gTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTgg             >  1:2149415/1‑93 (MQ=255)

                                                                             gTACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTgg             >  2:1341778/1‑93 (MQ=255)

                                                                              tACAGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGgc            <  2:1352379/93‑1 (MQ=255)

                                                                               tcAGGTATCGACATCGAGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAAACTTGCCAGTGgcg           <  2:4363013/92‑1 (MQ=255)

                                                                                 agatatcggcatcggACTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCt         <  2:2282308/90‑1 (MQ=255)

                                                                                 aGGTATCGGCATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCt         <  2:2282307/93‑1 (MQ=255)

                                                                                     atcggcatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCggg     >  2:879529/1‑92 (MQ=255)

                                                                                     atcggcatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGt     <  2:2368754/93‑1 (MQ=255)

                                                                                      tcggcatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTg    <  2:2769924/93‑1 (MQ=255)

                                                                                      tcggcatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTg    >  2:3841116/1‑93 (MQ=255)

                                                                                      tcggcatcggGCTGTTGCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGTCAGTGGCGCTCGGTg    >  2:4205680/1‑93 (MQ=255)

                                                                                       cggcatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCTATTGAAGGCTTGCCAGTGGCGCTCGGTGc   <  2:3233778/93‑1 (MQ=255)

                                                                                        ggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCg  <  1:2079844/93‑1 (MQ=255)

CGTCATCTTCACCGCCATTTCCGTAACCGGTGTGCGTACCTGGATCTTACGTAATTTGCCGATGGGTATCGCTCACGGTACAGGTATCGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTG‑CAGTGGCGCTCGGTGCG  >  NC_000913/3025513‑3025692

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