BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	727,318	C→T	8.2%	W6* (TGG→TAG) ‡	PP_0621 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	727,318	0	C	T	8.2%	273.0 / 12.3	98	W6* (TGG→TAG) ‡	PP_0621	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (45/45); new base T (2/6); total (47/51)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

Only 100 of 103 total aligned reads displayed.

GGTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGG  >  NC_002947/727169‑727464

ggTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGccc                                                                                                                                                    <  1:618635/150‑1 (MQ=255)

ggTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGccc                                                                                                                                                    <  2:524944/150‑1 (MQ=255)

 gTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCGAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGcaca                                                                                                                                                   >  1:170795/1‑150 (MQ=255)

     gACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGc                                                                                                                                               <  2:1438952/150‑1 (MQ=255)

        tgctgGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTc                                                                                                                                            >  2:1873523/1‑150 (MQ=255)

          ctggctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTccc                                                                                                                                          <  2:1731688/150‑1 (MQ=255)

           tggctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGAACTCTGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCt                                                                                                                                         <  1:776047/150‑1 (MQ=255)

             gctgACTGTGTATGTGCGCAACAGAACGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCttt                                                                                                                                       <  2:2135399/150‑1 (MQ=255)

                gACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCtttttt                                                                                                                                    <  1:1710702/150‑4 (MQ=255)

                 aCTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCtttttta                                                                                                                                   <  1:2251234/150‑5 (MQ=255)

                   tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttcattca                                                                                                                                 >  1:1606108/1‑150 (MQ=255)

                   tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttcattca                                                                                                                                 >  1:218569/1‑150 (MQ=255)

                   tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttcattca                                                                                                                                 <  1:1866431/150‑1 (MQ=255)

                    ttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCAc                                                                                                                                <  2:218569/150‑1 (MQ=255)

                     tGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCAcc                                                                                                                               <  1:265689/150‑1 (MQ=255)

                         tttGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTg                                                                                                                           <  2:790321/150‑1 (MQ=255)

                            gcgcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAg                                                                                                                        <  1:2190252/150‑1 (MQ=255)

                              gcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAGGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGt                                                                                                                      <  1:1647636/150‑1 (MQ=255)

                              gcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGt                                                                                                                      >  2:1964661/1‑150 (MQ=255)

                              gcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGt                                                                                                                      >  1:1789214/1‑150 (MQ=255)

                                cacaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCt                                                                                                                                         >  1:1390534/1‑129 (MQ=255)

                                cacaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCt                                                                                                                                         <  2:1390534/129‑1 (MQ=255)

                                cacaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                    >  2:2067173/1‑150 (MQ=255)

                                 acaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAg                                                                                                                   >  1:1887620/1‑150 (MQ=255)

                                     ccAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttca                                                                                                                                     >  1:157151/1‑128 (MQ=255)

                                          aCGCGGCGGCGATCAGGGCGATCCAGAAAAGGAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGa                                                                                                          <  1:75659/150‑1 (MQ=255)

                                               gcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCaaa                                                                                                     <  2:1380407/150‑1 (MQ=255)

                                                cggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAg                                                                                                    >  1:1771238/1‑150 (MQ=255)

                                                cggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAg                                                                                                    <  2:1235348/150‑1 (MQ=255)

                                                   cgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTc                                                                                                 <  2:1606108/150‑1 (MQ=255)

                                                    gATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAGGGTCAGTCGAACAGACCAAAGATcg                                                                                                >  1:488671/1‑149 (MQ=255)

                                                     aTCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGCACAAACCaaaaaacac                                                                                               <  2:1203199/150‑7 (MQ=255)

                                                     aTCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCAt                                                                                               <  2:1789214/150‑1 (MQ=255)

                                                      tCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAACGGGAGGCCCAAGGCCTCCCttt                                                                                                                                       >  2:1851209/1‑109 (MQ=255)

                                                      tCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCttt                                                                                                                                       <  1:1851209/109‑1 (MQ=255)

                                                      tCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                    <  1:2113129/128‑1 (MQ=255)

                                                      tCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                    >  2:2113129/1‑128 (MQ=255)

                                                       cAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGt                                                                                             <  2:2284368/150‑1 (MQ=255)

                                                         gggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGTCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAg                                                                                                                   <  1:623097/126‑1 (MQ=255)

                                                         gggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGTCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAg                                                                                                                   >  2:623097/1‑126 (MQ=255)

                                                         gggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATCCACGGCTGAACGTCAGTCGAACAGACCAAAGGTCATGTAg                                                                                           >  1:2065320/1‑150 (MQ=255)

                                                          ggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTc                                                                                                                 <  2:1184849/127‑1 (MQ=255)

                                                          ggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTc                                                                                                                 >  1:1184849/1‑127 (MQ=255)

                                                          ggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGc                                                                                          <  1:72244/150‑1 (MQ=255)

                                                          ggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGACGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGc                                                                                          >  1:2221847/1‑150 (MQ=255)

                                                               tCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAAc                                                                                     >  1:1764173/1‑150 (MQ=255)

                                                               tCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAAc                                                                                     <  2:1795001/150‑1 (MQ=255)

                                                                 cAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCa                                                                                   >  2:116658/1‑150 (MQ=255)

                                                                  aGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCAAACAGACCAAAGGTCATGTAGCTGAACCAg                                                                                  >  1:361365/1‑150 (MQ=255)

                                                                   gAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAgg                                                                                 >  2:863448/1‑150 (MQ=255)

                                                                    aaaaGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGa                                                                                <  1:116658/150‑1 (MQ=255)

                                                                     aaaGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGaa                                                                               <  2:2221847/150‑1 (MQ=255)

                                                                           ggCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCt                                                                         <  2:1764173/150‑1 (MQ=255)

                                                                                           cAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGt                                                                                                                  <  2:1328957/93‑1 (MQ=255)

                                                                                           cAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGt                                                                                                                  >  1:1328957/1‑93 (MQ=255)

                                                                                            aGCTTGTCGATTGAAGACGCCAACGGTCTAGGAGGTGGCGATTAATAAATGGGGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                        <  2:1424202/150‑1 (MQ=255)

                                                                                            aGCGTGTCGATTCAAGAAGCCAACTGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAACGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGCCTTTCTCGGCTTCTTTCgg                                                        <  2:488671/150‑1 (MQ=255)

                                                                                             gCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCggg                                                       >  1:16416/1‑150 (MQ=255)

                                                                                              cTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCg                                                         <  2:1279273/147‑1 (MQ=255)

                                                                                              cTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCg                                                         >  1:1279273/1‑147 (MQ=255)

                                                                                                 gTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGa                                                                                                          <  2:2137193/95‑1 (MQ=255)

                                                                                                 gTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGa                                                                                                          >  1:2137193/1‑95 (MQ=255)

                                                                                                 gTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCtt                                                   >  2:1135238/1‑150 (MQ=255)

                                                                                                  tCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTc                                                  <  1:1135238/150‑1 (MQ=255)

                                                                                                  tCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTc                                                  <  2:903715/150‑1 (MQ=255)

                                                                                                   cGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTcc                                                 <  1:2067173/150‑1 (MQ=255)

                                                                                                     aTTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCgg                                               >  2:1492809/1‑150 (MQ=255)

                                                                                                         aagaagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACg                                           >  1:783086/1‑150 (MQ=255)

                                                                                                          agaagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGa                                          <  2:1887620/150‑1 (MQ=255)

                                                                                                           gaagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAt                                         >  2:454033/1‑150 (MQ=255)

                                                                                                            aagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGc                                                                 <  1:1739341/125‑1 (MQ=255)

                                                                                                            aagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGc                                                                 >  2:1739341/1‑125 (MQ=255)

                                                                                                            aagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAtt                                        >  2:2002150/1‑150 (MQ=255)

                                                                                                            aagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTACGGCACGAtt                                        >  1:1645625/1‑150 (MQ=255)

                                                                                                               cgaaCAGTCTAGGGCGCGGCCAATAAAAAAAGGAGGCCCAAGGCTTCCCTTTCATGCACCTCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTggg                                     <  2:170795/148‑1 (MQ=255)

                                                                                                               ccAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTggg                                     <  1:863448/150‑1 (MQ=255)

                                                                                                                 aaCAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTc                                   >  1:303684/1‑150 (MQ=255)

                                                                                                                    aGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGcc                                >  2:2241076/1‑150 (MQ=255)

                                                                                                                      tCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAg                                                                                           >  1:582428/1‑89 (MQ=255)

                                                                                                                      tCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAg                                                                                           <  2:582428/89‑1 (MQ=255)

                                                                                                                         aGGACGCGGCCAATAACAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTtct                           >  1:361165/1‑150 (MQ=255)

                                                                                                                          ggACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTtctc                          >  1:655853/1‑150 (MQ=255)

                                                                                                                              gcgGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCGCGTCt                      <  2:1645625/150‑1 (MQ=255)

                                                                                                                                ggCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                        >  1:1577259/1‑114 (MQ=255)

                                                                                                                                ggCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                        <  2:1577259/114‑1 (MQ=255)

                                                                                                                                  ccAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGct                                                                <  1:1052041/104‑1 (MQ=255)

                                                                                                                                  ccAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGct                                                                >  2:1052041/1‑104 (MQ=255)

                                                                                                                                     aTAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgca               >  1:237281/1‑150 (MQ=255)

                                                                                                                                        aaaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTtctc                          <  2:1450552/136‑1 (MQ=255)

                                                                                                                                        aaaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTtctc                          >  1:1450552/1‑136 (MQ=255)

                                                                                                                                        aaaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcagca            <  2:2065320/150‑1 (MQ=255)

                                                                                                                                         gaaaGGGAGGCCTTGGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCt                      <  2:1394325/138‑1 (MQ=255)

                                                                                                                                         gaaaGGGAGGCCTTGGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCt                      >  1:1394325/2‑139 (MQ=255)

                                                                                                                                         aaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTTTTGAGGcggcag           >  2:1299250/1‑150 (MQ=255)

                                                                                                                                         aaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGcagcag           >  1:509003/1‑150 (MQ=255)

                                                                                                                                              ggAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcag              <  1:752839/142‑1 (MQ=255)

                                                                                                                                              ggAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcag              >  2:752839/1‑142 (MQ=255)

                                                                                                                                                  gCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCgg  >  2:1169571/1‑150 (MQ=255)

                                                                                                                                                    ccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTc                                                  <  2:343733/100‑1 (MQ=255)

                                                                                                                                                    ccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTc                                                  >  1:343733/1‑100 (MQ=255)

GGTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGG  >  NC_002947/727169‑727464

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