BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,464,993	A→T	8.8%	L58Q (CTG→CAG)	lipB ←	octanoyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,464,993	0	A	T	8.8%	305.8 / 17.7	102	L58Q (CTG→CAG)	lipB	octanoyltransferase
Reads supporting (aligned to +/- strand): ref base A (44/49); new base T (4/5); total (48/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 105 total aligned reads displayed.

CAGCTCGATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCAAGCTCGC  >  NC_002947/5464845‑5465140

cAGCTCGATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACcagcag                                                                                                                                                    >  2:1171217/1‑150 (MQ=255)

  gCTCGATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGt                                                                                                                                                  <  1:951573/150‑1 (MQ=255)

    tCGATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGc                                                                                                                                                <  2:1908772/150‑1 (MQ=255)

      gATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTACTc                                                                                                                                              >  2:485626/1‑150 (MQ=255)

       aTACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCATGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCg                                                                                                                                             <  2:1158900/150‑1 (MQ=255)

         aCGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCGAGCACCTGCTGGTGccgggc                                                                                                                                           <  1:141590/150‑6 (MQ=255)

           ggCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCTGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCt                                                                                                                                         >  2:1449017/1‑150 (MQ=255)

           ggCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCt                                                                                                                                         >  1:1200023/1‑150 (MQ=255)

            gCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCtt                                                                                                                                        >  2:1207517/1‑150 (MQ=255)

               gACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGcc                                                                                                                                     >  1:1512888/1‑150 (MQ=255)

                aCCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCg                                                                                                                                    <  2:931973/150‑1 (MQ=255)

                aCCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCg                                                                                                                                    <  1:1527154/150‑1 (MQ=255)

                  cAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGc                                                                                                                                  >  2:256095/1‑150 (MQ=255)

                  cAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGc                                                                                                                                  >  2:304605/1‑150 (MQ=255)

                  cAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGc                                                                                                                                  >  1:1827009/1‑150 (MQ=255)

                     cTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTg                                                                                                                               >  2:1180349/1‑150 (MQ=255)

                      tCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGc                                                                                                                              >  1:391292/1‑150 (MQ=255)

                           cacgcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGTGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTg                                                                                                                         >  2:616480/1‑150 (MQ=255)

                           cacgcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTg                                                                                                                         >  2:618432/1‑150 (MQ=255)

                            acgcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTgg                                                                                                                        <  1:100671/150‑1 (MQ=255)

                            acgcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTgg                                                                                                                        <  1:1207517/150‑1 (MQ=255)

                              gcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCTGCTGGTGccgggc                                                                                                                                           <  2:283663/129‑6 (MQ=255)

                              gcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCTGCTGGTGccgggc                                                                                                                                           >  1:283663/1‑124 (MQ=255)

                                cAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCTGCTGGTGccgg                                                                                                                                             <  1:1250091/125‑4 (MQ=255)

                                cAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCTGCTGGTGccgg                                                                                                                                             <  1:1250040/125‑4 (MQ=255)

                                cAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCTGCTGGTGccgg                                                                                                                                             >  2:1250091/1‑122 (MQ=255)

                                cAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGa                                                                                                                    >  2:721852/1‑150 (MQ=255)

                                cAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGa                                                                                                                    <  1:372598/150‑1 (MQ=255)

                                cAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGa                                                                                                                    >  2:96108/1‑150 (MQ=255)

                                 aaaGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGaa                                                                                                                   <  2:1526186/150‑1 (MQ=255)

                                     cccAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCACCGCGGTCGGTCTGCACCACCGGAATGACCCCAGGCACCGCCAGATGATTGGACAAGCCTGCCGTCCCCTGGGTGAAGAcc                                                                                                               >  2:681894/1‑150 (MQ=255)

                                      ccAGCCGGCGCACATCCAGCAGCAGGTAGGCCCACAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCCCCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGCCCTTGCCGGCCTGCCCCTGGGTGAAGACCg                                                                                                              <  2:877911/150‑1 (MQ=255)

                                      ccAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCg                                                                                                              >  1:996323/1‑150 (MQ=255)

                                        aGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCa                                                                                                            <  2:755495/150‑1 (MQ=255)

                                        aGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCa                                                                                                            >  1:536328/1‑150 (MQ=255)

                                          ccGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGTGTGAAGACCGCAgg                                                                                                          <  1:1013574/150‑1 (MQ=255)

                                          ccGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAgg                                                                                                          <  2:1682351/150‑1 (MQ=255)

                                                cacaTCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTc                                                                                                    <  1:1018207/150‑1 (MQ=255)

                                                cacaTCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCACTGATTCGCGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTc                                                                                                    >  1:1279901/1‑150 (MQ=255)

                                                       agcagcagGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGc                                                                                                                                  <  1:1274745/113‑1 (MQ=255)

                                                       agcagcagGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGc                                                                                                                                  >  2:1274745/1‑113 (MQ=255)

                                                           gcagGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTGGGCCTGGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGa                                                                                         >  2:990754/1‑150 (MQ=255)

                                                             agGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGAtt                                                                                       >  2:786192/1‑150 (MQ=255)

                                                                 aGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTNCAt                                                                                   <  2:1200023/150‑1 (MQ=255)

                                                                  ggCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATc                                                                                  <  1:307692/150‑1 (MQ=255)

                                                                  ggCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATc                                                                                  <  1:353100/150‑1 (MQ=255)

                                                                    ccaccaGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCt                                                                                >  2:428065/1‑150 (MQ=255)

                                                                     caccaGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCCTCCTg                                                                               >  2:1151397/1‑150 (MQ=255)

                                                                           cTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTc                                                                         <  1:1284378/150‑1 (MQ=255)

                                                                            tGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTcc                                                                        >  2:709890/1‑150 (MQ=255)

                                                                               ccGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGAcca                                                                                               >  1:1182877/1‑124 (MQ=255)

                                                                               ccGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGAcca                                                                                               <  2:1182877/124‑1 (MQ=255)

                                                                               ccGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCggg                                                                     >  1:340566/1‑150 (MQ=255)

                                                                                cgtggCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGgc                                                                    <  1:616480/147‑1 (MQ=255)

                                                                                cGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAAACTGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGgc                                                                    <  1:618432/150‑1 (MQ=255)

                                                                                    gCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGgctgcg                                                                >  2:1581744/1‑150 (MQ=255)

                                                                                    gCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGgctgcg                                                                <  2:979215/150‑1 (MQ=255)

                                                                                       aTGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGaa                                                                                                                   >  2:1915203/1‑96 (MQ=255)

                                                                                       aTGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGaa                                                                                                                   <  1:1915203/96‑1 (MQ=255)

                                                                                       aTGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGgctgcgctg                                                             >  1:846890/1‑150 (MQ=255)

                                                                                       aTGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGgctgcgctg                                                             <  1:24656/150‑1 (MQ=255)

                                                                                       aTGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTATCCGGgctgcgctg                                                             >  2:494884/1‑150 (MQ=255)

                                                                                        tGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGgctgcgctgc                                                            <  1:1188439/150‑1 (MQ=255)

                                                                                         ggTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTATCCGGGCTGCgctgct                                                           <  1:494884/150‑1 (MQ=255)

                                                                                           tAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCg                                                         <  1:410748/150‑1 (MQ=255)

                                                                                                 aCCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGAAGGTGCTCGGCCTTGCCGGCCTGCCCGTGGGTGAAGACCGCAGGGTGCTCGAACAGCCAGCTTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTaaaa                                                   <  1:1382468/150‑1 (MQ=255)

                                                                                                   cTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGAACCGAAAGAAATCTTTCTAACTTTCTGGGCTGCGCTGCTCGGTAAGacg                                                 >  2:496215/1‑150 (MQ=255)

                                                                                                   cTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCCGCCTGCGCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATCTCATGCTGGCTGTCCGCGCTGCGCTGCTCGGTAAAacg                                                 >  1:1039022/1‑150 (MQ=255)

                                                                                                     ggCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGGTTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAacgac                                               <  1:96108/150‑1 (MQ=255)

                                                                                                         gccgcGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCAt                                           <  2:391292/150‑1 (MQ=255)

                                                                                                          ccgcGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATg                                          <  2:536328/150‑1 (MQ=255)

                                                                                                           cgcgGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCATGCCCGTTGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATgg                                         <  2:1492268/150‑1 (MQ=255)

                                                                                                               gtcggtcTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCCTCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTc                                     <  1:1581744/150‑1 (MQ=255)

                                                                                                                 cggtcTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCa                                   <  1:990754/150‑1 (MQ=255)

                                                                                                                    tcTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGcc                                >  2:486915/1‑149 (MQ=255)

                                                                                                                    tcTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCa                                <  1:428065/150‑1 (MQ=255)

                                                                                                                    tcTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCa                                >  1:674386/1‑150 (MQ=255)

                                                                                                                     cTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCAc                               <  1:256095/150‑1 (MQ=255)

                                                                                                                      tGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGgctgcgct                                                              >  1:992036/1‑118 (MQ=255)

                                                                                                                        cacnacCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCgg                            >  2:1197232/1‑150 (MQ=255)

                                                                                                                        caccacCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCgg                            <  2:272063/150‑1 (MQ=255)

                                                                                                                         accaccGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCtt                                                                                                                                        <  1:1016330/41‑1 (MQ=255)

                                                                                                                         accaccGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCtt                                                                                                                                        >  2:1016330/1‑41 (MQ=255)

                                                                                                                         accaccGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCtt                                                                                                                                        <  1:75391/41‑1 (MQ=255)

                                                                                                                         accaccGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCtt                                                                                                                                        >  2:75391/1‑41 (MQ=255)

                                                                                                                            accGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGAccagcca                                                                                           <  1:218901/83‑1 (MQ=255)

                                                                                                                            accGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGAccagcca                                                                                           >  2:218901/1‑83 (MQ=255)

                                                                                                                              cGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATa                      <  2:492921/150‑1 (MQ=255)

                                                                                                                              cGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGATCAATTGTCCAGATTTCATCATCGCTAACCATACAGCGCATGTTTTTAAAACGACGCATGGCCTCCAGCACCGGTTCATa                      >  1:425361/1‑150 (MQ=255)

                                                                                                                               ggAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAg                     <  2:1827009/150‑1 (MQ=255)

                                                                                                                               ggAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAg                     <  2:1512888/150‑1 (MQ=255)

                                                                                                                                     tCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGc               <  1:721852/150‑1 (MQ=255)

                                                                                                                                       gCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAg             >  1:1487813/1‑150 (MQ=255)

                                                                                                                                            gCACCTGCTGGTGC‑CGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAgg            <  1:1817857/145‑1 (MQ=255)

                                                                                                                                            gCACCTGCTGGTGC‑CGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAgg            >  2:1817857/1‑145 (MQ=255)

                                                                                                                                            gCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCaa        >  1:954959/1‑150 (MQ=255)

                                                                                                                                              aCCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCAAGc      <  1:674664/150‑1 (MQ=255)

                                                                                                                                               ccAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCAAGCt     <  2:954959/150‑1 (MQ=255)

                                                                                                                                                 agcagGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTACTCCTATCTGGCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCAAGCTcg   >  2:1807975/1‑150 (MQ=255)

                                                                                                                                                  gcagGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCAAGCTcgc  >  1:1865131/1‑150 (MQ=255)

CAGCTCGATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCCATGGTAGGTGACCTGGCCGCCGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTGCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCAAGCTCGC  >  NC_002947/5464845‑5465140

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