BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,184,165	C→T	6.9%	A85V (GCC→GTC)	PP_4564 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,184,165	0	C	T	6.9%	338.4 / 11.1	116	A85V (GCC→GTC)	PP_4564	lipoprotein
Reads supporting (aligned to +/- strand): ref base C (54/54); new base T (2/6); total (56/60)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 118 total aligned reads displayed.

GCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGG  >  NC_002947/5184019‑5184311

gCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCcga                                                                                                                                                 <  2:1725685/150‑1 (MQ=255)

gCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCcga                                                                                                                                                 >  1:1107520/1‑150 (MQ=255)

  cGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCcgacg                                                                                                                                               >  1:1244057/1‑150 (MQ=255)

   gTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCcgacga                                                                                                                                              <  1:1461307/150‑1 (MQ=255)

   gTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCcgacga                                                                                                                                              >  1:986670/1‑150 (MQ=255)

    tGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCgacgac                                                                                                                                             >  2:155186/1‑150 (MQ=255)

    tGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCgacgac                                                                                                                                             >  2:636866/1‑150 (MQ=255)

      gTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCa                                                                                                                                           >  2:1122102/1‑150 (MQ=255)

       tCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAg                                                                                                                                          <  2:1244057/150‑1 (MQ=255)

       tCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGAACGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAg                                                                                                                                          <  1:1218999/150‑1 (MQ=255)

         gatGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                        >  1:1348253/3‑150 (MQ=255)

         gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCTTGGCATTGCAAGACTCAACTGCCCGGCTACCGCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                        <  2:153955/150‑1 (MQ=255)

         gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGCCCGGCGACGTTACATGGTCGCCTACGACCAggg                                                                                                                                        >  2:1361560/1‑149 (MQ=255)

         gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGt                                                                                                                                        <  2:1061433/150‑1 (MQ=255)

         gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                        <  1:1271014/150‑1 (MQ=255)

           tGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCa                                                                                                                                      >  2:1566369/1‑150 (MQ=255)

            gAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCAc                                                                                                                                     <  2:901903/150‑1 (MQ=255)

              gCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCg                                                                                                                                   >  2:1812749/1‑150 (MQ=255)

                      cTGGACATGACCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGcc                                                                                                                           <  1:770702/150‑1 (MQ=255)

                      cTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGcc                                                                                                                           >  1:1209856/1‑150 (MQ=255)

                      cTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGcc                                                                                                                           >  1:651761/1‑150 (MQ=255)

                       tGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGccc                                                                                                                          >  1:1015059/1‑150 (MQ=255)

                           caggAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTgca                                                                                                                      <  2:1180170/147‑1 (MQ=255)

                           cATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTgca                                                                                                                      >  2:504189/1‑150 (MQ=255)

                            aTGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGATCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTgcag                                                                                                                     >  2:1156034/1‑150 (MQ=255)

                                aCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                                 <  1:804620/150‑1 (MQ=255)

                                 ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAg                                                                                                                <  2:1348253/150‑1 (MQ=255)

                                 ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAg                                                                                                                >  2:1914070/1‑150 (MQ=255)

                                  cgcgAGCAAAGGGCGCTGCTGGATGATCGTTAGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCat                                                                                                                                          <  1:498601/123‑2 (MQ=255)

                                  cgcgAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCat                                                                                                                                          >  2:498601/1‑122 (MQ=255)

                                   gcgAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCa                                                                                                              >  1:1633108/1‑150 (MQ=255)

                                     gAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCAcc                                                                                                            <  1:1122102/150‑1 (MQ=255)

                                               cgcTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTcc                                                                                                  <  1:522496/150‑1 (MQ=255)

                                                 ctgctgGATGATCGGTTGCGCGCGCACCTGGCCTGGCCTTCTACGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGt                                                                                                <  1:679848/150‑1 (MQ=255)

                                                 ctgctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGt                                                                                                <  1:1385615/150‑1 (MQ=255)

                                                  tgctgGATGATCGGTTGCGCGAGAACGTTGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTg                                                                                               >  1:1697525/1‑150 (MQ=255)

                                                   gctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGATCGGCTACGGGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGa                                                                                              <  1:1156034/150‑1 (MQ=255)

                                                   gctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGa                                                                                              <  1:1742468/150‑1 (MQ=255)

                                                     tgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAgg                                                                                            <  1:545318/150‑1 (MQ=255)

                                                     tgGATGAGCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAgg                                                                                            <  2:1209856/150‑1 (MQ=255)

                                                      ggATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGc                                                                                           >  2:1500779/1‑150 (MQ=255)

                                                           atCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTg                                                                                      <  1:991264/150‑1 (MQ=255)

                                                            tCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATCgg                                                                                     >  2:244199/1‑150 (MQ=255)

                                                             cGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGTAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGt                                                                                    <  2:1015059/150‑1 (MQ=255)

                                                              ggTTGCGCGAGCACCTGGCCTGGGATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTc                                                                                   >  2:949107/1‑150 (MQ=255)

                                                               gTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCg                                                                                  >  2:422901/1‑150 (MQ=255)

                                                               gTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGATCg                                                                                  >  1:27098/1‑150 (MQ=255)

                                                                  gcgcGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGACTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                               >  2:175175/1‑150 (MQ=255)

                                                                        gCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTAACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGa                                                                         >  1:1584892/1‑150 (MQ=255)

                                                                         cACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTAACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAg                                                                        <  2:1584892/150‑1 (MQ=255)

                                                                          aCCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAgg                                                                       >  2:1168349/1‑150 (MQ=255)

                                                                               gcctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCAGCGGAAATc                                                                  <  2:27098/150‑1 (MQ=255)

                                                                                cctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCCAGGAAATCa                                                                 >  2:767166/1‑150 (MQ=255)

                                                                                 ctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAc                                                                >  1:1697135/1‑150 (MQ=255)

                                                                                 ctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAc                                                                >  1:176787/1‑150 (MQ=255)

                                                                                   ggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGACCGTCAGGCAATCGGTCGATTGGCCGAGGAAATCAccc                                                              >  2:627666/1‑150 (MQ=255)

                                                                                          cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGc                                                       >  1:1098236/1‑150 (MQ=255)

                                                                                          cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGc                                                       <  1:1914070/150‑1 (MQ=255)

                                                                                           aaGACTCAACTGCCCGGCTACCTCGACTTGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGcc                                                      >  2:351458/1‑150 (MQ=255)

                                                                                           aaGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCa                                                                                                                               >  2:1885496/1‑77 (MQ=255)

                                                                                           aaGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCa                                                                                                                               <  1:1885496/77‑1 (MQ=255)

                                                                                               cTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCg                                                  >  1:1320245/1‑150 (MQ=255)

                                                                                                  aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCg                                                         <  2:620154/140‑1 (MQ=255)

                                                                                                  aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCg                                                         >  1:620154/1‑140 (MQ=255)

                                                                                                  aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAAc                                               >  2:859053/1‑150 (MQ=255)

                                                                                                      gCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                                 >  1:461947/1‑80 (MQ=255)

                                                                                                      gCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                                 <  2:461947/80‑1 (MQ=255)

                                                                                                          ggCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCg                                                         >  1:674850/1‑132 (MQ=255)

                                                                                                          ggCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTgcgc                                       <  1:859053/150‑1 (MQ=255)

                                                                                                                cTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATg                                 >  1:1557871/1‑150 (MQ=255)

                                                                                                                 tCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGt                                                                                    <  1:1150486/98‑1 (MQ=255)

                                                                                                                 tCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGt                                                                                    >  2:1150486/1‑98 (MQ=255)

                                                                                                                 tCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGa                                >  2:1726655/1‑150 (MQ=255)

                                                                                                                 tCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGa                                >  1:1184178/1‑150 (MQ=255)

                                                                                                                 tCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGa                                <  2:1107520/150‑1 (MQ=255)

                                                                                                                  cGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCGCAGGCCACCGAACTGCTGCGCGGCATGAg                               <  1:989194/150‑1 (MQ=255)

                                                                                                                   gACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGc                              <  2:1633108/150‑1 (MQ=255)

                                                                                                                   gACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGc                              <  1:1582360/150‑1 (MQ=255)

                                                                                                                     cTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGa                            <  1:627666/150‑1 (MQ=255)

                                                                                                                        gCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACaa                         <  2:1697135/150‑1 (MQ=255)

                                                                                                                         cTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTc                                                                                   <  2:671488/91‑1 (MQ=255)

                                                                                                                         cTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTc                                                                                   >  1:671488/1‑91 (MQ=255)

                                                                                                                          tCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGaa                                                >  1:585627/1‑125 (MQ=255)

                                                                                                                          tCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGaa                                                <  2:585627/125‑1 (MQ=255)

                                                                                                                           cGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCa                      <  1:472314/150‑1 (MQ=255)

                                                                                                                               cGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCACGGCATGGACGACAACCAGGTa                  >  2:439112/1‑150 (MQ=255)

                                                                                                                                 gCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAActgct                                           >  1:1637843/1‑123 (MQ=255)

                                                                                                                                 gCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCACTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgc                >  1:1618371/1‑150 (MQ=255)

                                                                                                                                      gTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGa           <  2:1320245/150‑1 (MQ=255)

                                                                                                                                         aCATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGc        >  2:1319387/1‑150 (MQ=255)

                                                                                                                                         aCATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGc        <  1:410919/150‑1 (MQ=255)

                                                                                                                                         aCATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGAAGACGGGCCGAGGAGAAGACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGc        <  1:1178716/150‑1 (MQ=255)

                                                                                                                                           ctGGTCGTCGGCGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAAt                                                                   >  2:1910750/2‑89 (MQ=255)

                                                                                                                                           aTGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATNCGt      <  2:1194136/150‑1 (MQ=255)

                                                                                                                                           aTGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGt      <  2:12390/150‑1 (MQ=255)

                                                                                                                                             ggTCGCCGACGCACAGTTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGa    <  2:986670/150‑1 (MQ=255)

                                                                                                                                             ggTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGa    <  1:1361560/150‑1 (MQ=255)

                                                                                                                                               tCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAgg  >  2:1331195/1‑150 (MQ=255)

                                                                                                                                                  ccGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGa                                                                                 >  2:666883/1‑68 (MQ=255)

                                                                                                                                                  ccGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGa                                                                                 <  1:666883/68‑1 (MQ=255)

GCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGG  >  NC_002947/5184019‑5184311

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