BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,572,272	G→A	G190E (GGG→GAG)	yjiU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,572,272	0	G	A	85.0%	981.5 / 111.2	353	G190E (GGG→GAG)	yjiU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (27/26); new base A (147/153); total (174/179)
Fisher's exact test for biased strand distribution p-value = 8.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 361 total aligned reads displayed.

CTTTACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCAAAGG  >  NC_000913/4572180‑4572363

  ttACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAgaga                                                                                           <  2:5762727/93‑1 (MQ=255)

  ttACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAgaga                                                                                           <  2:6447947/93‑1 (MQ=255)

   tACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGAc                                                                                          <  2:6065612/93‑1 (MQ=255)

   tACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGAc                                                                                          <  1:4942390/93‑1 (MQ=255)

   tACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGAc                                                                                          <  2:420990/93‑1 (MQ=255)

   tACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGAc                                                                                          >  1:5493110/1‑93 (MQ=255)

   tACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGAc                                                                                          <  2:6246676/93‑1 (MQ=255)

    aCCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACa                                                                                         >  1:1825509/1‑93 (MQ=255)

     ccAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGACAg                                                                                        >  1:2069399/1‑93 (MQ=255)

     ccAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAg                                                                                        >  2:1396367/1‑93 (MQ=255)

      cAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGc                                                                                       <  2:6530414/93‑1 (MQ=255)

      cAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGc                                                                                       >  2:1619367/1‑93 (MQ=255)

       aGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGcc                                                                                      <  2:244048/93‑1 (MQ=255)

            tgaAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAAt                                                                                 <  1:3776411/93‑1 (MQ=255)

             gaAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATg                                                                                <  1:193136/93‑1 (MQ=255)

             gaAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATg                                                                                <  2:2181442/93‑1 (MQ=255)

               aGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGc                                                                              >  1:4180250/1‑93 (MQ=255)

                  cAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGt                                                                           <  2:5345546/93‑1 (MQ=255)

                   aaaGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTGCAGACAGTCGACAAGAGACAGCCCGAATGGCGGTg                                                                          >  2:4399012/1‑93 (MQ=255)

                   aaaGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTg                                                                          >  2:6841371/1‑93 (MQ=255)

                   aaaGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTg                                                                          >  1:6304777/1‑93 (MQ=255)

                    aaGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGc                                                                         <  1:7120100/93‑1 (MQ=255)

                      gagGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGt                                                                       >  1:2707910/1‑93 (MQ=255)

                       agGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTa                                                                      <  2:5428883/93‑1 (MQ=255)

                        gggCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTAt                                                                     <  1:874865/93‑1 (MQ=255)

                           cGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCGCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATgca                                                                  <  2:2849146/93‑1 (MQ=255)

                           cGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATgca                                                                  <  1:3076399/93‑1 (MQ=255)

                                  ccAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTGCAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaagaa                                                           >  1:7184753/1‑93 (MQ=255)

                                  ccAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaagaa                                                           >  2:6029862/1‑93 (MQ=255)

                                  ccAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaagaa                                                           <  2:4948973/93‑1 (MQ=255)

                                  ccAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaagaa                                                           <  2:3565836/93‑1 (MQ=255)

                                    aGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGc                                                         >  2:5224960/1‑93 (MQ=255)

                                    aGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGc                                                         <  1:5619178/93‑1 (MQ=255)

                                     gATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCa                                                        <  1:4261581/93‑1 (MQ=255)

                                      aTTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAg                                                       <  2:6365852/93‑1 (MQ=255)

                                      aTTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAg                                                       >  1:172825/1‑93 (MQ=255)

                                        tCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGaa                                                     <  1:5104251/93‑1 (MQ=255)

                                         cGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAAc                                                    <  1:7356820/93‑1 (MQ=255)

                                         cGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAAc                                                    <  2:1791225/93‑1 (MQ=255)

                                            tGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCa                                                 <  2:4571240/93‑1 (MQ=255)

                                            tGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCa                                                 <  1:1568924/93‑1 (MQ=255)

                                            tGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCa                                                 <  2:4793660/93‑1 (MQ=255)

                                            tGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCa                                                 <  1:4669812/93‑1 (MQ=255)

                                             gCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAg                                                >  2:5543219/1‑93 (MQ=255)

                                             gCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAg                                                <  2:169089/93‑1 (MQ=255)

                                             gCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAg                                                <  1:6129024/93‑1 (MQ=255)

                                              ccAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGt                                               <  2:2276407/93‑1 (MQ=255)

                                               cAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTc                                              <  1:3218359/93‑1 (MQ=255)

                                                 ggACGAGGTCGTCGACTCATCACCTTTGCAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGc                                            >  2:6725366/1‑93 (MQ=255)

                                                 ggACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGc                                            >  1:1337420/1‑93 (MQ=255)

                                                 ggACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGc                                            >  2:936177/1‑93 (MQ=255)

                                                  gACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCt                                           >  2:4668832/1‑93 (MQ=255)

                                                   aCGAGGTCGTCGACTCATCACCTTTGCAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCtt                                          >  2:6999539/1‑93 (MQ=255)

                                                    cGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTa                                         <  2:3134583/93‑1 (MQ=255)

                                                      aGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTAcg                                       >  2:5842630/1‑93 (MQ=255)

                                                      aGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTAcg                                       <  1:3524484/93‑1 (MQ=255)

                                                       ggTCGTCGACTCATCACCTTTGCAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTAcgc                                      <  1:6415875/93‑1 (MQ=255)

                                                       ggTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTAcgc                                      >  1:619794/1‑93 (MQ=255)

                                                        gtcgtcGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTAcgcg                                     <  2:6606672/93‑1 (MQ=255)

                                                         tcgtcgACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCgg                                    >  1:4785433/1‑93 (MQ=255)

                                                          cgtcgACTCATCACCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGa                                   <  1:2860645/93‑1 (MQ=255)

                                                                cTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTa                             <  1:5984271/93‑1 (MQ=255)

                                                                cTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTa                             >  1:4538014/1‑93 (MQ=255)

                                                                 tcatcaCCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAg                            <  1:4391906/93‑1 (MQ=255)

                                                                 tcatcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAg                            <  1:357742/93‑1 (MQ=255)

                                                                 tcatcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAg                            >  1:3267302/1‑93 (MQ=255)

                                                                 tcatcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAg                            >  1:6955341/1‑93 (MQ=255)

                                                                 tcatcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAg                            <  2:1251248/93‑1 (MQ=255)

                                                                  catcaCCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGt                           <  2:3415622/93‑1 (MQ=255)

                                                                  catcaCCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGt                           >  2:5904080/1‑93 (MQ=255)

                                                                  catcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGt                           >  1:6644879/1‑93 (MQ=255)

                                                                  catcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGt                           >  1:1141331/1‑93 (MQ=255)

                                                                     caCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGa                        >  2:7490849/1‑93 (MQ=255)

                                                                       ccTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGaaa                      >  1:6063356/1‑93 (MQ=255)

                                                                        cTTTGCAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAAt                     >  2:4160063/1‑93 (MQ=255)

                                                                          ttACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATtc                   <  1:6826403/93‑1 (MQ=255)

                                                                           tACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATtct                  <  2:3899104/93‑1 (MQ=255)

                                                                           tACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATtct                  >  2:3813155/1‑93 (MQ=255)

                                                                           tACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATtct                  >  1:5784458/1‑93 (MQ=255)

                                                                            acagacagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATtctc                 >  1:1036569/1‑93 (MQ=255)

                                                                              agacagTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAg               >  1:2459874/1‑93 (MQ=255)

                                                                               gacagTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGc              >  1:2272496/1‑93 (MQ=255)

                                                                               gacagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGc              >  1:3907099/1‑93 (MQ=255)

                                                                               gacagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGc              >  1:1411135/1‑93 (MQ=255)

                                                                                 cagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTg            >  1:506704/1‑93 (MQ=255)

                                                                                 cagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTg            <  2:6328731/93‑1 (MQ=255)

                                                                                  agTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTgg           <  2:7287483/93‑1 (MQ=255)

                                                                                   gTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGc          >  2:7125893/1‑93 (MQ=255)

                                                                                   gTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGc          >  1:6717090/1‑93 (MQ=255)

                                                                                   gTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGc          >  1:5681361/1‑93 (MQ=255)

                                                                                   gTCGACAAGAGACAGCCCGAATGGCGGGGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGc          >  1:2369157/1‑93 (MQ=255)

                                                                                     cGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCAt        <  2:5086196/93‑1 (MQ=255)

                                                                                       aCAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCa      >  2:1078632/1‑93 (MQ=255)

                                                                                       aCAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCa      >  1:2588158/1‑93 (MQ=255)

                                                                                       aCAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCa      >  1:7228740/1‑93 (MQ=255)

                                                                                        cAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCaa     <  2:2704677/93‑1 (MQ=255)

                                                                                         aaGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCaaa    <  1:1906484/93‑1 (MQ=255)

                                                                                         aaGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCaaa    <  2:6955341/93‑1 (MQ=255)

                                                                                          aGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCAAAg   <  1:2873057/93‑1 (MQ=255)

                                                                                          aGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCAAAg   >  2:5131624/1‑93 (MQ=255)

CTTTACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCAAAGG  >  NC_000913/4572180‑4572363

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