BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,764,907	A→T	5.5%	N62I (AAC→ATC)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,764,907	0	A	T	5.5%	664.4 / 20.5	202	N62I (AAC→ATC)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (74/115); new base T (7/4); total (83/119)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 210 total aligned reads displayed.

CCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCA  >  NC_002947/3764758‑3765056

ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGaa                                                                                                                                                        <  1:591889/150‑1 (MQ=255)

ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGaa                                                                                                                                                        <  1:4410827/150‑1 (MQ=255)

   gATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACAt                                                                                                                                                     >  1:4521478/1‑150 (MQ=255)

      ggTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGAGAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGAt                                                                                                                                                  <  2:3699963/150‑1 (MQ=255)

       gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATc                                                                                                                                                 <  1:2855374/150‑1 (MQ=255)

       gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATc                                                                                                                                                 >  1:1962106/1‑150 (MQ=255)

       gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATc                                                                                                                                                 <  1:3878943/150‑1 (MQ=255)

        tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATcc                                                                                                                                                <  1:1835955/150‑1 (MQ=255)

        tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATcc                                                                                                                                                >  1:1236015/1‑150 (MQ=255)

          aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCg                                                                                                                                              >  2:1665379/1‑150 (MQ=255)

          aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCg                                                                                                                                              <  2:2212449/150‑1 (MQ=255)

          aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCg                                                                                                                                              <  1:3401732/150‑1 (MQ=255)

          aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCg                                                                                                                                              <  1:2545648/150‑1 (MQ=255)

           gCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCgg                                                                                                                                             <  2:3852903/150‑1 (MQ=255)

           gCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCgg                                                                                                                                             >  2:3706547/1‑150 (MQ=255)

            cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGATCATGTTCCCGGt                                                                                                                                            <  1:1248565/150‑1 (MQ=255)

            cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGt                                                                                                                                            <  1:3323957/150‑1 (MQ=255)

              tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACACGGCGGAAAAGATCTAGCGCGCCTGGGCGTTCGACGTGACCGGGACCATGATCCCGGTcc                                                                                                                                          >  2:3098515/1‑149 (MQ=255)

              tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCa                                                                                                                                          >  2:4444620/1‑150 (MQ=255)

               cAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCAc                                                                                                                                         >  1:3931570/1‑150 (MQ=255)

               cAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCAc                                                                                                                                         >  1:3186023/1‑150 (MQ=255)

               cAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCAc                                                                                                                                         >  1:917184/1‑150 (MQ=255)

               cAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCAc                                                                                                                                         >  1:2112241/1‑150 (MQ=255)

                aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACg                                                                                                                                        <  2:3052327/150‑1 (MQ=255)

                 aCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAGGCGCGCCTGGGCGTTCGACGTGACCGGGACCATGATCCCGGTCACgg                                                                                                                                       >  2:635875/1‑150 (MQ=255)

                 aCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACgg                                                                                                                                       <  2:3502628/150‑1 (MQ=255)

                     aGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGa                                                                                                                                   <  1:2134929/150‑1 (MQ=255)

                       gcgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAgc                                                                                                                                 <  1:3027849/150‑1 (MQ=255)

                        cgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAgcg                                                                                                                                <  1:2128240/150‑1 (MQ=255)

                           ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCt                                                                                                                            <  2:1018241/150‑1 (MQ=255)

                              cGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTggn                                                                                                                         <  1:2908054/150‑2 (MQ=255)

                              cGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGc                                                                                                                         <  2:4054068/150‑1 (MQ=255)

                               ggCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCg                                                                                                                        <  2:1189822/150‑1 (MQ=255)

                                 ccTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGa                                                                                                                      <  2:3925537/150‑1 (MQ=255)

                                 ccTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGa                                                                                                                      <  2:688776/150‑1 (MQ=255)

                                 ccTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGa                                                                                                                      <  2:1394741/150‑1 (MQ=255)

                                   tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATc                                                                                                                    >  1:4457877/1‑150 (MQ=255)

                                    gtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCa                                                                                                                   <  1:1422214/150‑1 (MQ=255)

                                       gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGc                                                                                                                >  2:1930817/1‑150 (MQ=255)

                                         ccGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCg                                                                                                              >  1:860252/1‑150 (MQ=255)

                                             ggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGtgct                                                                                                          >  2:3423885/1‑150 (MQ=255)

                                               agCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGtgctgc                                                                                                        >  1:4253182/1‑150 (MQ=255)

                                                gCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCa                                                                                                       <  2:1853942/150‑1 (MQ=255)

                                                 cGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCaa                                                                                                      <  1:412397/150‑1 (MQ=255)

                                                  gACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAAc                                                                                                     >  2:3182668/1‑150 (MQ=255)

                                                    cGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCa                                                                                                   <  1:1465450/150‑1 (MQ=255)

                                                        gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTa                                                                                               >  2:1466058/1‑150 (MQ=255)

                                                         cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTAc                                                                                              >  1:1066701/1‑150 (MQ=255)

                                                                gAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCt                                                                                       <  2:2077956/150‑1 (MQ=255)

                                                                 aGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTg                                                                                      <  2:860252/150‑1 (MQ=255)

                                                                  ggCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGc                                                                                     >  1:537927/1‑150 (MQ=255)

                                                                   gCCGTGTACGGCCACGGCGCCGAGGATACCGCCTACAAGGCGGAATAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGAGCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCg                                                                                    <  2:3313533/150‑1 (MQ=255)

                                                                   gCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCg                                                                                    <  2:1265653/150‑1 (MQ=255)

                                                                    ccGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCgg                                                                                   <  1:1085323/150‑1 (MQ=255)

                                                                     cgtacaCGGCCTCGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGc                                                                                  >  1:4029674/6‑150 (MQ=255)

                                                                      gtgtACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGGCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCGGCGGCa                                                                                 >  2:1277305/1‑150 (MQ=255)

                                                                        gtACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACt                                                                               <  1:827465/150‑1 (MQ=255)

                                                                          aCGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATACCCGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                                             >  2:185883/1‑150 (MQ=255)

                                                                             gCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                          <  2:2106571/150‑1 (MQ=255)

                                                                              ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                                         <  1:4334888/150‑1 (MQ=255)

                                                                              ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                                         >  2:475250/1‑150 (MQ=255)

                                                                                aCGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgc                                                                       <  1:1653923/150‑1 (MQ=255)

                                                                                      ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTgg                                                                 >  2:3233523/1‑150 (MQ=255)

                                                                                            aTAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                                           <  1:1281962/150‑1 (MQ=255)

                                                                                            aTAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGACCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                                           <  1:1675503/150‑1 (MQ=255)

                                                                                              aaCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACggg                                                         <  1:3510946/150‑1 (MQ=255)

                                                                                                cGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGcc                                                       <  2:1075131/150‑1 (MQ=255)

                                                                                                   cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                    <  2:2620290/150‑1 (MQ=255)

                                                                                                   cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                    <  1:3290118/150‑1 (MQ=255)

                                                                                                   cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                    <  2:2968219/150‑1 (MQ=255)

                                                                                                   cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                    <  1:4126210/150‑1 (MQ=255)

                                                                                                    tACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTgg                                                   >  2:2753384/1‑150 (MQ=255)

                                                                                                      cAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACCTGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACGTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTc                                                 <  2:3199182/150‑1 (MQ=255)

                                                                                                          gCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAg                                             <  1:2608796/150‑1 (MQ=255)

                                                                                                           cGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgc                                            <  2:3763130/150‑1 (MQ=255)

                                                                                                                   aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGAc                                    <  2:2442929/150‑1 (MQ=255)

                                                                                                                      aaGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAg                                 <  2:2990580/150‑1 (MQ=255)

                                                                                                                      aaGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAAACACTTACCCCGCCTGCGGCACATGCCCTGCTGCTGGACAACGGGCCTTGGTAACAGCGCCCGGACAAg                                 >  1:2151905/1‑150 (MQ=255)

                                                                                                                         cgcgcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCaa                              <  1:564547/150‑1 (MQ=255)

                                                                                                                          gcgcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAt                             <  1:2192793/150‑1 (MQ=255)

                                                                                                                            gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                           <  2:4486048/150‑1 (MQ=255)

                                                                                                                            gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                           <  1:844206/150‑1 (MQ=255)

                                                                                                                            gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                           <  2:1962106/150‑1 (MQ=255)

                                                                                                                            gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                           <  2:2929142/150‑1 (MQ=255)

                                                                                                                            gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                           <  2:4388952/150‑1 (MQ=255)

                                                                                                                               tGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTc                        >  2:568823/1‑150 (MQ=255)

                                                                                                                               tGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTc                        >  2:739487/1‑150 (MQ=255)

                                                                                                                               tGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTc                        <  1:3628950/150‑1 (MQ=255)

                                                                                                                                    gTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCa                   >  1:2719617/1‑150 (MQ=255)

                                                                                                                                        gACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt               >  1:4368142/1‑150 (MQ=255)

                                                                                                                                        accGTGACCGGGATCATGTTCCCGGTCAC‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATc              >  1:3737310/3‑150 (MQ=255)

                                                                                                                                         ccGTGACCGGGATCATGTTCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCACTGGGTCAATCAGGATc              >  2:3621535/2‑150 (MQ=255)

                                                                                                                                           gTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGc            <  1:635875/150‑1 (MQ=255)

                                                                                                                                             gACCGGGATCATGTTCCCGGTCAC‑GTCGAACG‑CCCGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAcc         >  1:3775097/1‑150 (MQ=255)

                                                                                                                                             gACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAc          >  1:2745155/1‑150 (MQ=255)

                                                                                                                                             gACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAc          >  1:2981882/1‑150 (MQ=255)

                                                                                                                                              aCCGGGATCATGTTCCCGGTCAC‑GTCGAACGCCCAGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAcc         >  1:761583/1‑150 (MQ=38)

                                                                                                                                                cGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                                                    <  1:2289875/89‑1 (MQ=255)

                                                                                                                                                cGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                                                    >  2:2289875/1‑89 (MQ=255)

                                                                                                                                                    aaCATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTc   <  1:1172757/150‑1 (MQ=255)

CCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCA  >  NC_002947/3764758‑3765056

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: `‑`