BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	431,405	C→T	intergenic (‑168/+131)	tsx ← / ← yajI	Tsx/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	431,405	0	C	T	99.6%	817.6 / ‑5.5	221	intergenic (‑168/+131)	tsx/yajI	Tsx/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (162/58); total (162/59)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

Only 100 of 231 total aligned reads displayed.

GAGCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGG  >  NC_000913/431316‑431494

gagCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATAtt                                                                                        <  2:2116547/93‑1 (MQ=255)

gagCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATAtt                                                                                        <  1:3804031/93‑1 (MQ=255)

gagCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATAtt                                                                                        >  1:2323290/1‑93 (MQ=255)

 agCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTc                                                                                       <  1:3098783/93‑1 (MQ=255)

 agCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTc                                                                                       >  2:419297/1‑93 (MQ=255)

  gCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTcg                                                                                      >  1:1295102/1‑93 (MQ=255)

   ccTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTcgg                                                                                     <  1:2719757/93‑1 (MQ=255)

    cTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTcggg                                                                                    >  2:2801950/1‑93 (MQ=255)

     tattatAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTcgggc                                                                                   >  2:3256863/1‑93 (MQ=255)

       ttatAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTcgggcgg                                                                                 >  1:831522/1‑93 (MQ=255)

            aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGAtt                                                                            >  1:3531120/1‑93 (MQ=255)

            aaaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGAtt                                                                            >  2:2244950/1‑93 (MQ=255)

             aaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTg                                                                           <  1:4490752/93‑1 (MQ=255)

             aaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTg                                                                           >  1:2457483/1‑93 (MQ=255)

             aaaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTg                                                                           >  1:4244087/1‑93 (MQ=255)

              aaTCATTGCTCACAAAAGTATGTTTCGTTTCAAAGTTCTCCCATTTACGTAATTTATTGCTTTCACGTTTACACATATTGGGGGGGGGATTGt                                                                          >  2:3647010/1‑93 (MQ=255)

              aaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGt                                                                          >  2:2481170/1‑93 (MQ=255)

              aaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGt                                                                          >  2:558162/1‑93 (MQ=255)

              aaTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGt                                                                          >  1:434274/1‑93 (MQ=255)

               aTCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTa                                                                         >  2:2280814/1‑93 (MQ=255)

                  aTTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTg                                                                      >  2:4303961/1‑93 (MQ=255)

                  aTTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTg                                                                      >  2:1490912/1‑93 (MQ=255)

                     gCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACt                                                                   >  1:2994096/1‑93 (MQ=255)

                       tCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACttt                                                                 >  2:4244869/1‑93 (MQ=255)

                       tCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACttt                                                                 >  2:3862284/1‑93 (MQ=255)

                         acaAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCg                                                               >  1:3263721/1‑93 (MQ=255)

                         acaAAAATATGTTTCGTTTCACAGTTCTATCATTAACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCg                                                               >  1:2219668/1‑93 (MQ=255)

                          caAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGc                                                              >  2:1191695/1‑93 (MQ=255)

                           aaaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGcc                                                             >  1:218546/1‑93 (MQ=255)

                            aaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGccc                                                            >  1:2933305/1‑93 (MQ=255)

                            aaaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGCTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGccc                                                            >  2:2221049/1‑93 (MQ=255)

                             aaaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCt                                                           >  1:588874/1‑93 (MQ=255)

                              aaTATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTTGCCCtt                                                          >  2:2619835/1‑93 (MQ=255)

                                tatGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCttgt                                                        >  2:3427086/1‑93 (MQ=255)

                                tatGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCttgt                                                        >  2:18093/1‑93 (MQ=255)

                                 atGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCttgtt                                                       >  2:2410730/1‑93 (MQ=255)

                                 atGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCttgtt                                                       <  2:4653063/93‑1 (MQ=255)

                                  tGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCttgttg                                                      >  2:1746457/1‑93 (MQ=255)

                                    tttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCa                                                    >  2:2591189/1‑93 (MQ=255)

                                    tttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCa                                                    >  2:3092568/1‑93 (MQ=255)

                                    tttcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCa                                                    >  1:4588162/1‑93 (MQ=255)

                                      tcgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCaaa                                                  >  2:2962468/1‑93 (MQ=255)

                                       cgtttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAAc                                                 <  2:3835248/93‑1 (MQ=255)

                                          ttcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACttt                                              <  2:1018705/93‑1 (MQ=255)

                                           tcACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACtttt                                             >  1:1000161/1‑93 (MQ=255)

                                            cacaGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTa                                            <  2:3489205/93‑1 (MQ=255)

                                             acaGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTAc                                           <  1:1622310/93‑1 (MQ=255)

                                               aGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTa                                         >  2:3728321/1‑93 (MQ=255)

                                               aGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTa                                         >  1:2006381/1‑93 (MQ=255)

                                                gTTCTATCATTTACGTAATCGATTGCGTTCGCGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTAt                                        <  2:122920/93‑1 (MQ=255)

                                                gTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTAt                                        >  2:1350328/1‑93 (MQ=255)

                                                 ttCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATc                                       >  2:2434460/1‑93 (MQ=255)

                                                  tCTATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCa                                      >  2:2993954/1‑93 (MQ=255)

                                                    tATCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAAt                                    >  2:3142908/1‑93 (MQ=255)

                                                     aTCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATc                                   >  1:3626080/1‑93 (MQ=255)

                                                      tCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATcc                                  <  2:1904738/93‑1 (MQ=255)

                                                      tCATTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATcc                                  >  1:3374331/1‑93 (MQ=255)

                                                        aTTTACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAg                                >  2:579551/1‑93 (MQ=255)

                                                         tttACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                               >  1:1137397/1‑93 (MQ=255)

                                                         tttACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                               <  2:1583823/93‑1 (MQ=255)

                                                         tttACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                               >  1:4000769/1‑93 (MQ=255)

                                                         tttACGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAgg                               >  2:4369076/1‑93 (MQ=255)

                                                             cGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGAt                           <  1:47966/93‑1 (MQ=255)

                                                             cGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGAt                           <  2:3544507/93‑1 (MQ=255)

                                                             cGTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGAt                           <  2:4379844/93‑1 (MQ=255)

                                                              gTAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATg                          >  2:141179/1‑93 (MQ=255)

                                                               tAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGt                         >  2:3765268/1‑93 (MQ=255)

                                                               tAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGt                         >  1:2645627/1‑93 (MQ=255)

                                                               tAATCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGt                         >  1:1573945/1‑93 (MQ=255)

                                                                aaTCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTa                        <  1:4455959/93‑1 (MQ=255)

                                                                aaTCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTa                        <  2:3984878/93‑1 (MQ=255)

                                                                aaTCGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTa                        >  2:4129755/1‑93 (MQ=255)

                                                                   cGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGc                     <  2:727310/93‑1 (MQ=255)

                                                                   cGATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGc                     >  2:3366149/1‑93 (MQ=255)

                                                                    gATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCa                    >  2:1726369/1‑93 (MQ=255)

                                                                    gATTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCa                    >  1:3863641/1‑93 (MQ=255)

                                                                     aTTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCag                   >  1:2396098/1‑93 (MQ=255)

                                                                     aTTGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCag                   >  1:786577/1‑93 (MQ=255)

                                                                      ttGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCaga                  >  1:2255003/1‑93 (MQ=255)

                                                                      ttGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCaga                  <  1:1978973/93‑1 (MQ=255)

                                                                      ttGCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCaga                  >  2:2028363/1‑93 (MQ=255)

                                                                        gCGTTCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGt                >  1:4144247/1‑93 (MQ=255)

                                                                           ttCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAAc             >  1:703504/1‑93 (MQ=255)

                                                                           ttCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAAc             >  2:3603817/1‑93 (MQ=255)

                                                                           ttCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAAc             <  2:2857009/93‑1 (MQ=255)

                                                                            tCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACa            >  1:2623882/1‑93 (MQ=255)

                                                                            tCACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACa            >  2:631792/1‑93 (MQ=255)

                                                                             cACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaa           >  1:1432306/1‑93 (MQ=255)

                                                                             cACGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaa           <  1:1731520/93‑1 (MQ=255)

                                                                              aCGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaa          <  2:1116540/93‑1 (MQ=255)

                                                                               cGTTTACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACaaaa         >  1:2438476/1‑93 (MQ=255)

                                                                                 tttACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATg       >  2:389918/1‑93 (MQ=255)

                                                                                  ttACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGa      <  2:2722376/93‑1 (MQ=255)

                                                                                  ttACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGa      <  2:4272081/93‑1 (MQ=255)

                                                                                  ttACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGa      >  1:3694119/1‑93 (MQ=255)

                                                                                  ttACACATATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGa      >  2:2663861/1‑93 (MQ=255)

                                                                                     cacataTTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACgg   >  2:513843/1‑93 (MQ=255)

                                                                                     cacataTTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACgg   >  1:1246133/1‑93 (MQ=255)

                                                                                     cacataTTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACgg   >  2:841323/1‑93 (MQ=255)

                                                                                      acataTTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACggg  >  2:4068400/1‑93 (MQ=255)

GAGCCTATTATAAAAATCATTGCTCACAAAAATATGTTTCGTTTCACAGTTCTATCATTTACGTAATCGATTGCGTTCACGTTTACACACATTCGGGCGGGGATTGTACTGACTTTCGCCCTTGTTGCAAACTTTTACTATCAATCCAGGTGATGTAAGCAGAGTAACAAAATGACGGG  >  NC_000913/431316‑431494

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