BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,395,567	G→C	5.3%	A84G (GCC→GGC)	PP_3871 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,395,567	0	G	C	5.3%	337.6 / 7.5	107	A84G (GCC→GGC)	PP_3871	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (49/51); new base C (3/3); total (52/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 111 total aligned reads displayed.

GAGAAGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAG  >  NC_002947/4395423‑4395715

gagaAGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGc‑c                                                                                                                                                 <  2:2612867/150‑1 (MQ=255)

    aGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGCCGCCGAC‑ctgct                                                                                                                                             <  2:195106/150‑1 (MQ=255)

      cTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTgg                                                                                                                                           <  2:2924884/150‑1 (MQ=255)

       tGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGc                                                                                                                                          >  2:1826685/1‑150 (MQ=255)

        gAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGcc                                                                                                                                         >  2:834886/1‑150 (MQ=255)

        gAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGcc                                                                                                                                         >  2:1091987/1‑150 (MQ=255)

        gAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGcc                                                                                                                                         >  1:986613/1‑150 (MQ=255)

         aaGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGAAGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCcg                                                                                                                                        >  1:848094/1‑150 (MQ=255)

           gCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCGGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCcgcg                                                                                                                                      >  2:2669209/1‑150 (MQ=255)

            cGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACTATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCgcgc                                                                                                                                     <  1:2928448/150‑1 (MQ=255)

               tAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCg                                                                                                                                  >  2:2663714/1‑150 (MQ=255)

                     accCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTg                                                                                                                            >  2:1547185/1‑150 (MQ=255)

                       ccGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCgcgc                                                                                                                                     <  1:2023721/139‑1 (MQ=255)

                       ccGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGgt                                                                                                                          >  2:2797931/1‑150 (MQ=255)

                       ccGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGgt                                                                                                                          >  1:1672869/1‑150 (MQ=255)

                       ccGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCGGGCCgcgc                                                                                                                                     >  2:2023721/1‑139 (MQ=255)

                       ccGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCAGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGgt                                                                                                                          >  1:985768/1‑150 (MQ=255)

                         gcgcACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGgtgt                                                                                                                        <  2:848094/150‑1 (MQ=255)

                                   gATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGaa                                                                                                              <  2:986613/150‑1 (MQ=255)

                                    ttcaacaCCAAACCGCCACCGTCGTAGACGATTGGGTCGTACTCATCACCAGTCCGCCATCCGACCAAGGCCCGCCACAAGA‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CCGCTGCCCGCGCTCGTCCTTGGTGTCCATGGCGAac                                                                                                             <  1:196737/144‑1 (MQ=255)

                                    aTCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAac                                                                                                             <  1:2797931/150‑1 (MQ=255)

                                     tCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAaca                                                                                                            >  1:1335448/1‑150 (MQ=255)

                                     tCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAaca                                                                                                            >  2:3010150/1‑150 (MQ=255)

                                     tCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGCCGCCGAC‑ctgct                                                                                                                                             >  2:2137146/1‑117 (MQ=255)

                                     tCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGCCGCCGAC‑ctgct                                                                                                                                             <  1:2137146/117‑1 (MQ=255)

                                      cAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAacac                                                                                                           <  2:1420460/150‑1 (MQ=255)

                                        gTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGgt                                                                                                                          <  2:811236/133‑1 (MQ=255)

                                        gTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGgt                                                                                                                          >  1:811236/1‑133 (MQ=255)

                                         tACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAAcaccac                                                                                                        <  1:2669209/150‑1 (MQ=255)

                                         tACCAACCCGCCACCGTCGTAGACGATGGGGTCGCACTCATCACCAGGCCGCCATCCGCCAAACGCCCGCCACCACTGACCATTC‑ATGGTTG‑ACATGTGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTCTCCAGTGCGAActccac                                                                                                        >  1:2228534/1‑150 (MQ=255)

                                           ccAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTg                                                                                                                            >  2:2452078/1‑128 (MQ=255)

                                           ccAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTg                                                                                                                            <  1:2452078/128‑1 (MQ=255)

                                            cAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCAGGGCGAACACCACCTc                                                                                                     >  1:114780/1‑150 (MQ=255)

                                            cAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGCCGCCGAC‑ctgctgcac                                                                                                                                         <  2:909697/114‑4 (MQ=255)

                                            cAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGCCGCCGAC‑ctgctgcac                                                                                                                                         >  1:909697/1‑111 (MQ=255)

                                               cccGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGCCCTTGGTGTCCATGGCGAACACCACCTCaaa                                                                                                  >  2:984968/1‑150 (MQ=255)

                                               cccGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCAGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCaaa                                                                                                  <  2:985768/150‑1 (MQ=255)

                                                      ccGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCAt                                                                                                                    <  2:138640/125‑1 (MQ=255)

                                                      ccGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCAt                                                                                                                    >  1:138640/1‑125 (MQ=255)

                                                            tAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGCCCCGCCACAACT‑‑CCGTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATg                                                                                     <  1:1437887/150‑1 (MQ=255)

                                                                 gATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTg                                                                                <  1:2166375/150‑1 (MQ=255)

                                                                 gATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTg                                                                                >  2:1626585/1‑150 (MQ=255)

                                                                    ggggTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGt                                                                             >  1:1642030/1‑150 (MQ=255)

                                                                       gTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGt                                                                             <  2:1848544/147‑1 (MQ=255)

                                                                       gTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGt                                                                             >  1:1848544/1‑147 (MQ=255)

                                                                            aCTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCtt                                                                                                                             >  2:2604624/1‑94 (MQ=255)

                                                                             cTCATGACCAGGCCGCCATCCGGTTATTGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTTTCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTTCTGCGTAACTGCCtt                                                                    <  2:1256171/150‑1 (MQ=255)

                                                                              tcatcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTg                                                                                >  2:2795002/1‑137 (MQ=255)

                                                                              tcatcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTg                                                                                <  1:2795002/137‑1 (MQ=255)

                                                                              tcatcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTc                                                                   <  1:1091987/150‑1 (MQ=255)

                                                                               catcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCt                                                                  >  1:510867/1‑150 (MQ=255)

                                                                                atcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGa                                                                                       <  1:2789893/128‑1 (MQ=255)

                                                                                atcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGa                                                                                       >  2:2789893/1‑128 (MQ=255)

                                                                                atcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTggg                                                                              >  1:1634679/1‑137 (MQ=255)

                                                                                atcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTggg                                                                              <  2:1634679/137‑1 (MQ=255)

                                                                                atcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTg                                                                 >  1:627630/1‑150 (MQ=255)

                                                                                  caCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTgcg                                                               >  2:2960711/1‑150 (MQ=255)

                                                                                   aCCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTgcgc                                                              >  1:1827846/1‑150 (MQ=255)

                                                                                     cAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAAcaccac                                                                                                        >  1:115230/1‑106 (MQ=255)

                                                                                     cAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAAcaccac                                                                                                        <  2:115230/106‑1 (MQ=255)

                                                                                      cggacgccATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGT‑CAGCCGCGAGTTTTCGCCGCACTCGGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCtt                                                           <  1:2663714/146‑1 (MQ=255)

                                                                                       ggCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTg                                                          <  2:1642030/150‑1 (MQ=255)

                                                                                         ccgccATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGtt                                                        <  2:1356435/150‑1 (MQ=255)

                                                                                          cgccATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTtc                                                       >  2:591963/1‑150 (MQ=255)

                                                                                              aTCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTgc                                                                >  1:41144/1‑137 (MQ=255)

                                                                                              aTCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTgc                                                                <  2:41144/137‑1 (MQ=255)

                                                                                               tCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGcc                                                  >  1:733910/1‑150 (MQ=255)

                                                                                                ccGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAAcaccac                                                                                                        >  2:1503488/1‑95 (MQ=255)

                                                                                                ccGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAAcaccac                                                                                                        <  1:1503488/95‑1 (MQ=255)

                                                                                                 cGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGc                                                                       <  1:1970139/127‑1 (MQ=255)

                                                                                                 cGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGc                                                                       >  2:1969669/1‑127 (MQ=255)

                                                                                                 cGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGc                                                                       >  2:1970139/1‑127 (MQ=255)

                                                                                                  gACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGc                                               >  1:573644/1‑150 (MQ=255)

                                                                                                    ccAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCtt                                             <  1:3010150/150‑1 (MQ=255)

                                                                                                     cAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCt                                                                                          <  1:1199372/104‑1 (MQ=255)

                                                                                                     cAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCt                                                                                          >  2:1199372/1‑104 (MQ=255)

                                                                                                         gCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGgtg                                                                                                                         >  2:1161014/1‑69 (MQ=255)

                                                                                                         gCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGgtg                                                                                                                         <  1:1161014/69‑1 (MQ=255)

                                                                                                               cacaACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTgcgc                                                              >  1:2834270/1‑122 (MQ=255)

                                                                                                                acaACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCg                                                                      <  2:2033832/113‑1 (MQ=255)

                                                                                                                acaACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCg                                                                      >  1:2033832/1‑113 (MQ=255)

                                                                                                                 caACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGTTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGt                                <  2:1202208/150‑1 (MQ=255)

                                                                                                                 caACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGTGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGt                                <  2:358718/150‑1 (MQ=255)

                                                                                                                 caACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGt                                <  1:1826685/150‑1 (MQ=255)

                                                                                                                  aaCT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTg                               <  2:1335448/150‑1 (MQ=255)

                                                                                                                  aaCT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTg                               >  1:2180011/1‑150 (MQ=255)

                                                                                                                   aCT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGa                              <  2:733910/150‑1 (MQ=255)

                                                                                                                    cT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGAt                             <  1:1626585/150‑1 (MQ=255)

                                                                                                                     t‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGTTCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATc                            <  2:696230/150‑1 (MQ=255)

                                                                                                                         cTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATcac                          <  2:1827846/150‑1 (MQ=255)

                                                                                                                          ttGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATcaca                         >  2:484728/1‑150 (MQ=255)

                                                                                                                          ttGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATcaca                         >  1:2136680/1‑150 (MQ=255)

                                                                                                                              gAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAgg                     <  2:573644/150‑1 (MQ=255)

                                                                                                                                  gTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATcaca                         >  1:2743514/1‑142 (MQ=255)

                                                                                                                                  gTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATcaca                         <  2:2743514/142‑1 (MQ=255)

                                                                                                                                        cagcagGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTg                                                                                <  2:604846/81‑1 (MQ=255)

                                                                                                                                        cagcagGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTg                                                                                >  1:604846/1‑81 (MQ=255)

                                                                                                                                                gcgcggc‑CAGCAGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGAt                             <  2:1659580/122‑1 (MQ=255)

                                                                                                                                                gcgcggc‑CAGCAGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGAt                             >  1:1659580/3‑124 (MQ=255)

                                                                                                                                                 ggcggc‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAg  >  1:1330108/1‑150 (MQ=255)

GAGAAGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACT‑‑CCTTGCGAAAGTCGTGCAGCAGGTCGGCGGC‑CTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAG  >  NC_002947/4395423‑4395715

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