BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,741,946	C→T	5.1%	E88E (GAG→GAA)	PP_5039 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,741,946	0	C	T	5.1%	292.3 / 3.2	97	E88E (GAG→GAA)	PP_5039	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (44/48); new base T (2/3); total (46/51)
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.90e-01

Only 100 of 101 total aligned reads displayed.

GATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGC  >  NC_002947/5741798‑5742095

gATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCt                                                                                                                                                      >  1:515429/1‑150 (MQ=255)

 aTAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCtctc                                                                                                                                                     >  1:1103625/1‑150 (MQ=255)

    aCGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGaa                                                                                                                                                  <  2:1787830/150‑1 (MQ=255)

     cGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCgcaa                                                                                                                                                 <  1:1664782/150‑4 (MQ=255)

     cGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                 <  2:1232743/150‑1 (MQ=255)

     cGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                 <  1:1444970/150‑1 (MQ=255)

      gCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCGGCAGCACGAACAGGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGCCCACCTAGGCCCGCCAGCCGCCCTCGGAgg                                                                                                                                                >  2:1321800/1‑150 (MQ=255)

        tcgaCGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCATCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcg                                                                                                                                              >  1:550785/4‑150 (MQ=255)

        tCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcg                                                                                                                                              >  2:1130811/1‑150 (MQ=255)

        cccACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGACCTGCATGGCGACCCCCTGGGCACGCCGGTCTCCCTCGAAGgcg                                                                                                                                              <  2:1728411/149‑1 (MQ=255)

          cACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGATGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcgcg                                                                                                                                            >  1:1000639/1‑150 (MQ=255)

           aCGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGACACGCCGCGCGCCTTCGAGGgcgacc                                                                                                                                           >  2:633628/1‑147 (MQ=255)

              cacaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGACACGCCGCGCGCCTTCGAGGgcgaccggc                                                                                                                                        <  1:841074/150‑7 (MQ=255)

              cacaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGACACGCCGCGCGCCTTCGAGGgcgaccggc                                                                                                                                        <  1:2343799/150‑7 (MQ=255)

               acaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGACCGGGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCAAGCCGGTCGCCCTCGAAGGCGCgcggcg                                                                                                                                       >  1:1108058/1‑150 (MQ=255)

                caGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCgt                                                                                                                                      <  2:1000639/150‑1 (MQ=255)

                    ttCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGACACGCCGCGCGCCTTCGAGGGCGACCGGCgtgccc                                                                                                                                  <  1:633628/150‑4 (MQ=38)

                       tGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                                               >  2:2292150/1‑150 (MQ=255)

                        gccaggccagCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGt                                                                                                                              <  2:131368/150‑1 (MQ=255)

                        gCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGt                                                                                                                              <  1:1274651/150‑1 (MQ=255)

                         ccAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTa                                                                                                                             <  1:1527657/150‑1 (MQ=255)

                         ccAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTa                                                                                                                             >  2:1350219/1‑150 (MQ=255)

                          cAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGt                                                                                                                              <  1:1341836/148‑1 (MQ=255)

                          cAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGt                                                                                                                              >  2:1341836/1‑148 (MQ=255)

                           aGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGa                                                                                                                           >  1:446785/1‑150 (MQ=255)

                               cAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                                               <  2:2039517/142‑1 (MQ=255)

                               cAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                                               >  1:2039517/1‑142 (MQ=255)

                                 gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGc                                                                                                                     >  1:77016/1‑150 (MQ=255)

                                   aCTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGCCGCGCGGCGTGTCGCGGTAGATCATGCGt                                                                                                                   >  2:180315/1‑150 (MQ=255)

                                   aCTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTCCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCCCCCCGCGTGTGGGGGTAGGGCATGCGt                                                                                                                   >  2:180336/1‑150 (MQ=255)

                                    cTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCgt                                                                                                                                      <  1:1928844/130‑1 (MQ=255)

                                    cTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCgt                                                                                                                                      >  2:1928844/1‑130 (MQ=255)

                                     ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTgg                                                                                                                 >  2:1828237/1‑150 (MQ=255)

                                          cATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGt                                                                                                                   >  2:2343832/1‑143 (MQ=255)

                                          cATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGt                                                                                                                   <  1:2343832/143‑1 (MQ=255)

                                              cTCCAGTACTAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCCGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTAGCGGTAGTTCATGCGTCGCAGATGGTc                                                                                                        <  2:1446956/150‑1 (MQ=255)

                                               tCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCGGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCg                                                                                                       >  2:1594211/1‑150 (MQ=255)

                                                ccAGTACCAGCAGCCCGACCACGAGCTGGATCGACAGCGGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                      <  2:857279/150‑1 (MQ=255)

                                                ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                      <  1:1386592/150‑1 (MQ=255)

                                                ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                      <  1:1321800/150‑1 (MQ=255)

                                                ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCCGCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                      <  1:1208121/150‑1 (MQ=255)

                                                         gcagcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATca                                                                                             >  2:1141318/1‑150 (MQ=255)

                                                          cagcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcag                                                                                            <  2:828561/150‑1 (MQ=255)

                                                           agcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCTCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagc                                                                                           <  1:662540/150‑1 (MQ=255)

                                                             caCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagcag                                                                                         >  2:614310/1‑150 (MQ=255)

                                                              aCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAgg                                                                                        <  1:454406/150‑1 (MQ=255)

                                                                gAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGcc                                                                                      >  1:375446/1‑150 (MQ=255)

                                                                         cTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCagccagcc                                                                             >  2:102826/1‑150 (MQ=255)

                                                                            gATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTg                                                                          <  2:1467010/150‑1 (MQ=255)

                                                                              tCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCg                                                                        >  1:628250/1‑150 (MQ=255)

                                                                              tCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGGCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCg                                                                        >  2:527034/1‑150 (MQ=255)

                                                                               cGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgg                                                                       <  1:1350219/150‑1 (MQ=255)

                                                                                ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                      <  1:2078463/150‑1 (MQ=255)

                                                                                ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                      <  1:1951345/150‑1 (MQ=255)

                                                                                ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                      <  1:2004858/150‑1 (MQ=255)

                                                                                 gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                                     >  2:809133/1‑150 (MQ=255)

                                                                                 gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                      <  1:867465/149‑1 (MQ=255)

                                                                                      ttcTGGTAGTCCCAGCCCTTCTCCAGGGCGTTCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgc                                                                <  2:1096776/148‑1 (MQ=255)

                                                                                       gcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                 <  2:1116517/68‑1 (MQ=255)

                                                                                       gcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                 >  1:1116517/1‑68 (MQ=255)

                                                                                       gcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCccagc                                                                                  <  1:1553203/131‑1 (MQ=255)

                                                                                       gcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCccagc                                                                                  >  2:1553203/1‑131 (MQ=255)

                                                                                       gcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTAGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                                               <  2:417615/150‑1 (MQ=255)

                                                                                           gTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgctagat                                                           >  1:1096776/1‑146 (MQ=255)

                                                                                               tCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCActct                                                       >  1:1231460/1‑150 (MQ=255)

                                                                                                 ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                                     >  2:1850905/1‑150 (MQ=255)

                                                                                                  cAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCg                                                    >  1:2184823/1‑150 (MQ=255)

                                                                                                  cAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCg                                                    >  2:1574394/1‑150 (MQ=255)

                                                                                                   aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGa                                                   <  1:809133/150‑1 (MQ=255)

                                                                                                       cTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagca                                                                                          <  1:278433/107‑1 (MQ=255)

                                                                                                       cTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagca                                                                                          >  2:278433/1‑107 (MQ=255)

                                                                                                       cTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGGCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCCGCGGGACTTGCTGCAACTCTAGGAGcgc                                               >  2:1091963/1‑148 (MQ=255)

                                                                                                           tGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCAccc                                           <  1:9965/150‑1 (MQ=255)

                                                                                                           tGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCAccc                                           >  1:2116003/1‑150 (MQ=255)

                                                                                                             cAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                                     >  1:1086259/1‑122 (MQ=255)

                                                                                                             cAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                                     <  2:1086259/122‑1 (MQ=255)

                                                                                                               ggCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                      <  1:142694/87‑1 (MQ=255)

                                                                                                               ggCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                      >  2:142694/1‑87 (MQ=255)

                                                                                                                 cccTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAgg                                     <  1:1574394/150‑1 (MQ=255)

                                                                                                                   cTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGc                                   <  1:2407715/150‑1 (MQ=255)

                                                                                                                   cTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGc                                   <  2:1231460/150‑1 (MQ=255)

                                                                                                                     gCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                      >  1:1174827/1‑113 (MQ=255)

                                                                                                                     gCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                      <  2:1174827/113‑1 (MQ=255)

                                                                                                                       aTGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTGGTTCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCATGTTTCCCGTGGGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGc                               <  1:1636885/150‑1 (MQ=255)

                                                                                                                         ggCGACCACCTGGGCACGCCGGCCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTc                             <  1:1531851/150‑1 (MQ=255)

                                                                                                                            gACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGAt                                                                                                                          <  1:673730/54‑1 (MQ=255)

                                                                                                                            gACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGAt                                                                                                                          >  2:673730/1‑54 (MQ=255)

                                                                                                                            gACCACCTGGGCACGCCGGTCGCACTCGAAGGCGCGGGGGGTGTGGCGGTAGATCATGCGTGGCAGCTGGTGGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGAc                          <  2:1103625/150‑1 (MQ=255)

                                                                                                                              ccaccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCa                        <  2:2184823/150‑1 (MQ=255)

                                                                                                                                accTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAgg                      >  2:742815/1‑150 (MQ=255)

                                                                                                                                accTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAgg                      >  2:489335/1‑150 (MQ=255)

                                                                                                                                  cTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTc                    >  1:699944/1‑150 (MQ=255)

                                                                                                                                   tGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCg                   <  1:614310/150‑1 (MQ=255)

                                                                                                                                      acACGCCGCGCGCCTTCGAGGGCGACCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGGTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTCCTCGACCAGGTCGCCa                >  2:135457/2‑150 (MQ=255)

                                                                                                                                            cGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGTGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAgc          <  1:1594211/150‑1 (MQ=255)

                                                                                                                                            cGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAgc          <  1:705939/150‑1 (MQ=255)

                                                                                                                                                 gCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATg     >  2:1083705/1‑150 (MQ=255)

                                                                                                                                                  cccTCGAAGGCGCGCGGCGTGGCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATgg    <  2:550785/150‑1 (MQ=255)

                                                                                                                                                    cTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGACAGAGCTTGCTCGACCAGGTCGCCAAACAGCGCATGGGc  >  2:1145413/1‑150 (MQ=255)

                                                                                                                                                    cTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCATGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGc  >  2:444082/1‑150 (MQ=255)

GATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGC  >  NC_002947/5741798‑5742095

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