BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,100,361	C→T	5.5%	R334Q (CGG→CAG) ‡	acsA‑I ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,100,361	0	C	T	5.5%	923.5 / 27.2	289	R334Q (CGG→CAG) ‡	acsA‑I	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base C (124/149); new base T (4/12); total (128/161)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 302 total aligned reads displayed.

GGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTT  >  NC_002947/5100213‑5100509

ggTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCg                                                                                                                                                     >  2:178328/1‑150 (MQ=255)

     tGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCAt                                                                                                                                                <  2:3912350/150‑1 (MQ=255)

     tGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCAt                                                                                                                                                >  2:76213/1‑150 (MQ=255)

     tGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCAt                                                                                                                                                <  1:5394136/150‑1 (MQ=255)

      gTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATg                                                                                                                                               >  2:2402968/1‑150 (MQ=255)

      gTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCACGGACGGTAATCGAACACCCGTTCATg                                                                                                                                               >  2:1001339/1‑150 (MQ=255)

       tCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATgg                                                                                                                                              >  1:7999666/1‑150 (MQ=255)

       tCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATgg                                                                                                                                              >  1:4122289/1‑150 (MQ=255)

        ccGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATggg                                                                                                                                             <  1:70856/150‑1 (MQ=255)

         cGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGt                                                                                                                                            >  1:4320873/1‑150 (MQ=255)

           ggTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCa                                                                                                                                          <  2:9072970/150‑1 (MQ=255)

                ttGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGcgcg                                                                                                                                     >  2:6513570/1‑150 (MQ=255)

                ttGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGcgcg                                                                                                                                     >  2:3551984/1‑150 (MQ=255)

                  gggCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCgcggc                                                                                                                                   >  1:7485870/1‑150 (MQ=255)

                  gggCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCgcggc                                                                                                                                   >  1:3601494/1‑150 (MQ=255)

                         ccTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATgg                                                                                                                                              >  2:246854/1‑132 (MQ=255)

                          cTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCa                                                                                                                           >  1:3925648/1‑150 (MQ=255)

                           ttCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAg                                                                                                                          >  2:6126225/1‑150 (MQ=255)

                           ttCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAg                                                                                                                          >  1:1042766/1‑150 (MQ=255)

                                    cAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGcc                                                                                                                 >  1:3138531/1‑150 (MQ=255)

                                        gtggtgGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCACCCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCggtg                                                                                                             >  2:8509426/1‑150 (MQ=255)

                                        gtggtgGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCtgtg                                                                                                             <  2:3997213/150‑1 (MQ=255)

                                        gtggtgGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCtgtg                                                                                                             >  1:1544532/1‑150 (MQ=255)

                                        gtggtgGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCtgtg                                                                                                             <  2:1922563/150‑1 (MQ=255)

                                         tggtggCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTgtgg                                                                                                            >  2:3805846/1‑150 (MQ=255)

                                             ggCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGtgtg                                                                                                        <  1:3420696/150‑1 (MQ=255)

                                              gcgcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgc                                                                                                       >  1:6619224/1‑150 (MQ=255)

                                               cgcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgca                                                                                                      <  2:2035722/150‑1 (MQ=255)

                                               cgcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgca                                                                                                      <  2:4636303/150‑1 (MQ=255)

                                                gcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgcag                                                                                                     >  2:6961068/1‑150 (MQ=255)

                                                gcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgcag                                                                                                     >  2:8107541/1‑150 (MQ=255)

                                                 ccGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgcagc                                                                                                    >  1:5315436/1‑150 (MQ=255)

                                                  cGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgcagca                                                                                                   >  2:5680342/1‑150 (MQ=255)

                                                  cGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGAACCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAAACAGGTAACCGCCTGTGGTGTgcagca                                                                                                   >  2:1134101/1‑150 (MQ=255)

                                                     tGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGc                                                                                                >  1:7563552/1‑150 (MQ=255)

                                                                ccGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTg                                                                                     >  1:911357/1‑150 (MQ=255)

                                                                     gACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGt                                                                                <  1:7374498/150‑1 (MQ=255)

                                                                     gACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGt                                                                                >  1:5167475/1‑150 (MQ=255)

                                                                     gACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGt                                                                                <  2:618431/150‑1 (MQ=255)

                                                                      aCGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTg                                                                               <  2:7279008/150‑1 (MQ=255)

                                                                      aCGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTg                                                                               >  1:8419889/1‑150 (MQ=255)

                                                                        gATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGa                                                                             >  1:6202321/1‑150 (MQ=255)

                                                                               cTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAg                                                                      <  2:8155786/150‑1 (MQ=255)

                                                                               cTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAg                                                                      >  2:4463089/1‑150 (MQ=255)

                                                                               cTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAg                                                                      <  2:3206765/150‑1 (MQ=255)

                                                                                      cTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATaaa                                                               <  2:1331417/150‑1 (MQ=255)

                                                                                      cTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATaaa                                                               <  2:7218060/150‑1 (MQ=255)

                                                                                      cTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATaaa                                                               <  1:194780/150‑1 (MQ=255)

                                                                                      cTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATaaa                                                               <  1:6074443/150‑1 (MQ=255)

                                                                                       tgtgACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAg                                                              <  2:4904539/150‑1 (MQ=255)

                                                                                        gtgACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAgg                                                             >  1:9074602/1‑150 (MQ=255)

                                                                                         tgACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGa                                                            <  1:5912790/150‑1 (MQ=255)

                                                                                         tcaCCCAGCCCACATCCGCGGTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGa                                                            <  2:5859896/148‑1 (MQ=255)

                                                                                         tcaCCCAGCCCACATCCGCGGTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGa                                                            <  2:8815713/148‑1 (MQ=255)

                                                                                          caCCCAGCCCACATCCGCGGTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGAt                                                           <  2:1636829/149‑1 (MQ=255)

                                                                                             ccAGCCCACATCCGCGGTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGaa                                                        >  2:668923/1‑150 (MQ=255)

                                                                                             ccAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGaa                                                        >  2:578933/1‑150 (MQ=255)

                                                                                               aaCCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACa                                                      >  2:519624/1‑150 (MQ=255)

                                                                                                  cGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAgcg                                                   >  1:3497943/1‑150 (MQ=255)

                                                                                                   gACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAgcgc                                                  >  2:5270854/1‑150 (MQ=255)

                                                                                                   gACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAgcgc                                                  <  1:7051819/150‑1 (MQ=255)

                                                                                                     cGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCct                                                <  2:7410584/150‑1 (MQ=255)

                                                                                                        cggcggTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCtt                                             <  2:5602057/150‑1 (MQ=255)

                                                                                                         ggcggTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTc                                            <  2:7960860/150‑1 (MQ=255)

                                                                                                         ggcggTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTc                                            <  1:8222855/150‑1 (MQ=255)

                                                                                                          gcggTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGAACAgcgg                                                  >  1:4786666/1‑142 (MQ=255)

                                                                                                           cggTTCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCTGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAg                                          <  2:4562476/150‑1 (MQ=255)

                                                                                                                 aCCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCa                                    >  1:2026602/1‑150 (MQ=255)

                                                                                                                   cAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCAtt                                  >  1:2767946/1‑150 (MQ=255)

                                                                                                                      tACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCt                               >  1:8054702/1‑150 (MQ=255)

                                                                                                                        cacCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                             <  1:505496/150‑1 (MQ=255)

                                                                                                                        cacCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                             <  2:6849324/150‑1 (MQ=255)

                                                                                                                        cacCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                             <  1:7627014/150‑1 (MQ=255)

                                                                                                                         acCTCGCCCGGACGGTCATCGAACACCCGTTCCTGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCGCTTCAGCACCCCTTTCtggg                            >  2:5908434/1‑147 (MQ=255)

                                                                                                                         acCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTg                                                                               >  2:772616/1‑99 (MQ=255)

                                                                                                                          ccTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTgg                           <  2:7506799/150‑1 (MQ=255)

                                                                                                                          ccTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTgg                           <  1:3760077/150‑1 (MQ=255)

                                                                                                                           cTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTggg                          <  1:881595/150‑1 (MQ=255)

                                                                                                                            tCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAg                                                                      <  2:2232232/105‑1 (MQ=255)

                                                                                                                            tCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAg                                                                      >  1:2232232/1‑105 (MQ=255)

                                                                                                                                ccGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                             <  2:1725366/142‑1 (MQ=255)

                                                                                                                                ccGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                             >  1:1725366/1‑142 (MQ=255)

                                                                                                                                 cGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCa                    <  1:452761/150‑1 (MQ=255)

                                                                                                                                  ggACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAg                   >  1:2378322/1‑150 (MQ=255)

                                                                                                                                        tAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTg             >  2:3349171/1‑150 (MQ=255)

                                                                                                                                        tAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTg             <  1:7746919/150‑1 (MQ=255)

                                                                                                                                         aaTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTgg            >  2:1229791/1‑150 (MQ=255)

                                                                                                                                         aaTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTgg            >  2:4733241/1‑150 (MQ=255)

                                                                                                                                         aaTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTgg            >  1:2475093/1‑150 (MQ=255)

                                                                                                                                         aaTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTgg            <  1:2403157/150‑1 (MQ=255)

                                                                                                                                          aTCTAACACCCGTTCATTGTTCATCGGGGCTTAAACAATGTAAGCTCCTGTGGTTTGCATCTCTCCCTTCGGCTTTCCGGTGGTGAGGGAGGTATAATGGTTGATCAGCGCCTCTCGATCACCCATTTCTTTGGGGGCGCAGTGGCTgta           <  1:4677743/150‑3 (MQ=255)

                                                                                                                                          aTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGTGGGCGCAGTGGCTGGa           <  1:1057766/150‑1 (MQ=255)

                                                                                                                                          aTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGa           >  1:6964271/1‑150 (MQ=255)

                                                                                                                                          aTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGa           >  1:295769/1‑150 (MQ=255)

                                                                                                                                            cGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgc                                                                                                       <  1:8284413/56‑1 (MQ=255)

                                                                                                                                              aaCACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGcc       >  2:4156167/1‑150 (MQ=255)

                                                                                                                                              aaCACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGcc       >  2:4064600/1‑150 (MQ=255)

                                                                                                                                               acacCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCa      <  2:2853755/150‑1 (MQ=255)

                                                                                                                                                   ccGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACttt  >  2:8327118/1‑150 (MQ=255)

                                                                                                                                                   ccGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACttt  >  1:1314730/1‑150 (MQ=255)

GGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTT  >  NC_002947/5100213‑5100509

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