BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,150,120	T→C	26.0%	F153F (TTT→TTC)	yghA →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,150,120	0	T	C	26.0%	299.1 / 121.8	188	F153F (TTT→TTC)	yghA	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (53/85); new base C (23/26); total (77/111)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

CGAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTC  >  NC_000913/3150021‑3150218

cGAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCAATCCGGATATTGCAGACCTCACCAGCGAACAGttt                                                                                                    >  2:2176219/1‑100 (MQ=255)

 gAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTc                                                                                                   <  2:3120214/100‑1 (MQ=255)

 gAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACACTTTc                                                                                                   <  1:3621442/100‑1 (MQ=255)

   aGCGTACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCaa                                                                                                 <  2:3215612/100‑1 (MQ=255)

    gcgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCaaa                                                                                                <  2:2020317/100‑1 (MQ=255)

     cgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCaaaa                                                                                               <  1:2215144/100‑1 (MQ=255)

     cgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCaaaa                                                                                               >  1:3537576/1‑100 (MQ=255)

     cgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCaaaa                                                                                               >  2:3822545/1‑100 (MQ=255)

     cgcACAAGGCGTTAGGCGGGCTGGAGATTATGTCGATGTTCGAGGGGAAACAGGTTGCCATTCCGCATATTGCAGAACTCACCATCGAACAGTTCCaaaa                                                                                               <  2:548361/100‑1 (MQ=255)

      gcACAAGGCGTTAGGCGGGCTGGATATTATGTCGCTGGTCGCCGGGAAACAGGTTGCCATTGCGGATATGGCAGACCTCACCAGCGAACAGTTTCAAAAg                                                                                              >  1:188862/1‑100 (MQ=255)

      gcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAg                                                                                              >  2:1790655/1‑100 (MQ=255)

      gcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAg                                                                                              >  1:3886046/1‑100 (MQ=255)

       cacaAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGa                                                                                             <  2:277834/100‑1 (MQ=255)

        acaAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGAc                                                                                            <  2:1219212/100‑1 (MQ=255)

        acaAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGAc                                                                                            >  1:256538/1‑100 (MQ=255)

        acaAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGAc                                                                                            <  1:1835056/100‑1 (MQ=255)

        acaAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGAc                                                                                            <  1:262025/100‑1 (MQ=255)

         caAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGAcc                                                                                           <  1:860079/100‑1 (MQ=255)

           aGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCtt                                                                                         >  2:3746349/1‑100 (MQ=255)

           aGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCtt                                                                                         >  2:2338325/1‑100 (MQ=255)

           aGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCtt                                                                                         <  2:795302/100‑1 (MQ=255)

             gCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCttcg                                                                                       <  1:3822570/100‑3 (MQ=255)

                ttAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCa                                                                                    <  1:3002427/100‑1 (MQ=255)

                 tAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCAt                                                                                   <  1:2024482/100‑1 (MQ=255)

                   ggCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTa                                                                                 >  2:1400450/1‑100 (MQ=255)

                    gCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATCaa                                                                                >  2:3371814/1‑100 (MQ=255)

                      gggCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACg                                                                              <  2:1093981/100‑1 (MQ=255)

                       ggCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGt                                                                             >  1:992251/1‑100 (MQ=255)

                        gCTGGATATTATGGCGCTGGTCGTCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCAt                                                                                   <  2:695387/93‑1 (MQ=255)

                        gCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGtt                                                                            >  1:1015386/1‑100 (MQ=255)

                         cTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGttt                                                                           >  2:1040646/1‑100 (MQ=255)

                         cTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGttt                                                                           >  1:3913313/1‑100 (MQ=255)

                          tGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATCTTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTCTGCCATTAACGtatt                                                                          >  2:3309998/1‑100 (MQ=255)

                           ggATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCGGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTc                                                                         >  1:1309621/1‑100 (MQ=255)

                           ggAGATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATGCAGGATATTGCAGAGGTCACGAGCGAAGAGGTGCAAAAGACCTTTGCCATTAACGTTTTc                                                                         >  2:2248227/1‑100 (MQ=255)

                            gATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAATACCTTTGCCATTAACGTTTTcg                                                                        <  2:731172/100‑1 (MQ=255)

                            gATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTcg                                                                        <  1:2592496/100‑1 (MQ=255)

                            gATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATGAACGTTTTcg                                                                        <  2:2876231/100‑1 (MQ=255)

                             atatTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTcgc                                                                       <  2:188861/100‑1 (MQ=255)

                             atatTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTcgc                                                                       <  2:3886021/100‑1 (MQ=255)

                              tattatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTcgcg                                                                      <  2:353385/100‑1 (MQ=255)

                               attacGGCGCTGGTCGCCCGGAAACAGGTAGGCATGCCGGATACTGCAGACCTCACCAGCGAACAGTTCCAAACGACCTTTGCCATTAACGTTTTCgcgc                                                                     <  2:1927229/100‑1 (MQ=255)

                                ttatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGttt                                                                           >  2:3481482/1‑93 (MQ=255)

                                 tatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTg                                                                   <  2:528871/100‑1 (MQ=255)

                                 tatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTg                                                                   <  1:1790667/100‑1 (MQ=255)

                                 tatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTg                                                                   >  1:773232/1‑100 (MQ=255)

                                  atGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGt                                                                  >  2:3800683/1‑100 (MQ=255)

                                  atGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGt                                                                  >  2:3985356/1‑100 (MQ=255)

                                     gcgcTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCt                                                               >  2:3677699/1‑100 (MQ=255)

                                        cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGTAGACCTAACCAGCGAACAGTTTCAAAAGACCTATGCCATTAACGTTTTCGCGCTGttgaggc                                                            <  2:3930329/100‑6 (MQ=255)

                                        cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGc                                                            <  1:2429441/100‑1 (MQ=255)

                                        cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGc                                                            <  2:2361173/100‑1 (MQ=255)

                                               ccGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCa                                                     <  1:2284036/100‑1 (MQ=255)

                                               ccGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCa                                                     >  2:1421701/1‑100 (MQ=255)

                                                cGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAg                                                    <  1:1522967/100‑1 (MQ=255)

                                                     aaCAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGc                                               >  1:3250436/1‑100 (MQ=255)

                                                     aaCAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCGAGGAAGc                                               <  1:2159087/100‑1 (MQ=255)

                                                         ggTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATc                                           <  2:3137133/100‑1 (MQ=255)

                                                           ttGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCACAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATccc                                         <  2:403306/100‑1 (MQ=255)

                                                           ttGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATccc                                         >  1:2514478/1‑100 (MQ=255)

                                                            tGCCATTCCGGATATTGCAGACCTCGCCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATccc                                         <  1:2234156/99‑1 (MQ=255)

                                                             gCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATccccc                                       <  1:3685964/100‑1 (MQ=255)

                                                              ccATTCCGGCTATTGCAGACCTCACGCGCGAACAGTTTCAAAAGACGTTTGCCATTCACGTTTTCGCGCTGTTATGGCTAACCCAGTAAGCGATCCCCCt                                      >  2:3594874/1‑100 (MQ=255)

                                                              ccATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCt                                      <  1:560867/100‑1 (MQ=255)

                                                               cATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTg                                     <  1:1170297/100‑1 (MQ=255)

                                                                aTTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGc                                    <  1:4079972/100‑1 (MQ=255)

                                                                aTTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGc                                    <  1:909311/100‑1 (MQ=255)

                                                                aTTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGc                                    <  1:2268877/100‑1 (MQ=255)

                                                                 ttCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCt                                   <  2:1199758/100‑1 (MQ=255)

                                                                   ccGGATATTGCAGACCTCACCGGCGAACAGTTTCAAAAGACCGTTGCCATTAACGGTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGc                                    <  2:3253763/97‑1 (MQ=255)

                                                                    cGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAcc                                <  1:2728300/100‑1 (MQ=255)

                                                                           tGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGt                         <  1:625358/100‑1 (MQ=255)

                                                                           tGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGt                         >  2:4057016/1‑100 (MQ=255)

                                                                           tGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTCCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGt                         >  2:2213539/1‑100 (MQ=255)

                                                                            gCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTg                        <  1:3588589/100‑1 (MQ=255)

                                                                            gCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTg                        <  2:773226/100‑1 (MQ=255)

                                                                            gCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTg                        >  1:271635/1‑100 (MQ=255)

                                                                             cAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGc                       <  2:1665874/100‑1 (MQ=255)

                                                                              aGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCa                      <  2:2624176/100‑1 (MQ=255)

                                                                                aCCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                    <  1:3370162/100‑1 (MQ=255)

                                                                                aCCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                    <  1:3456259/100‑1 (MQ=255)

                                                                                aCCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                    <  1:3746373/100‑1 (MQ=255)

                                                                                aCCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                    <  1:3737170/100‑1 (MQ=255)

                                                                                 cttCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGt                   >  1:2804138/3‑100 (MQ=255)

                                                                                   tCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTa                  >  2:1955340/1‑99 (MQ=255)

                                                                                    cACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAATCGATCCCCCTGCTACCGAAAGGTGCAAGTatc                <  2:2586240/100‑1 (MQ=255)

                                                                                      ccAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGCTTTCGCGCTGTTCTGGCCAACCCAGGAAGCGATCCCCCTGCTACCGAAAGCTGCAAGTatcat              >  2:3585289/1‑100 (MQ=255)

                                                                                      ccAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatcat              >  2:1467479/1‑100 (MQ=255)

                                                                                        aGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTAtcatca            >  2:1598990/1‑100 (MQ=255)

                                                                                         gCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcac           >  1:1408345/1‑100 (MQ=255)

                                                                                          cGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacc          <  1:334288/100‑1 (MQ=255)

                                                                                          cGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacc          <  1:3985382/100‑1 (MQ=255)

                                                                                          cGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacc          <  1:2913216/100‑1 (MQ=255)

                                                                                             aCAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACt       <  2:3008187/100‑1 (MQ=255)

                                                                                              cAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACt       >  2:426498/1‑99 (MQ=255)

                                                                                                gTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCg    >  1:948614/1‑100 (MQ=255)

                                                                                                gTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCg    >  1:3197683/1‑100 (MQ=255)

                                                                                                gTTCCAAAAGACCTTTGCCATTGACGTTTTCGCGCGATGCTAGCTAACCCCGGCAGCGGTGCCACTGCTACCCAACGGCCAAAATCCCATCCCCACTTc     >  2:3395794/1‑99 (MQ=255)

                                                                                                  ttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTc  >  2:1569763/1‑100 (MQ=255)

                                                                                                  ttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTc  <  1:1019981/100‑1 (MQ=255)

CGAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTC  >  NC_000913/3150021‑3150218

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