BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,991,301	C→T	80.3%	P50L (CCG→CTG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,301	0	C	T	80.3%	263.7 / 62.7	132	P50L (CCG→CTG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base C (15/11); new base T (55/51); total (70/62)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

Only 100 of 134 total aligned reads displayed.

AATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAAT  >  NC_000913/3991202‑3991400

aaTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGcc                                                                                                     <  2:1618138/100‑1 (MQ=255)

   tGCGTGTGGATCTCGGGGTTGCCGGTATCGGTACAGCATTCAAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTggg                                                                                                  >  1:627187/1‑100 (MQ=255)

   tGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTggg                                                                                                  <  2:876795/100‑1 (MQ=255)

    gCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGt                                                                                                 <  1:623824/100‑1 (MQ=255)

     cGTGTGGATCGCGAGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGtt                                                                                                >  2:1955402/1‑100 (MQ=255)

      gtgtGGATCGCGGGCTTGCCGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTa                                                                                               >  1:1947826/1‑100 (MQ=255)

      gtgtGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTa                                                                                               >  1:1947976/1‑100 (MQ=255)

      gtgtGGATCGCGCGCGTCCCCCTATGGTGCCAGCATACCAACAGGTCTACAGGCTACTGCCGACATTGTTACACTATCACCATCCGCTAATGCTAGGTTa                                                                                               >  1:140765/1‑100 (MQ=255)

        gtgGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTAcc                                                                                             >  1:2527189/1‑100 (MQ=255)

          ggATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCtt                                                                                           <  1:21922/100‑1 (MQ=255)

           gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTg                                                                                          <  1:1406972/100‑1 (MQ=255)

           gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTg                                                                                          <  1:2228910/100‑1 (MQ=255)

            aTCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGa                                                                                         >  2:192688/1‑100 (MQ=255)

            aTCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGa                                                                                         >  2:17796/1‑100 (MQ=255)

            aTCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGa                                                                                         >  1:1537650/1‑100 (MQ=255)

              cgcgcgCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATg                                                                                       <  2:1576376/100‑1 (MQ=255)

              cgcgcgCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATg                                                                                       <  1:17795/100‑1 (MQ=255)

               gcgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATg                                                                                       >  2:495352/1‑99 (MQ=255)

                cgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGt                                                                                     >  2:787202/1‑100 (MQ=255)

                   gcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAAc                                                                                  >  2:1443385/1‑100 (MQ=255)

                    cTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAAc                                                                                  <  2:104169/99‑1 (MQ=255)

                    cTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACg                                                                                 >  2:865381/1‑100 (MQ=255)

                    cTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACg                                                                                 >  2:865418/1‑100 (MQ=255)

                      tgctgcTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGtt                                                                               >  2:2192382/1‑100 (MQ=255)

                       gctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTc                                                                              <  2:1330917/100‑1 (MQ=255)

                       gctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTc                                                                              >  1:890957/1‑100 (MQ=255)

                        ctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTcc                                                                             >  2:455142/1‑100 (MQ=255)

                        ctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTcc                                                                             >  2:701037/1‑100 (MQ=255)

                        ctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTcc                                                                             <  1:1132557/100‑1 (MQ=255)

                          ggtATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCa                                                                           >  2:1957103/3‑100 (MQ=255)

                          gctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCa                                                                           >  2:325806/1‑100 (MQ=255)

                          gctACGGGGTCGGCATTCCAAGAGGTGTAGATTCTACTGCCGACATTGTTGCACTATCACCAGCCGCTAATGCTGGGTTACCTTGATGGTAACGTGCCCa                                                                           <  1:2254568/100‑1 (MQ=255)

                              tGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAgg                                                                       <  2:2333195/100‑1 (MQ=255)

                              tGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAgg                                                                       <  2:426307/100‑1 (MQ=255)

                              tGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAgg                                                                       <  2:2451801/100‑1 (MQ=255)

                                 ggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAt                                                                    >  2:114148/1‑100 (MQ=255)

                                   ccTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAttt                                                                  >  1:1233927/1‑100 (MQ=255)

                                     tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGc                                                                >  2:628364/1‑100 (MQ=255)

                                      gCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTCCCCAAAGGCATTTGcc                                                               <  2:888287/100‑1 (MQ=255)

                                       cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCt                                                              >  2:2257703/1‑100 (MQ=255)

                                         ttCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                            <  2:18380/100‑1 (MQ=255)

                                           ccAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTac                                                          >  1:1865264/1‑100 (MQ=255)

                                            cAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTaca                                                         >  2:1117962/1‑100 (MQ=255)

                                            cAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttggg                                                         <  1:2030130/100‑4 (MQ=255)

                                             aaCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAACGCATTTGCCTTTacac                                                        >  1:1000209/1‑100 (MQ=255)

                                              aCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttgggaa                                                       <  2:1676913/100‑6 (MQ=255)

                                              aCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACg                                                       <  1:1373852/100‑1 (MQ=255)

                                                aGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGcc                                                     <  2:2665884/100‑1 (MQ=255)

                                                aGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttgggaacg                                                     <  1:1373479/100‑8 (MQ=255)

                                                 ggTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCt                                                    >  1:2784293/1‑100 (MQ=255)

                                                  gTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                   >  1:265581/1‑100 (MQ=255)

                                                   tCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                  <  1:147381/100‑1 (MQ=255)

                                                    cTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgat                                                 >  1:372129/1‑100 (MQ=255)

                                                      aCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgatga                                               <  2:2179918/100‑1 (MQ=255)

                                                      aCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgatga                                               <  1:562318/100‑1 (MQ=255)

                                                         ntCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTATCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACATGCCTGATGAAAc                                            >  1:2582269/2‑100 (MQ=255)

                                                         gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                            <  1:1096949/100‑1 (MQ=255)

                                                         gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                            <  2:13146/100‑1 (MQ=255)

                                                         gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                            >  1:1147704/1‑100 (MQ=255)

                                                         gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                            <  2:814237/100‑1 (MQ=255)

                                                         gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                            >  2:387619/1‑100 (MQ=255)

                                                              cTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                       >  2:2865236/1‑100 (MQ=255)

                                                                gCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGc                                     >  1:77277/1‑100 (MQ=255)

                                                                gCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGc                                     >  2:1965688/1‑100 (MQ=255)

                                                                 ccGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGcc                                    <  1:2192374/100‑1 (MQ=255)

                                                                 ccGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGcc                                    >  1:440921/1‑100 (MQ=255)

                                                                  cGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCa                                   <  2:941055/100‑1 (MQ=255)

                                                                   gACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCAc                                  <  2:626355/100‑1 (MQ=255)

                                                                    aCATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACt                                 >  1:541781/1‑100 (MQ=255)

                                                                     cATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTa                                >  2:12994/1‑100 (MQ=255)

                                                                     cATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTa                                >  1:244608/1‑100 (MQ=255)

                                                                       ttgttgCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTTATGAAACTCAACGCCACTAcc                              <  1:671906/100‑1 (MQ=255)

                                                                       ttgttgCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTAcc                              <  2:253970/100‑1 (MQ=255)

                                                                       ttgttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTAcc                              >  1:2778567/1‑100 (MQ=255)

                                                                        tgttgCACGATCACGATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTGACACGCGTGATGAGACTCAACGCCACTAcc                              <  1:565940/99‑1 (MQ=255)

                                                                         gttgCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTg                            <  2:1617318/100‑1 (MQ=255)

                                                                            gCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAAc                         <  2:1184438/100‑1 (MQ=255)

                                                                            gCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAAc                         >  1:2364872/1‑100 (MQ=255)

                                                                               cTATCACCATCCGCTAATCCTGGGTTACCTTCATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAAc                         >  1:809664/1‑97 (MQ=255)

                                                                                tATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                     <  1:1366247/100‑1 (MQ=255)

                                                                                tATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                     >  1:1097917/1‑100 (MQ=255)

                                                                                 aTCACCATCCGCTAATGCTTGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt                    <  1:693810/100‑1 (MQ=255)

                                                                                  tCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                     <  1:2739303/98‑1 (MQ=255)

                                                                                   cACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg                  >  2:764676/1‑100 (MQ=255)

                                                                                   cACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg                  >  2:250302/1‑100 (MQ=255)

                                                                                   cACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg                  >  2:1754092/1‑100 (MQ=255)

                                                                                    aCCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa                 >  1:1526716/1‑100 (MQ=255)

                                                                                    aCCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa                 >  1:753982/1‑100 (MQ=255)

                                                                                    aCCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACAAGCCTGATGAAACTCAACGCCACTAACTGAACGAGCTTGa                 >  1:2251634/1‑100 (MQ=255)

                                                                                     ccATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                >  1:2648670/1‑100 (MQ=255)

                                                                                      cATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc               <  2:1867415/100‑1 (MQ=255)

                                                                                       aTCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACt              >  2:2589972/1‑100 (MQ=255)

                                                                                        tCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTg             <  2:1733329/100‑1 (MQ=255)

                                                                                          cGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTa           >  2:1026199/1‑100 (MQ=255)

                                                                                          cGCTAATGCCGGGTTACCTTGATGATAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGCGCTTGAACTGTa           >  1:1576448/1‑100 (MQ=255)

                                                                                             tcttGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCg        <  1:732846/97‑1 (MQ=255)

                                                                                             tAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCg        <  1:1754086/100‑1 (MQ=255)

                                                                                               aTGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTg      <  1:807488/100‑1 (MQ=255)

                                                                                                tGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTgg     <  2:1368211/100‑1 (MQ=255)

                                                                                                   cGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAAt  >  1:211582/1‑100 (MQ=255)

AATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAAT  >  NC_000913/3991202‑3991400

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