BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,542,191	A→C	49.2%	D52A (GAC→GCC)	fimE →	regulator for fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,191	0	A	C	49.2%	27.7 / 282.1	124	D52A (GAC→GCC)	fimE	regulator for fimA
Reads supporting (aligned to +/- strand): ref base A (32/31); new base C (23/38); total (55/69)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

TTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG  >  NC_000913/4542093‑4542286

ttGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGga                                                                                                 <  1:1401410/99‑1 (MQ=255)

   ttACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgacct                                                                                              <  1:891362/99‑1 (MQ=255)

   ttACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgacct                                                                                              <  2:2003004/99‑1 (MQ=255)

   ttACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgacct                                                                                              <  2:1453/99‑1 (MQ=255)

        gggCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgcccttgacc                                                                                         <  1:2678956/99‑1 (MQ=255)

        gggCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgcccttgacc                                                                                         >  1:994380/1‑99 (MQ=255)

            aaCGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTaa                                                                                     <  2:973408/99‑1 (MQ=255)

              cGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATg                                                                                   <  1:2003051/99‑1 (MQ=255)

               gggAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGa                                                                                  <  1:1881566/99‑1 (MQ=255)

                ggAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGaa                                                                                 <  2:1328580/99‑1 (MQ=255)

                ggAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGaa                                                                                 <  1:2734854/99‑1 (MQ=255)

                 ggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAg                                                                                <  1:938572‑M2/48‑1 (MQ=255)

                  aGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAgg                                                                               <  1:1004802/99‑1 (MQ=255)

                    ccAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTa                                                                             >  2:1024123/1‑99 (MQ=255)

                      agagATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGa                                                                           >  1:1539612/1‑99 (MQ=255)

                      agagATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGa                                                                           <  1:1138748/99‑1 (MQ=255)

                       gagaTTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGaa                                                                          <  2:487959/99‑1 (MQ=255)

                       gagaTTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaa                                                                          <  1:2702029/99‑1 (MQ=255)

                       gagaTTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaa                                                                          >  1:460315/1‑99 (MQ=255)

                        agaTTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGaat                                                                         <  1:822987/99‑1 (MQ=255)

                        agaTTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaat                                                                         <  2:856301/99‑1 (MQ=255)

                           ttattGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGaataaa                                                                      <  1:460175/99‑1 (MQ=255)

                           ttattGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGaataaa                                                                      <  2:954986/99‑1 (MQ=255)

                             attGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaataaata                                                                    <  2:2359183/99‑1 (MQ=255)

                              ttGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAatat                                                                   <  2:2338176/99‑1 (MQ=255)

                               tGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATAtt                                                                  >  2:2158988/1‑99 (MQ=255)

                                gTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTc                                                                 <  1:2636702/99‑1 (MQ=255)

                                 tCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCg                                                                <  2:2175409/99‑1 (MQ=255)

                                  cTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGc                                                               <  2:22675/99‑1 (MQ=255)

                                  cTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGc                                                               <  2:737852/99‑1 (MQ=255)

                                    tATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCg                                                             >  2:1445325/1‑99 (MQ=255)

                                     atctcatatcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGa                                                            >  2:2822162‑M2/32‑99 (MQ=255)

                                      ttCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGAc                                                           <  2:2667057/99‑1 (MQ=255)

                                       tCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACt                                                          <  1:902554/99‑1 (MQ=255)

                                       tCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACt                                                          <  1:193857/99‑1 (MQ=255)

                                         tGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTga                                                        >  1:1189169/1‑99 (MQ=255)

                                         tGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTga                                                        >  1:670431/1‑99 (MQ=255)

                                             tcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTgaagaa                                                    <  1:96295‑M2/76‑1 (MQ=255)

                                                atatCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTACAATAAATATTCGCCGACTGAAGAACgg                                                 >  2:2740487/1‑99 (MQ=255)

                                                 tatCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGa                                                >  1:547910/1‑99 (MQ=255)

                                                 tatCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGa                                                >  2:623618/1‑99 (MQ=255)

                                                   tCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGAtt                                              <  1:689389/99‑1 (MQ=255)

                                                   tCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGAtt                                              >  1:1115051/1‑99 (MQ=255)

                                                    cttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGAttt                                             >  1:1971635‑M2/17‑99 (MQ=255)

                                                    cGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGAttt                                             <  2:2831987/99‑1 (MQ=255)

                                                     ggCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGAtttt                                            >  1:2806523/1‑99 (MQ=255)

                                                      gCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTc                                           <  1:2106793/99‑1 (MQ=255)

                                                       cATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCt                                          >  1:1989261/1‑99 (MQ=255)

                                                        aTGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTa                                         <  1:2638885/99‑1 (MQ=255)

                                                        aTGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTa                                         <  2:2021474/99‑1 (MQ=255)

                                                         tGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTAc                                        >  2:2013954/1‑99 (MQ=255)

                                                         tGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTAc                                        >  1:1092959/1‑99 (MQ=255)

                                                         tGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTAc                                        <  1:2339965/99‑1 (MQ=255)

                                                         tGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTAc                                        >  2:1331574/1‑99 (MQ=255)

                                                          gggATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTAcc                                       <  1:1286352/99‑1 (MQ=255)

                                                           ggATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCg                                      >  2:698322/1‑99 (MQ=255)

                                                           ggATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCg                                      >  1:727386/1‑99 (MQ=255)

                                                           ggATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCg                                      <  1:764061/99‑1 (MQ=255)

                                                           ggATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCg                                      <  1:689042/99‑1 (MQ=255)

                                                            gATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGt                                     >  2:2179445/1‑99 (MQ=255)

                                                            gATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGt                                     <  2:488140/99‑1 (MQ=255)

                                                             aTGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGtt                                    >  1:134235/1‑99 (MQ=255)

                                                             aTGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGtt                                    >  1:2351867/1‑99 (MQ=255)

                                                              tGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTc                                   >  1:2250608/1‑99 (MQ=255)

                                                              tGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTc                                   >  1:816596/1‑99 (MQ=255)

                                                                cGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCAc                                 <  1:2197811/99‑1 (MQ=255)

                                                                 gTATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCAcc                                <  2:1054028/99‑1 (MQ=255)

                                                                  tATTAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCAccc                               <  2:465948/99‑1 (MQ=255)

                                                                   aTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCg                              <  1:2678969/99‑1 (MQ=255)

                                                                    ttAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGt                             <  2:1588612/99‑1 (MQ=255)

                                                                    ttAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGt                             >  2:1174920/1‑99 (MQ=255)

                                                                    ttAGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGt                             <  1:981577/99‑1 (MQ=255)

                                                                    ttAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGt                             >  2:1469024/1‑99 (MQ=255)

                                                                    ttAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGt                             <  2:322088/99‑1 (MQ=255)

                                                                      aGTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTa                           <  2:1814494/99‑1 (MQ=255)

                                                                       gTGAACTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTAc                          >  2:862296/1‑99 (MQ=255)

                                                                          aaCTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGtt                       >  2:177050/1‑99 (MQ=255)

                                                                           aCTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGttt                      >  1:1271950/1‑99 (MQ=255)

                                                                           aCTGCTTGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGttt                      <  1:2122144/99‑1 (MQ=255)

                                                                            cTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGtttt                     >  1:2500637/1‑99 (MQ=255)

                                                                             tGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTtg                    >  2:2357395/1‑99 (MQ=255)

                                                                             tGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTtg                    >  1:1506211/1‑99 (MQ=255)

                                                                                ttGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTtgatg                 >  1:157079/1‑99 (MQ=255)

                                                                                 tGATCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTtgatga                <  1:1473411/99‑1 (MQ=255)

                                                                                 tGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTtgatga                <  1:1912086/99‑1 (MQ=255)

                                                                                   aTCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGc              >  1:2838318/1‑99 (MQ=255)

                                                                                    tCTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCg             <  1:2236137/99‑1 (MQ=255)

                                                                                    tCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCg             <  1:586339/99‑1 (MQ=255)

                                                                                     cTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGt            <  1:56355/99‑1 (MQ=255)

                                                                                     cTGCATTATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGt            >  1:742908/1‑99 (MQ=255)

                                                                                     cTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGt            <  1:119987/99‑1 (MQ=255)

                                                                                      tGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTg           >  1:2318048/1‑99 (MQ=255)

                                                                                          ttATCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGc       >  2:576147/1‑99 (MQ=255)

                                                                                          ttATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGc       >  1:2762916/1‑99 (MQ=255)

                                                                                           tATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGcc      <  1:2130823/99‑1 (MQ=255)

                                                                                           tATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGcc      >  2:1227945/1‑99 (MQ=255)

                                                                                            aTCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCg     >  2:2635179/1‑99 (MQ=255)

                                                                                            aTCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCg     <  1:1433844/99‑1 (MQ=255)

                                                                                             tCAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGt    >  2:2149625/1‑99 (MQ=255)

                                                                                              cAGGCCCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTg   <  2:1671873/99‑1 (MQ=255)

TTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG  >  NC_000913/4542093‑4542286

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