BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,938,723	C→T	R28R (CGG→CGA)	yiePs ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,938,723	0	C	T	99.5%	2150.4 / ‑5.4	584	R28R (CGG→CGA)	yiePs	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/3); new base T (316/265); total (316/268)
Fisher's exact test for biased strand distribution p-value = 9.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 598 total aligned reads displayed.

AAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAA  >  NC_000913/3938631‑3938812

  gTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCg                                                                                         >  1:5717950/1‑93 (MQ=255)

   tGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGa                                                                                        >  1:4705334/1‑93 (MQ=255)

      tCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTa                                                                                     <  1:664603/93‑1 (MQ=255)

        gATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTAcc                                                                                   >  2:9467083/1‑93 (MQ=255)

        gATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTAcc                                                                                   >  2:7588848/1‑93 (MQ=255)

         aTGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCa                                                                                  >  2:5224868/1‑93 (MQ=255)

           gACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAAt                                                                                <  1:1137961/93‑1 (MQ=255)

              cTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcg                                                                             >  1:931898/1‑93 (MQ=255)

               tgatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgc                                                                            >  1:2849789/1‑93 (MQ=255)

               tgatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgc                                                                            >  1:3747270/1‑93 (MQ=255)

                ggtgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcg                                                                           <  2:8522200/91‑1 (MQ=255)

                gatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcg                                                                           >  1:2686532/1‑93 (MQ=255)

                gatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcg                                                                           >  2:2182785/1‑93 (MQ=255)

                 atgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcgg                                                                          <  1:8779030/93‑1 (MQ=255)

                  tgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggt                                                                         >  1:867484/1‑93 (MQ=255)

                  tgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggt                                                                         >  1:6500978/1‑93 (MQ=255)

                  tgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCCATTGTGGCATGACTCGAGTACCAATtcgcggt                                                                         <  2:8415008/93‑1 (MQ=255)

                   gaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggtc                                                                        >  2:6633062/1‑93 (MQ=255)

                       gTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggtcgcgg                                                                    >  2:3363849/1‑93 (MQ=255)

                         tctcTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTa                                                                  <  1:5331834/93‑1 (MQ=255)

                         tctcTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTa                                                                  >  1:4452631/1‑93 (MQ=255)

                         tctcTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTa                                                                  >  2:411615/1‑93 (MQ=255)

                          ctctTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaa                                                                 >  2:4622941/1‑93 (MQ=255)

                          ctctTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaa                                                                 <  1:3991862/93‑1 (MQ=255)

                          ctctTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaa                                                                 >  1:3002420/1‑93 (MQ=255)

                            cttctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaaaa                                                               >  1:6095040/1‑93 (MQ=255)

                            cttctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaaaa                                                               >  2:9593421/1‑93 (MQ=255)

                             ttctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAAc                                                              >  2:3573185/1‑93 (MQ=255)

                              tctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAAcc                                                             <  1:8261931/93‑1 (MQ=255)

                                tGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCAt                                                           >  1:8005516/1‑93 (MQ=255)

                                 gTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATc                                                          <  2:9383114/93‑1 (MQ=255)

                                   cATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATccc                                                        <  2:8197420/93‑1 (MQ=255)

                                   cATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATccc                                                        >  1:2285911/1‑93 (MQ=255)

                                    aTCCCCCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCt                                                       <  2:1909522/93‑1 (MQ=255)

                                    aTCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCt                                                       <  1:3466581/93‑1 (MQ=255)

                                     tCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCtt                                                      >  2:2309940/1‑93 (MQ=255)

                                     tCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCtt                                                      <  1:4347587/93‑1 (MQ=255)

                                      ccaccaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCttt                                                     <  2:8037599/93‑1 (MQ=255)

                                      ccaccaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCttt                                                     <  1:7666733/93‑1 (MQ=255)

                                        accaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTg                                                   >  1:4507285/1‑93 (MQ=255)

                                        accaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTg                                                   >  1:10005800/1‑93 (MQ=255)

                                         ccaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGc                                                  <  2:7792630/93‑1 (MQ=255)

                                          caGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGcc                                                 <  1:5545900/93‑1 (MQ=255)

                                           aGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCg                                                <  2:6405699/93‑1 (MQ=255)

                                             gTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGtt                                              <  1:1556637/93‑1 (MQ=255)

                                             gTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGtt                                              >  1:6296704/1‑93 (MQ=255)

                                              tAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTa                                             <  1:4613549/93‑1 (MQ=255)

                                              tAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTa                                             >  1:2184294/1‑93 (MQ=255)

                                               aaGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTaa                                            >  2:6729944/1‑93 (MQ=255)

                                                aGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAAc                                           >  2:6142823/1‑93 (MQ=255)

                                                  cAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGt                                         <  1:9770282/93‑1 (MQ=255)

                                                  cAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGt                                         >  2:5492127/1‑93 (MQ=255)

                                                   aaTTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGtt                                        >  1:6949301/1‑93 (MQ=255)

                                                    aTTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGttt                                       <  2:486297/93‑1 (MQ=255)

                                                    aTTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGttt                                       <  1:362294/93‑1 (MQ=255)

                                                     ttCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGTCGTTAACGtttt                                      >  2:7504668/1‑93 (MQ=255)

                                                     ttCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGtttt                                      <  2:8798381/93‑1 (MQ=255)

                                                     ttCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGtttt                                      >  2:6953190/1‑93 (MQ=255)

                                                     ttCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGCTAACGtttt                                      >  1:1869979/1‑93 (MQ=255)

                                                      tCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTg                                     <  2:3010159/93‑1 (MQ=255)

                                                      tCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTg                                     >  2:6536529/1‑93 (MQ=255)

                                                       ccTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGa                                    >  2:1193817/1‑93 (MQ=255)

                                                       ccTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGa                                    <  2:4271760/93‑1 (MQ=255)

                                                          gATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCg                                 <  2:2474167/93‑1 (MQ=255)

                                                          gATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCg                                 >  2:8664388/1‑93 (MQ=255)

                                                           aTCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGc                                <  2:9723735/93‑1 (MQ=255)

                                                            tCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCt                               >  1:4897723/1‑93 (MQ=255)

                                                            tCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCt                               <  1:4114342/93‑1 (MQ=255)

                                                             cAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCtt                              <  1:7572438/93‑1 (MQ=255)

                                                             cAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCtt                              >  2:5656040/1‑93 (MQ=255)

                                                              aaGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTc                             <  2:3785213/93‑1 (MQ=255)

                                                                gAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTcgc                           <  1:2007213/93‑1 (MQ=255)

                                                                   aaTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTa                        <  1:5948812/93‑1 (MQ=255)

                                                                    aTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAc                       >  1:9096699/1‑93 (MQ=255)

                                                                    aTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAc                       <  2:5130259/93‑1 (MQ=255)

                                                                    aTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAc                       <  1:1672620/93‑1 (MQ=255)

                                                                     ttCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAcc                      <  1:7123098/93‑1 (MQ=255)

                                                                      tCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCg                     >  2:9017936/1‑93 (MQ=255)

                                                                       ccAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGc                    <  2:1355945/93‑1 (MQ=255)

                                                                       caaGTTCAATTGTGGCATGATTAGTGAACCAATTCTCAGTCGCGGTAAAACCATCCCTTTTGTCGCTAACGCTTTGCCAGATTAGCGTATCgt                    >  2:1604039/3‑92 (MQ=255)

                                                                        cAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCt                   <  1:9529637/93‑1 (MQ=255)

                                                                          gTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTACAACCATCCCTTTTGCCGTTCACGTTTTGACCGCTTCGCGTACCGCTGt                 >  2:7217270/1‑93 (MQ=255)

                                                                          gTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGt                 <  1:2568991/93‑1 (MQ=255)

                                                                           ttCGATTGTGGCATGGCTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTa                <  1:3239293/93‑1 (MQ=255)

                                                                           ttCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTa                <  2:7276896/93‑1 (MQ=255)

                                                                            tCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTAc               <  2:5884315/93‑1 (MQ=255)

                                                                            tCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTAc               <  2:1656718/93‑1 (MQ=255)

                                                                              gATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGa             <  1:7274986/93‑1 (MQ=255)

                                                                                ttGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGact           >  2:2859341/1‑93 (MQ=255)

                                                                                ttGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACCact           >  1:8450410/1‑93 (MQ=255)

                                                                                  gtgGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactca         <  1:7425126/93‑1 (MQ=255)

                                                                                   tgGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactcac        >  2:4184857/1‑93 (MQ=255)

                                                                                    ggCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactccct       >  1:5133099/1‑93 (MQ=255)

                                                                                    ggCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactcact       >  2:9132553/1‑93 (MQ=255)

                                                                                       aTGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCa    >  2:4756050/1‑93 (MQ=255)

                                                                                        tGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCaa   <  1:7785166/93‑1 (MQ=255)

                                                                                        tGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCaa   >  2:451539/1‑93 (MQ=255)

                                                                                        tGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCaa   >  2:9810457/1‑93 (MQ=255)

                                                                                         gACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCaaa  <  2:2719398/93‑1 (MQ=255)

                                                                                         gACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCaaa  <  2:4002582/93‑1 (MQ=255)

AAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAA  >  NC_000913/3938631‑3938812

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