BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,261,198	T→C	T205A (ACC→GCC)	qor ←	quinone oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,198	0	T	C	99.4%	1081.3 / ‑3.9	348	T205A (ACC→GCC)	qor	quinone oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (193/153); total (194/154)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

Only 100 of 356 total aligned reads displayed.

TTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCG  >  NC_000913/4261108‑4261287

ttAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGa                                                                                         >  1:3616173/1‑93 (MQ=255)

 tAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAt                                                                                        <  1:5938555/93‑1 (MQ=255)

  aaGCCGCGGCGTTGCAGGCAATCCAGCGTCCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAtc                                                                                       <  1:3352156/93‑1 (MQ=255)

    gCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAtctc                                                                                     >  1:2741240/1‑93 (MQ=255)

    gCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAtctc                                                                                     >  1:1044358/1‑93 (MQ=255)

    gCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAtctc                                                                                     <  2:6175492/93‑1 (MQ=255)

      cgcgGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCtt                                                                                   <  2:1286981/93‑1 (MQ=255)

      cgcgGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCtt                                                                                   >  2:7287101/1‑93 (MQ=255)

       gcggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTTCCGCCGGCGATCTCttt                                                                                  <  2:2149361/93‑1 (MQ=255)

       gcggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCttt                                                                                  <  1:4383739/93‑1 (MQ=255)

       gcggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCttt                                                                                  >  1:7168842/1‑93 (MQ=255)

        cggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCtttt                                                                                 <  2:1728113/93‑1 (MQ=255)

        cggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCtttt                                                                                 <  1:6698257/93‑1 (MQ=255)

         ggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTa                                                                                <  2:4157648/93‑1 (MQ=255)

         ggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTa                                                                                >  1:3013792/1‑93 (MQ=255)

          gcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTaa                                                                               >  1:3557820/1‑93 (MQ=255)

           cgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAAc                                                                              <  2:3060165/93‑1 (MQ=255)

             ttGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCg                                                                            >  2:657871/1‑93 (MQ=255)

             ttGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCg                                                                            >  2:6978528/1‑93 (MQ=255)

                 aGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCg                                                                        >  1:924865/1‑93 (MQ=255)

                   gCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGAc                                                                      <  1:5382130/93‑1 (MQ=255)

                     aaTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCa                                                                    >  1:2149835/1‑93 (MQ=255)

                     aaTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCa                                                                    <  1:1332195/93‑1 (MQ=255)

                      aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAg                                                                   >  2:2859264/1‑93 (MQ=255)

                      aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAg                                                                   >  2:4490964/1‑93 (MQ=255)

                       tCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGa                                                                  >  2:4688419/1‑93 (MQ=255)

                       tCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGa                                                                  >  1:1129755/1‑93 (MQ=255)

                        ccAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGAt                                                                 <  2:5882373/93‑1 (MQ=255)

                         cAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATc                                                                >  2:3276006/1‑93 (MQ=255)

                         cAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATc                                                                >  1:5221374/1‑93 (MQ=255)

                           gCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCct                                                              >  2:4840010/1‑93 (MQ=255)

                              aCCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCtt                                                           >  2:2395595/1‑93 (MQ=255)

                               ccGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTc                                                          >  2:141276/1‑93 (MQ=255)

                                cGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCa                                                         <  2:1816423/93‑1 (MQ=255)

                                 gtccccAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCAc                                                        >  1:2331605/4‑93 (MQ=255)

                                  tccccAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACg                                                       >  1:1285426/3‑93 (MQ=255)

                                  tccccAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACg                                                       >  1:7580904/3‑93 (MQ=255)

                                   ccccAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGa                                                      >  2:6007415/2‑93 (MQ=255)

                                   ccccAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGa                                                      >  1:6244122/2‑93 (MQ=255)

                                     ccAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATa                                                    >  1:287237/1‑93 (MQ=255)

                                     ccAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATa                                                    >  2:5016240/1‑93 (MQ=255)

                                      cAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAg                                                   <  2:3208679/93‑1 (MQ=255)

                                      cAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAg                                                   >  2:4550466/1‑93 (MQ=255)

                                      cAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAg                                                   >  1:2061716/1‑93 (MQ=255)

                                      cAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAg                                                   >  2:2458244/1‑93 (MQ=255)

                                       aGGTGTCTCTGCCCACTGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACTAGATCCTCTTCACGATAGt                                                  >  2:7017019/1‑93 (MQ=255)

                                       aGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGt                                                  <  1:1980333/93‑1 (MQ=255)

                                        ggTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGtt                                                 >  1:1755718/1‑93 (MQ=255)

                                          tgtCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTtaa                                               >  1:457220/1‑93 (MQ=255)

                                          tgtCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTtaa                                               <  1:3875283/93‑1 (MQ=255)

                                           gtCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGTCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTtaat                                              <  2:5297476/93‑1 (MQ=255)

                                           gtCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTtaat                                              <  1:3500314/93‑1 (MQ=255)

                                            tctcTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTtaata                                             <  1:7403713/93‑1 (MQ=255)

                                            tctcTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTtaata                                             <  2:5621733/93‑1 (MQ=255)

                                             gtctGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTtaataa                                            <  1:248369/92‑1 (MQ=255)

                                              tctGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAAc                                           <  1:4430167/93‑1 (MQ=255)

                                                tGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCt                                         <  1:6780015/93‑1 (MQ=255)

                                                tGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCt                                         <  1:4827045/93‑1 (MQ=255)

                                                 gCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTg                                        <  1:6084686/93‑1 (MQ=255)

                                                   ccACGGAATCGTACACCACGCGCACTTGCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGcc                                      >  1:816507/1‑93 (MQ=255)

                                                    cACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAATCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCa                                     >  2:4032381/1‑93 (MQ=255)

                                                     aCGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCAc                                    >  1:5348323/1‑93 (MQ=255)

                                                      cGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCAcg                                   <  1:382633/93‑1 (MQ=255)

                                                       ggAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCAcgc                                  >  1:3673248/1‑93 (MQ=255)

                                                        gAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCAcgcg                                 >  1:646752/1‑93 (MQ=255)

                                                        gAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCAcgcg                                 >  1:3453035/1‑93 (MQ=255)

                                                        gAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCAcgcg                                 >  1:3285605/1‑93 (MQ=255)

                                                         aaTCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACgcgc                                <  2:1805180/93‑1 (MQ=255)

                                                         aaTCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACgcgc                                >  2:1171580/1‑93 (MQ=255)

                                                         aaTCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACgcgc                                >  2:4376842/1‑93 (MQ=255)

                                                           tCGTACACCCCGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTGGACCAGATCCTCTTCGCGATCGTTAATACACTGCCACGCGccg                              >  1:6954356/1‑92 (MQ=255)

                                                           tCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGccc                              <  1:5731203/93‑1 (MQ=255)

                                                           tCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGccc                              >  1:166805/1‑93 (MQ=255)

                                                             gTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGc                            >  2:3031669/1‑93 (MQ=255)

                                                              tACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCt                           >  1:4632662/1‑93 (MQ=255)

                                                              tACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCt                           >  2:2269227/1‑93 (MQ=255)

                                                                ctccacGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCttt                         <  2:922997/91‑1 (MQ=255)

                                                                 accacGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCtttt                        >  2:1531448/1‑93 (MQ=255)

                                                                   cacGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAg                      >  1:6347175/1‑93 (MQ=255)

                                                                    acGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAgc                     >  2:514662/1‑93 (MQ=255)

                                                                       cgcACTTTCTTCCCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAgcgcg                  <  1:1565415/93‑1 (MQ=255)

                                                                       cgcACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAgcgcg                  >  1:7007150/1‑93 (MQ=255)

                                                                       cgcACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCGGCCACGCGCCCGCTTTTAgcgcg                  >  2:4298114/1‑93 (MQ=255)

                                                                        gcACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAgcgcgc                 >  2:3486424/1‑93 (MQ=255)

                                                                         cACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGct                >  1:2670782/1‑93 (MQ=255)

                                                                         cACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGct                <  2:7481886/93‑1 (MQ=255)

                                                                          aCTTTCTTACCGCCGGCGATCTCTTTTAGCCGCTCGACCAGACCCTCTGCACGATAGTTAAGAACCTGCCACGCGCCCGCTTTTAGCGCGctc               <  2:688772/93‑1 (MQ=255)

                                                                          aCTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGctc               >  1:261266/1‑93 (MQ=255)

                                                                          aCTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGctc               >  2:5495394/1‑93 (MQ=255)

                                                                           cTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGctct              >  1:7836461/1‑93 (MQ=255)

                                                                           cTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGctct              <  2:3532222/93‑1 (MQ=255)

                                                                            tttCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTg             >  1:2475673/1‑93 (MQ=255)

                                                                            tttCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTg             >  2:780223/1‑93 (MQ=255)

                                                                            tttCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTg             <  2:3685727/93‑1 (MQ=255)

                                                                                ttACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCt         >  1:4994918/1‑93 (MQ=255)

                                                                                  aCCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCttt       >  2:2001092/1‑93 (MQ=255)

                                                                                  aCCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCttt       >  2:2228931/1‑93 (MQ=255)

                                                                                   ccgccgGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCtttt      <  1:6167949/93‑1 (MQ=255)

                                                                                    cgccgGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCttttt     >  2:6102500/1‑93 (MQ=255)

                                                                                       cgGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTgcg  >  1:4919191/1‑93 (MQ=255)

TTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCG  >  NC_000913/4261108‑4261287

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