BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	571,696	T→C	L3P (CTC→CCC)	ybcN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	571,696	0	T	C	99.1%	1049.3 / ‑3.9	339	L3P (CTC→CCC)	ybcN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/3); new base C (191/145); total (191/148)
Fisher's exact test for biased strand distribution p-value = 8.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 354 total aligned reads displayed.

CAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTC  >  NC_000913/571607‑571785

cAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCcctc                                                                                        >  1:660084/1‑93 (MQ=255)

    ctctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACaa                                                                                    >  2:1531039/1‑93 (MQ=255)

    ctctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACaa                                                                                    <  1:1834157/93‑1 (MQ=255)

     tctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAg                                                                                   >  1:3792412/1‑93 (MQ=255)

      ctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGa                                                                                  <  1:5691943/93‑1 (MQ=255)

      ctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGa                                                                                  <  2:6497023/93‑1 (MQ=255)

        cttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACg                                                                                >  2:7722290/1‑93 (MQ=255)

        cttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACg                                                                                >  1:2905735/1‑93 (MQ=255)

         ttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACgg                                                                               >  2:5506476/1‑93 (MQ=255)

         ttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACgg                                                                               >  2:3998707/1‑93 (MQ=255)

         tactcACTTCGCGGAAGCAGAAATGCAAATGCTATTCTTTGCATCTAGCGATTGACCTCTGCGGAGGGGAGGCGTACACACCTAACAAGACga                                                                               <  2:7014211/91‑2 (MQ=255)

           ctcCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCa                                                                             <  2:5288132/93‑1 (MQ=255)

            tcCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCAt                                                                            >  1:7427473/1‑93 (MQ=255)

            tcCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCAt                                                                            <  2:3152817/93‑1 (MQ=255)

             cccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATc                                                                           >  2:6274368/1‑93 (MQ=255)

              ccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCa                                                                          >  2:5881/1‑93 (MQ=255)

              ccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCa                                                                          >  2:2818431/1‑93 (MQ=255)

                ttCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCaaa                                                                        >  2:285974/1‑93 (MQ=255)

                ttCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCaaa                                                                        >  1:1279201/1‑93 (MQ=255)

                 tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAAt                                                                       >  1:7214627/1‑93 (MQ=255)

                    caGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTAc                                                                    >  2:5820510/1‑93 (MQ=255)

                     aGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACa                                                                   <  2:3081542/93‑1 (MQ=255)

                     aGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACa                                                                   <  2:4442199/93‑1 (MQ=255)

                      gATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACAt                                                                  >  1:2352984/1‑93 (MQ=255)

                      gATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACAt                                                                  <  2:4156713/93‑1 (MQ=255)

                       aTGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATc                                                                 <  2:3209928/93‑1 (MQ=255)

                        tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATcg                                                                >  2:745734/1‑93 (MQ=255)

                         gCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCGAATTACATcgc                                                               <  2:1180196/93‑1 (MQ=255)

                           aGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCgg                                                             >  1:5607850/1‑93 (MQ=255)

                            gAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGc                                                            >  2:1393081/1‑93 (MQ=255)

                            gAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGc                                                            >  2:4292413/1‑93 (MQ=255)

                              aaTGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCaa                                                          >  2:1809369/1‑93 (MQ=255)

                              aaTGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCaa                                                          >  2:6768155/1‑93 (MQ=255)

                               aTGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAAc                                                         <  1:872614/93‑1 (MQ=255)

                                tGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACt                                                        >  1:4803055/1‑93 (MQ=255)

                                   aCCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCa                                                     >  1:4107135/1‑93 (MQ=255)

                                   aCCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCa                                                     >  1:7813648/1‑93 (MQ=255)

                                   aCCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCa                                                     <  1:169547/93‑1 (MQ=255)

                                    ccTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCAc                                                    <  1:7611927/93‑1 (MQ=255)

                                     cTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCAcc                                                   <  1:6940380/93‑1 (MQ=255)

                                      tGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCg                                                  >  2:4603868/1‑93 (MQ=255)

                                      tGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCg                                                  >  2:2999580/1‑93 (MQ=255)

                                      tGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCg                                                  >  1:3452972/1‑93 (MQ=255)

                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGc                                                 <  2:6345532/93‑1 (MQ=255)

                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGc                                                 <  2:389968/93‑1 (MQ=255)

                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGc                                                 <  1:2430262/93‑1 (MQ=255)

                                        cAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCt                                                >  1:6176946/1‑93 (MQ=255)

                                         aaTTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTa                                               >  1:4696227/1‑93 (MQ=255)

                                          aTTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAt                                              >  2:6581828/1‑93 (MQ=255)

                                           ttCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtc                                             >  2:5000057/1‑93 (MQ=255)

                                             cTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGGCGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcgg                                           >  2:6878762/1‑93 (MQ=255)

                                              tttGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggt                                          >  1:6441132/1‑93 (MQ=255)

                                              tttGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggt                                          >  1:2501518/1‑93 (MQ=255)

                                              tttGCATCTAGCGTTTGACCTCTGCGGAGGGAAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggt                                          <  1:2123503/93‑1 (MQ=255)

                                                tGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcg                                        <  2:4546624/93‑1 (MQ=255)

                                                   aTCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGTGGCAACTTCACCGCTATCGTTCGGCa                                     >  1:2188565/1‑93 (MQ=255)

                                                     cTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGa                                   <  1:316839/93‑1 (MQ=255)

                                                      tAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                  >  1:4343434/1‑93 (MQ=255)

                                                      tAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                  >  2:5360234/1‑93 (MQ=255)

                                                      tAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                  >  1:3052692/1‑93 (MQ=255)

                                                      tAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                  >  2:721741/1‑93 (MQ=255)

                                                       aGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATc                                 >  2:1050609/1‑93 (MQ=255)

                                                          gTTTGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAg                              >  1:3004702/1‑93 (MQ=255)

                                                           tttGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGc                             >  1:3558510/1‑93 (MQ=255)

                                                           tttGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGc                             <  1:3378420/93‑1 (MQ=255)

                                                           tttGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGc                             >  2:4730874/1‑93 (MQ=255)

                                                            ttGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                            <  2:5271579/93‑1 (MQ=255)

                                                             tGACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCt                           <  1:1350163/93‑1 (MQ=255)

                                                              gACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCCAATTACATCGCGGCAACTTCACCGCTATCGTACGGCAGATCCAGTCtt                          >  2:3748481/1‑93 (MQ=255)

                                                              gACCTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCtt                          >  1:6488166/1‑93 (MQ=255)

                                                               aCCTCTGCGGAGGGGAAGCGAGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCAGCGCTATCGGTCGGCAGATACAGCCTTa                         >  1:3941454/1‑93 (MQ=255)

                                                                ccTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTAt                        <  2:7908934/93‑1 (MQ=255)

                                                                ccTCTGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTAt                        <  2:705433/93‑1 (MQ=255)

                                                                 ctctGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATc                       >  1:2638668/1‑93 (MQ=255)

                                                                 ctctGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATc                       <  2:2139362/93‑1 (MQ=255)

                                                                 ctctGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATc                       <  2:623432/93‑1 (MQ=255)

                                                                 cgctGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATc                       <  2:6518390/91‑1 (MQ=255)

                                                                  tctGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCt                      >  1:4931780/1‑93 (MQ=255)

                                                                    tGCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTgg                    >  1:1454216/1‑93 (MQ=255)

                                                                     gCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTgga                   <  1:1107960/93‑1 (MQ=255)

                                                                     gCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTgga                   >  1:5774038/1‑93 (MQ=255)

                                                                     gCGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTgga                   <  2:6091341/93‑1 (MQ=255)

                                                                      cGGAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggag                  <  1:152619/93‑1 (MQ=255)

                                                                       ggAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggagg                 >  1:4584395/1‑93 (MQ=255)

                                                                        gAGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTACCTcgagga                >  1:1036802/1‑93 (MQ=255)

                                                                         aGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGgaggag               <  1:5252343/93‑1 (MQ=255)

                                                                         aGGGGAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGgaggag               >  1:3591495/1‑93 (MQ=255)

                                                                          ggggAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGaggagg              >  1:607936/1‑93 (MQ=255)

                                                                          ggggAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGaggagg              <  1:3697294/93‑1 (MQ=255)

                                                                          ggggAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGaggagg              <  2:2708348/93‑1 (MQ=255)

                                                                          ggggAAGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGaggagg              >  1:3267785/1‑93 (MQ=255)

                                                                               aGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCaaa         <  1:6138113/93‑1 (MQ=255)

                                                                               aGCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCaaa         <  2:2237565/93‑1 (MQ=255)

                                                                                gCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAAt        >  2:3345913/1‑93 (MQ=255)

                                                                                gCGTGAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAAt        <  2:2529220/93‑1 (MQ=255)

                                                                                    gAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCtt    >  1:5133283/1‑93 (MQ=255)

                                                                                    gAACCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCtt    <  2:3634947/93‑1 (MQ=255)

                                                                                     aaCCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCttt   <  1:3866799/93‑1 (MQ=255)

                                                                                     aaCCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCttt   <  2:737896/93‑1 (MQ=255)

                                                                                     aaCCCCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCttt   <  2:2397641/93‑1 (MQ=255)

CAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTC  >  NC_000913/571607‑571785

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