BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,804,224	A→G	T157T (ACT→ACC)	rfaG ←	glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,804,224	0	A	G	98.6%	673.4 / ‑1.1	213	T157T (ACT→ACC)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (117/93); total (119/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 240 total aligned reads displayed.

GCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCAGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTCGCTCGCTC  >  NC_000913/3804134‑3804313

gCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGt                                                                                         >  1:3539565/1‑93 (MQ=255)

gCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGt                                                                                         <  2:10587/93‑1 (MQ=255)

 cGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGtt                                                                                        >  1:1403708/1‑93 (MQ=255)

  gATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGttt                                                                                       <  2:860692/93‑1 (MQ=255)

  gATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGttt                                                                                       >  2:100663/1‑93 (MQ=255)

  gATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGttt                                                                                       >  2:2121783/1‑93 (MQ=255)

  gATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGttt                                                                                       >  2:3320387/1‑93 (MQ=255)

    tAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGa                                                                                     >  2:1794903/1‑93 (MQ=255)

    tAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGa                                                                                     >  1:1963521/1‑93 (MQ=255)

      aaTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATa                                                                                   <  2:1619364/93‑1 (MQ=255)

      aaTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATa                                                                                   <  2:3558680/93‑1 (MQ=255)

      aaTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATa                                                                                   >  1:4583459/1‑93 (MQ=255)

       aTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATaa                                                                                  <  2:4304466/93‑1 (MQ=255)

       aTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATaa                                                                                  >  2:3513970/1‑93 (MQ=255)

        tttCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAAt                                                                                 >  2:3626738/1‑93 (MQ=255)

            aCGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCtt                                                                             <  2:4105448/93‑1 (MQ=255)

            aCGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCtt                                                                             >  2:1447770/1‑93 (MQ=255)

            aCGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCtt                                                                             >  2:2394319/1‑93 (MQ=255)

            aCGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGCTCACGTTCGGCTTGATAATGCtt                                                                             >  1:2669621/1‑93 (MQ=255)

             cGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTc                                                                            >  2:1231427/1‑93 (MQ=255)

              ggCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCNNGTTCGGTTTGATAATGCTNCt                                                                           >  1:816429/1‑93 (MQ=255)

               gCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTg                                                                          >  2:802933/1‑93 (MQ=255)

               gCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTg                                                                          <  1:4299710/93‑1 (MQ=255)

                 tGTTTGGGATTTGCTCGCTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGa                                                                        >  2:2835613/1‑93 (MQ=255)

                 tGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGa                                                                        >  1:2415780/1‑93 (MQ=255)

                  gTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGaa                                                                       >  2:3200684/1‑93 (MQ=255)

                   tttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGaaa                                                                      >  1:2925698/1‑93 (MQ=255)

                   tttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGaaa                                                                      <  1:2025765/93‑1 (MQ=255)

                   tttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGaaa                                                                      >  1:4018962/1‑93 (MQ=255)

                    ttGGGATTTTCTCACTGTATTTTCTGTCCGAATAAATACCGGGAGAAAGAATTTGAAAACGTGCAGGTTCGGTTTGATAATGCTTCTGGAAAt                                                                     <  2:371764/93‑1 (MQ=255)

                    ttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCAGTTTGATAATGCTTCTGGAAAt                                                                     >  2:4408758/1‑93 (MQ=255)

                     tGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATc                                                                    >  1:4716300/1‑93 (MQ=255)

                          tttGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGa                                                               <  2:1033052/93‑1 (MQ=255)

                          tttGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGa                                                               <  2:1823646/93‑1 (MQ=255)

                           ttGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGAt                                                              <  1:4432923/93‑1 (MQ=255)

                              cTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTg                                                           >  2:1694254/1‑93 (MQ=255)

                               tCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGc                                                          >  2:1098501/1‑93 (MQ=255)

                                 aCTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCtt                                                        >  1:3106179/1‑93 (MQ=255)

                                 aCTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCtt                                                        >  2:4637942/1‑93 (MQ=255)

                                 aCTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCtt                                                        >  2:3164902/1‑93 (MQ=255)

                                   tGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTAt                                                      <  1:1132769/93‑1 (MQ=255)

                                    gTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATc                                                     <  1:3927457/93‑1 (MQ=255)

                                    gTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATc                                                     <  2:3342157/93‑1 (MQ=255)

                                    gTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGATTGATAATGCTTCTGGAAATCGGCGATTTGCTTATc                                                     >  1:858163/1‑93 (MQ=255)

                                       ttttCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGt                                                  >  2:535791/1‑93 (MQ=255)

                                       ttttCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGt                                                  <  2:4469121/93‑1 (MQ=255)

                                        tttCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTc                                                 <  1:2270832/93‑1 (MQ=255)

                                        tttCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTc                                                 <  1:1079693/93‑1 (MQ=255)

                                           cTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGc                                              <  2:3820910/93‑1 (MQ=255)

                                             gTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGcat                                            >  1:3319285/1‑93 (MQ=255)

                                             gTCCGGATAAATACCAGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGcat                                            <  2:4074661/93‑1 (MQ=255)

                                              tCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGcatc                                           <  1:9226/93‑1 (MQ=255)

                                              tCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGcatc                                           >  2:2190401/1‑93 (MQ=255)

                                                  gATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATaa                                       >  2:1296304/1‑93 (MQ=255)

                                                   aTAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAg                                      >  2:1984782/1‑93 (MQ=255)

                                                     aaaTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGtt                                    >  2:2022794/1‑93 (MQ=255)

                                                      aaTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGttt                                   >  2:4276663/1‑93 (MQ=255)

                                                      aaTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGttt                                   <  2:1154908/93‑1 (MQ=255)

                                                      aaTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGttt                                   <  1:3719440/93‑1 (MQ=255)

                                                       aTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATGGGTCAGCATCATAAGTTtc                                  <  1:1193770/93‑1 (MQ=255)

                                                       aTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtc                                  <  1:4183749/93‑1 (MQ=255)

                                                       aTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtc                                  >  1:1526394/1‑93 (MQ=255)

                                                            gggAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGa                             <  1:4368761/93‑1 (MQ=255)

                                                            gggAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGa                             <  1:287406/93‑1 (MQ=255)

                                                             ggaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAt                            >  2:761377/1‑93 (MQ=255)

                                                             ggaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAt                            >  2:4498955/1‑93 (MQ=255)

                                                              gaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAtt                           >  1:1486004/1‑93 (MQ=255)

                                                              gaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAtt                           <  1:1876459/93‑1 (MQ=255)

                                                              gaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAtt                           <  1:620516/93‑1 (MQ=255)

                                                               aggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAttt                          <  1:2535399/93‑1 (MQ=255)

                                                                 gaagaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAc                        >  1:126950/1‑93 (MQ=255)

                                                                  aagaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAcc                       <  2:710409/93‑1 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCNTNATAAGTTTCGTCGATTTAccc                      >  1:2428096/1‑93 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      >  2:913155/1‑93 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      >  2:3997481/1‑93 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      >  2:1498234/1‑93 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      <  1:3226871/93‑1 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      <  2:1177610/93‑1 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      >  2:2293101/1‑93 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      >  1:4031901/1‑93 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      >  1:346167/1‑93 (MQ=255)

                                                                    gaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCt                     >  2:3279051/1‑93 (MQ=255)

                                                                           aaaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGaa              <  2:4332173/93‑1 (MQ=255)

                                                                            aaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGaaa             >  2:2654469/1‑93 (MQ=255)

                                                                            aaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGaaa             >  2:4680068/1‑93 (MQ=255)

                                                                             aaCGTTCAGGTTCGGTTTGATAATGCTTCTGGTAATCGGCGATTTGCTTATCGGTCAGCATCATTAGTTTCGTTGATTTACCCTGCTCGAAAg            >  1:4222782/1‑93 (MQ=255)

                                                                             aaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAg            >  1:1028113/1‑93 (MQ=255)

                                                                              aCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGt           <  2:2207823/93‑1 (MQ=255)

                                                                               cGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtc          >  1:4500191/1‑93 (MQ=255)

                                                                                gttcaggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcg         >  2:590333/1‑93 (MQ=255)

                                                                                 ttcaggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTGAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgc        >  1:3065493/1‑93 (MQ=255)

                                                                                 ttcaggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgc        >  2:2345701/1‑93 (MQ=255)

                                                                                 ttcaggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgc        >  1:2592056/1‑93 (MQ=255)

                                                                                  tcaggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgct       <  2:622073/93‑1 (MQ=255)

                                                                                  tcaggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgct       <  2:4154183/93‑1 (MQ=255)

                                                                                   caggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctc      >  1:176227/1‑93 (MQ=255)

                                                                                   caggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctc      <  1:4726195/93‑1 (MQ=255)

                                                                                    aggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctcg     >  1:191724/1‑93 (MQ=255)

                                                                                    aggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctcg     <  1:2323202/93‑1 (MQ=255)

                                                                                       ttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTCgctcgctc  >  1:1632009/1‑93 (MQ=255)

GCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCAGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTCGCTCGCTC  >  NC_000913/3804134‑3804313

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