BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,375	+A	coding (458/516 nt)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,372	1	.	A	99.0%	346.5 / ‑3.8	96	coding (455/516 nt)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base A (32/63); total (33/63)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑AAAGGCCATGCGCTGGAA‑G‑GCTGGCG‑‑‑‑AA‑‑‑‑‑CT‑GGCTACAG‑CTGAAGGATTG‑ATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCC  >  NC_000913/4008280‑4008449

gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTc                                                                                                  <  2:341848/93‑1 (MQ=255)

gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTc                                                                                                  <  2:4402295/93‑1 (MQ=255)

 cTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑a                                                                                                >  2:2326490/1‑93 (MQ=255)

 cTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑a                                                                                                >  1:5244410/1‑93 (MQ=255)

  tGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑aa                                                                                               >  2:2027251/1‑93 (MQ=255)

  tGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑aa                                                                                               <  1:472365/93‑1 (MQ=255)

  tGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑aa                                                                                               <  1:3256673/93‑1 (MQ=255)

   agAGCAGGCGATTAACGCTCGTTGGTGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑aaa                                                                                              <  2:1074239/92‑1 (MQ=255)

    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑aaaa                                                                                             <  1:2285683/93‑2 (MQ=255)

    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑aaaa                                                                                             <  2:4163007/93‑2 (MQ=255)

    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GGTTC‑aaaa                                                                                             >  2:2686103/1‑92 (MQ=255)

    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GGTTC‑aaaa                                                                                             >  2:4295461/1‑92 (MQ=255)

     aGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTT‑aaaag                                                                                            >  1:3192920/1‑91 (MQ=255)

       cAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑aaaaggc                                                                                          >  2:1634344/1‑89 (MQ=255)

       cAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑aaaaggc                                                                                          <  2:5700231/93‑5 (MQ=255)

         ggCGATTAACGCTCGTTGGGGCTATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCa                                                                                         <  1:6112872/93‑1 (MQ=255)

         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCa                                                                                         <  1:1506873/93‑1 (MQ=255)

         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCa                                                                                         <  2:4951341/93‑1 (MQ=255)

         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCa                                                                                         <  1:5722761/93‑1 (MQ=255)

          gCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCAt                                                                                        <  2:1599366/93‑1 (MQ=255)

          gCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCAt                                                                                        >  2:4819971/1‑93 (MQ=255)

           cGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATg                                                                                       <  2:2169652/93‑1 (MQ=255)

           cGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATg                                                                                       <  1:4269990/93‑1 (MQ=255)

            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATgc                                                                                      <  1:2157699/93‑1 (MQ=255)

            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATgc                                                                                      <  2:5492824/93‑1 (MQ=255)

            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATgc                                                                                      >  1:5172308/1‑93 (MQ=255)

            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATgc                                                                                      <  1:3756419/93‑1 (MQ=255)

            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATgc                                                                                      <  1:2478711/93‑1 (MQ=255)

             aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCAcgcg                                                                                     >  2:1385819/1‑93 (MQ=255)

             aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATgcg                                                                                     >  1:1442960/1‑93 (MQ=255)

              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATgcgc                                                                                    >  1:5228684/1‑93 (MQ=255)

              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAAGCCATgcgc                                                                                    >  1:5071661/1‑93 (MQ=255)

              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAAGCCATgcgc                                                                                    >  2:5994215/1‑93 (MQ=255)

                aaCGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTg                                                                                  <  2:2092422/93‑1 (MQ=255)

                 aCGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTgg                                                                                 <  2:2390309/93‑1 (MQ=255)

                  cGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTTC‑AAAAGACATGCGCTTGa                                                                                >  2:557163/1‑93 (MQ=255)

                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGaa                                                                               <  2:1607153/93‑1 (MQ=255)

                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGaa                                                                               >  1:2885126/1‑93 (MQ=255)

                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGaa                                                                               <  1:2474865/93‑1 (MQ=255)

                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGaa                                                                               <  1:868213/93‑1 (MQ=255)

                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGaa                                                                               <  1:926211/93‑1 (MQ=255)

                    cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑g                                                                             <  1:5707178/93‑1 (MQ=255)

                      cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑gcg                                                                           <  2:4574996/93‑3 (MQ=255)

                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGG‑GGTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                          >  1:4813315/1‑90 (MQ=255)

                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                          <  1:4623304/93‑4 (MQ=255)

                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                          <  2:3083575/93‑4 (MQ=255)

                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                          <  2:3097571/93‑4 (MQ=255)

                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                          <  2:3282794/93‑4 (MQ=255)

                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GGTTCAAAAAGCCATGCGCTGGAA‑gcgg                                                                          >  1:4041485/1‑90 (MQ=255)

                        ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑gcggt                                                                         <  2:3738063/93‑5 (MQ=255)

                             gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑a                                                                >  1:882407/1‑93 (MQ=255)

                             gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑a                                                                <  2:722109/93‑1 (MQ=255)

                             gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑a                                                                >  2:2003960/1‑93 (MQ=255)

                             gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GGTTCAAAAGGGCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑a                                                                >  2:4689360/1‑93 (MQ=255)

                              cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aa                                                               <  2:2103880/93‑1 (MQ=255)

                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                             <  2:5143571/93‑3 (MQ=255)

                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                             <  2:219070/93‑3 (MQ=255)

                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                             <  2:4328059/93‑3 (MQ=255)

                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                             <  2:3925526/93‑3 (MQ=255)

                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                             >  1:4242722/1‑91 (MQ=255)

                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                             >  1:3268770/1‑91 (MQ=255)

                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCACTGGCTGGAAGG‑GGTGGCG‑‑‑‑aaaa                                                             >  1:1410130/1‑91 (MQ=255)

                                 tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aagaa                                                            >  1:2474435/1‑90 (MQ=255)

                                  cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aagaaa                                                           <  1:3617388/93‑5 (MQ=255)

                                  cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aagaaa                                                           <  2:2759649/93‑5 (MQ=255)

                                                                  tCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                        <  1:3067422/93‑1 (MQ=255)

                                                                  tCTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                        >  2:5367497/1‑93 (MQ=255)

                                                                   cTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                       <  1:3423220/93‑2 (MQ=255)

                                                                   cTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                       >  2:2690426/1‑92 (MQ=255)

                                                                   cTGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                       >  2:2607369/1‑92 (MQ=255)

                                                                    tGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      >  2:4708075/1‑91 (MQ=255)

                                                                    tGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      <  1:1051706/93‑3 (MQ=255)

                                                                    tGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      >  1:5951006/1‑91 (MQ=255)

                                                                    tGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      >  1:5995372/1‑91 (MQ=255)

                                                                    tGGAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      >  2:1156237/1‑91 (MQ=255)

                                                                     ggAAGTGATGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttgg                     >  2:630739/1‑90 (MQ=255)

                                                                             tGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa             <  1:100493/93‑1 (MQ=255)

                                                                             tGGCGGGCGGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa             <  2:349093/93‑1 (MQ=255)

                                                                             tGGCGGGCGGC‑GTTTCAAAAAGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa             >  2:6045798/1‑93 (MQ=255)

                                                                              ggcgggcgGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt            <  2:2148939/93‑1 (MQ=255)

                                                                              ggcgggcgGC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑AAGAAACCG‑GGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt            <  2:1686766/93‑1 (MQ=255)

                                                                               gcgggcggC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg           <  2:2236559/93‑1 (MQ=255)

                                                                               gcgggcggC‑GTGTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg           <  2:790947/93‑1 (MQ=255)

                                                                                cgggcggC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGa          <  1:5478423/93‑1 (MQ=255)

                                                                                 gggcggC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa         <  1:2635131/93‑1 (MQ=255)

                                                                                 gggcggC‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa         <  2:3445771/93‑1 (MQ=255)

                                                                                  ggcggc‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac        <  2:2772791/93‑1 (MQ=255)

                                                                                  ggcggc‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac        <  1:4957056/93‑1 (MQ=255)

                                                                                  ggcggc‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac        >  1:3993000/1‑93 (MQ=255)

                                                                                  ggcggc‑GTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac        <  2:331868/93‑1 (MQ=255)

                                                                                   gcggc‑gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg       >  1:4814006/1‑93 (MQ=255)

                                                                                   gcggc‑gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg       <  1:2957114/93‑1 (MQ=255)

                                                                                     ggc‑gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTT‑‑GTGACGGGATGAacgac     >  2:5773184/1‑93 (MQ=255)

                                                                                     ggc‑gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac     <  2:4498313/93‑1 (MQ=255)

                                                                                     ggc‑gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac     >  2:5281397/1‑93 (MQ=255)

                                                                                     ggc‑gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac     <  2:4279944/93‑1 (MQ=255)

                                                                                     ggc‑gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac     >  2:3119511/1‑93 (MQ=255)

                                                                                      gc‑gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg    >  2:4300379/1‑93 (MQ=255)

                                                                                       c‑gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGc   >  2:5938791/1‑93 (MQ=255)

                                                                                         gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGcc  <  1:4418546/93‑1 (MQ=255)

GCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTC‑AAAGGCCATGCGCTGGAA‑G‑GCTGGCG‑‑‑‑AA‑‑‑‑‑CT‑GGCTACAG‑CTGAAGGATTG‑ATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCC  >  NC_000913/4008280‑4008449

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