BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	805,313	A→C	16.4%	E604A (GAG→GCG) ‡	ybhJ →	aconitase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	805,313	0	A	C	16.4%	447.3 / 89.3	190	E604A (GAG→GCG)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand): ref base A (78/80); new base C (12/19); total (90/100)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-02

Only 100 of 192 total aligned reads displayed.

GGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACC  >  NC_000913/805215‑805412

ggCGGAGTTTGCCCGGTCGCGCCGCGATCCCGATTATGTTAGCAGAAGTAAAGCGCCTGCTCAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAg                                                                                                    <  1:1124385/100‑1 (MQ=255)

ggCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAg                                                                                                    >  2:2692565/1‑100 (MQ=255)

ggCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTATGGAGGGGAATGTCAGCGAg                                                                                                    >  1:4096491/1‑100 (MQ=255)

 gCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGc                                                                                                   <  1:1305326/100‑1 (MQ=255)

  cGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCATCGTCTGGCGGGGAATGTCAGCGAGCt                                                                                                  <  1:4300559/100‑1 (MQ=255)

     aGTTTACCCTGTCTCGGCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTCGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGAc                                                                                               <  1:2727272/100‑1 (MQ=255)

     aGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGAc                                                                                               >  1:3426424/1‑100 (MQ=255)

      gTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACa                                                                                              >  2:1876380/1‑100 (MQ=255)

        ttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACatt                                                                                            <  1:660066/100‑3 (MQ=255)

         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattc                                                                                           <  2:2245798/100‑4 (MQ=255)

         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattc                                                                                           <  2:53064/100‑4 (MQ=255)

         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACagag                                                                                           >  1:2049398/1‑100 (MQ=255)

         tACCCTGTCTCGCCGCGATCCCGGTTATGCTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattc                                                                                           <  2:4335296/100‑4 (MQ=255)

          aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattcc                                                                                          <  2:1066998/100‑5 (MQ=255)

          aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCCTCTGGCGGGGAATGTCAGCGAGCTGACAGAgg                                                                                          <  2:4119554/100‑1 (MQ=255)

           cccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattccc                                                                                         <  2:1311045/100‑6 (MQ=255)

           cccTGTCGCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAgg                                                                                          <  1:2911783/99‑1 (MQ=255)

            ccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattcccc                                                                                        <  2:157300/100‑7 (MQ=255)

            ccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattcccc                                                                                        <  1:2896353/100‑7 (MQ=255)

            ccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgtg                                                                                        <  1:1546257/100‑1 (MQ=255)

              tGTCTCGCTGCGGTCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAACATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattccccg                                                                                       <  1:2197372/99‑8 (MQ=255)

                tctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTg                                                                                    <  1:4262368/100‑1 (MQ=255)

                 ctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgc                                                                                   >  2:3666581/1‑100 (MQ=255)

                 ctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgc                                                                                   >  2:3025478/1‑100 (MQ=255)

                   cgccgcGATCCTGGTTATGTTAGCAGAAGTAAAGCGACTGCTCAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTacgc                                                                                 >  1:189090/1‑100 (MQ=255)

                   cgccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcgc                                                                                 >  1:4187528/1‑100 (MQ=255)

                   cgccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTg                                                                                    <  1:3655562/97‑1 (MQ=255)

                    gacgcgATCCCGGCTATGTTAGCTGAAGTAAAGCGCCTGCTGAGCTGGACAAACAGCCTCTGGCGGGGAATGTCAGGGAGCTGACAGAGGCGTTTgcgc                                                                                 >  1:3498868/3‑99 (MQ=255)

                     ccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcgcgc                                                                               >  1:3537639/1‑100 (MQ=255)

                      cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                                                                              <  2:293722/100‑1 (MQ=255)

                       gcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAt                                                                             <  2:1375954/100‑1 (MQ=255)

                        cgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                            >  2:1835759/1‑100 (MQ=255)

                         gATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTa                                                                           >  1:938531/1‑100 (MQ=255)

                          aTCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTaa                                                                          >  1:1315216/1‑100 (MQ=255)

                            cccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAg                                                                         <  1:1269803/99‑1 (MQ=255)

                            cccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGc                                                                        <  2:3930353/100‑1 (MQ=255)

                            cccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGc                                                                        <  1:802885/100‑1 (MQ=255)

                             ccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCa                                                                       <  1:2637757/100‑1 (MQ=255)

                             ccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCa                                                                       <  2:2415582/100‑1 (MQ=255)

                               ggTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGa                                                                     >  2:1131764/1‑100 (MQ=255)

                               ggTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGa                                                                     >  1:1147397/1‑100 (MQ=255)

                                    tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCt                                                                <  1:3760210/100‑1 (MQ=255)

                                    tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCt                                                                >  1:2378108/1‑100 (MQ=255)

                                    tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCt                                                                <  2:2049405/100‑1 (MQ=255)

                                     gTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCt                                                                >  2:24487/1‑99 (MQ=255)

                                     gTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTg                                                               <  2:743905/100‑1 (MQ=255)

                                      ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCTGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTgg                                                              <  2:3095455/100‑1 (MQ=255)

                                       tAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTAGCGGGGAATGTCAGCGAGCTGACAGAGATGTTTGCGCGCATTAAGCAGATTGCTGGt                                                             <  1:470114/100‑1 (MQ=255)

                                        ttccgatctAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCACATTAAGCAGATTGCTGGTc                                                            <  1:3102060/92‑1 (MQ=255)

                                         gCAGAAGTAAAGCGACTGATGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGCTGGCGGGCATTAAGCAGATTg                                                                  >  1:3670465/1‑93 (MQ=255)

                                          cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAg                                                          <  2:133437/100‑1 (MQ=255)

                                           aGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAgg                                                         <  2:3158531/100‑1 (MQ=255)

                                            gAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGa                                                        >  1:391483/1‑100 (MQ=255)

                                             aaGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                       <  2:3713754/100‑1 (MQ=255)

                                               gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa                                                     <  2:376716/100‑1 (MQ=255)

                                               gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa                                                     <  1:2782787/100‑1 (MQ=255)

                                                tgaaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                    >  2:2964165/3‑100 (MQ=255)

                                                   aGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCatat                                                  >  1:2059956/1‑99 (MQ=255)

                                                   aGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATAtt                                                 <  2:1147401/100‑1 (MQ=255)

                                                        ctgctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                            >  2:260013/1‑100 (MQ=255)

                                                        ctgctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                            <  1:3244876/100‑1 (MQ=255)

                                                         tgctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCg                                           <  2:2122033/100‑1 (MQ=255)

                                                          gctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGc                                          <  1:3277247/100‑1 (MQ=255)

                                                           ctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCt                                         >  2:4061844/1‑100 (MQ=255)

                                                            tgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTg                                        <  2:2383379/100‑1 (MQ=255)

                                                               gCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCaa                                     >  2:2082221/1‑100 (MQ=255)

                                                                cTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCaaa                                    >  2:3645619/1‑100 (MQ=255)

                                                                 tGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                   >  2:2781816/1‑100 (MQ=255)

                                                                 tGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                   <  1:2487298/100‑1 (MQ=255)

                                                                 tGGAAAACCAGCGTCTGGCGGGGACTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                   <  2:3081709/100‑1 (MQ=255)

                                                                  ggAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACt                                  >  2:1210309/1‑100 (MQ=255)

                                                                  ggAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCACATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACt                                  >  2:357113/1‑100 (MQ=255)

                                                                    aaaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGa                                >  1:2670691/1‑100 (MQ=255)

                                                                    aaaaTCAGCGCCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGTGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGa                                >  1:2388883/1‑100 (MQ=255)

                                                                        tCAGCGTCTGGCGGGGAATGTCATCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATCGATCCGCTGCCAACAGACAtt                            >  1:4293061/1‑100 (MQ=255)

                                                                        tCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                            >  2:3255540/1‑100 (MQ=255)

                                                                        tCAGCGTCTGGCGGGGAATGTCAGCGAGCGGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCCGGAGCATATTGATCCGCTGCAAACTGAAAtt                            >  1:718629/1‑100 (MQ=255)

                                                                         cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTg                           <  2:2684810/100‑1 (MQ=255)

                                                                          aGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTgg                          >  2:2158941/1‑100 (MQ=255)

                                                                           gCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgc                         >  1:1990434/1‑100 (MQ=255)

                                                                               cTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagca                     >  2:940905/1‑100 (MQ=255)

                                                                                caaacacctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                    >  2:1797271/11‑100 (MQ=38)

                                                                                caaacacctctGTCAGCTCGCTGACAGAGGTGGTTGCGCAACTTGAGCAGCTTGCTGGTCCGGAGAATCTTGATCCGCTGCAACCTGCAGTTGGAAGCAt                    >  1:3891429/11‑100 (MQ=12)

                                                                                 ggCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATg                   >  1:959384/1‑100 (MQ=255)

                                                                                   cggaggaTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATgg                  <  1:4305390/93‑1 (MQ=255)

                                                                                    cacctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                >  2:1519028/7‑100 (MQ=255)

                                                                                    cacctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                >  1:773422/7‑100 (MQ=255)

                                                                                      cctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa              >  2:3525984/5‑100 (MQ=255)

                                                                                      cctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa              >  2:2528136/5‑100 (MQ=255)

                                                                                        tctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGATGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATg            >  1:2450044/3‑100 (MQ=255)

                                                                                        aaTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTAt             <  1:3305470/99‑1 (MQ=255)

                                                                                         aTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc           >  2:2341034/1‑100 (MQ=255)

                                                                                         aTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc           <  1:3824300/100‑1 (MQ=255)

                                                                                          tGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCg          <  1:260012/100‑1 (MQ=255)

                                                                                             cAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg       <  1:4061814/100‑1 (MQ=255)

                                                                                             cAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg       <  2:4259921/100‑1 (MQ=255)

                                                                                             cAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg       <  1:962391/100‑1 (MQ=255)

                                                                                              aGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGa      <  1:3378765/100‑1 (MQ=255)

                                                                                                  aGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAcc  <  2:680410/100‑1 (MQ=255)

                                                                                                  aGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAcc  >  1:1624848/1‑100 (MQ=255)

GGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACC  >  NC_000913/805215‑805412

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