BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,463,997	A→T	12.2%	*139C (TGA→TGT)	PP_3956 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,463,997	0	A	T	12.2%	318.0 / 31.4	115	*139C (TGA→TGT)	PP_3956	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (58/43); new base T (5/9); total (63/52)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 116 total aligned reads displayed.

CCCGAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCA  >  NC_002947/4463848‑4464143

cccGAGGCTGGTGACTACCAGGTCTGCACGACAGACCGGTGCCAGGAATTGGGTGATTGGGTAGGACCGCCGAGCTAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGa                                                                                                                                                    >  1:646708/1‑150 (MQ=255)

   gAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTGAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACgg                                                                                                                                                 >  2:416208/1‑150 (MQ=255)

          gTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGCAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAgg                                                                                                                                          <  1:1326338/150‑1 (MQ=255)

            gACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGc                                                                                                                                        <  2:470794/150‑1 (MQ=255)

            gACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGc                                                                                                                                        >  1:57777/1‑150 (MQ=255)

            gACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGc                                                                                                                                        <  2:2249454/150‑1 (MQ=255)

              cTACCAGGTCTGCACGACAGATGGGTGCCAGCCACTGTCTTATTGGTTAGGACCGCCGAGCGATCGCATTGCTGCAGTAGACTTCTGTGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGTCCTGCGCTGTCGCCGTCAGGGCAg                                                                                                                                      <  2:1209316/150‑1 (MQ=255)

              cTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGCCCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGCAg                                                                                                                                      <  2:1333976/150‑1 (MQ=255)

              cTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGCAg                                                                                                                                      <  2:1296411/150‑1 (MQ=255)

              cTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAg                                                                                                                                      <  2:2127842/150‑1 (MQ=255)

              cTACCAGGTCTGCACGACAGACCGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAg                                                                                                                                      <  2:1202152/150‑1 (MQ=255)

               tACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTATGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGCAg                                                                                                                                      <  1:1684099/149‑1 (MQ=255)

               tACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTATGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGCAg                                                                                                                                      >  2:1684099/1‑149 (MQ=255)

                 ccAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGTCCTGCCCTGTCGCCGTCAGGGCAgggc                                                                                                                                   <  1:832408/150‑3 (MQ=255)

                 ccAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGAc                                                                                                                                   <  1:1529812/150‑1 (MQ=255)

                   aGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCg                                                                                                                                 >  1:276028/1‑150 (MQ=255)

                     gTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGaa                                                                                                                               <  2:646708/150‑1 (MQ=255)

                       cTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATc                                                                                                                             <  2:143322/150‑1 (MQ=255)

                         gCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGTAGGCGAACCGTCCTGCCCTGTCGCCGTCAGGGCa                                                                                                                                       >  1:407221/1‑138 (MQ=255)

                         gCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGTAGGCGAACCGTCCTGCCCTGTCGCCGTCAGGGCa                                                                                                                                       <  2:407221/138‑1 (MQ=255)

                         gCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCa                                                                                                                           <  2:276028/150‑1 (MQ=255)

                          cACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGa                                                                                                                                    <  2:1456348/140‑1 (MQ=255)

                          cACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGa                                                                                                                                    >  1:1456348/1‑140 (MQ=255)

                           aCGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACt                                                                                                                         >  2:1305374/1‑150 (MQ=255)

                            cGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACtt                                                                                                                        <  1:315575/150‑1 (MQ=255)

                             gACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTTTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGa                                                                                                                                    >  2:1991149/1‑137 (MQ=255)

                             gACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTTTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGa                                                                                                                                    <  1:1991149/137‑1 (MQ=255)

                             gACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGCAg                                                                                                                                      >  1:2042244/1‑135 (MQ=255)

                             gACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGCAg                                                                                                                                      <  2:2042244/135‑1 (MQ=255)

                             gACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGAc                                                                                                                                   <  1:1443612/138‑1 (MQ=255)

                             gACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGAc                                                                                                                                   >  2:1443612/1‑138 (MQ=255)

                                      gTGCCAGCCACTGTCCTATTGGTCAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCccacca                                                                                                              >  2:1183581/1‑150 (MQ=255)

                                       tGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCcaccac                                                                                                             <  1:2066989/150‑1 (MQ=255)

                                          cagccaCTGTCCTATTGGTTAGGACCGCAGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTg                                                                                                          >  1:2175500/1‑150 (MQ=255)

                                          cagccaCTGTCCTATTGGTTAGGACCGCAGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTg                                                                                                          >  1:2175997/1‑150 (MQ=255)

                                            gccaCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCa                                                                                                        >  1:456779/1‑150 (MQ=255)

                                             ccaCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCaa                                                                                                       <  2:768119/150‑1 (MQ=255)

                                             ccaCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCaa                                                                                                       >  1:303696/1‑150 (MQ=255)

                                                      tATTGGTTAGGACCGCCGAGCGAACGCATTTCTACAGGAGACTTCTGGTCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAg                                                                                              >  1:461076/1‑150 (MQ=255)

                                                        ttGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGt                                                                                            >  2:1183595/1‑150 (MQ=255)

                                                        ttGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGt                                                                                            >  2:2411610/1‑150 (MQ=255)

                                                                  cggccgAGCTACCGCATGGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTa                                                                                  <  2:1958083/148‑1 (MQ=255)

                                                                   cgccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTTAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                                                 <  2:905444/150‑1 (MQ=255)

                                                                   cgccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                                                 >  1:2403662/1‑150 (MQ=255)

                                                                   cgccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                                                 <  1:1305374/150‑1 (MQ=255)

                                                                      cgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTATCCCGACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGa                                                                              >  1:623744/1‑150 (MQ=255)

                                                                          cGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCa                                                                          <  2:57777/150‑1 (MQ=255)

                                                                            aaCGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATGCGCCTTTGGACAGAAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGCATCCACTTCAGCCCACCACGTGCAACTAATGCAGGTGGCGTTCGTACTGATCCAgg                                                                        >  2:649786/1‑150 (MQ=255)

                                                                             aCGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGt                                                                                                           >  2:2264877/1‑114 (MQ=255)

                                                                             aCGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGt                                                                                                           <  1:2264877/114‑1 (MQ=255)

                                                                                 aTTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAATTTGACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTc                                                                   >  1:1879093/1‑150 (MQ=255)

                                                                                  ttGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCg                                                                  >  1:959754/1‑150 (MQ=255)

                                                                                     ctgcAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCaa                                                                                                       <  1:2163608/110‑1 (MQ=255)

                                                                                      tgcAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCaa                                                              >  1:1856294/1‑150 (MQ=255)

                                                                                           tAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATcaca                                                         >  1:2424409/1‑150 (MQ=255)

                                                                                           tAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATcaca                                                         >  1:1119359/1‑150 (MQ=255)

                                                                                                ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAggg                                                                                                                                         >  1:1725076/1‑65 (MQ=255)

                                                                                                ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAggg                                                                                                                                         <  2:1725076/65‑1 (MQ=255)

                                                                                                ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCGAATCACATCGTc                                                    >  2:1330565/1‑150 (MQ=255)

                                                                                                ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATcaca                                                         <  1:1417200/145‑1 (MQ=255)

                                                                                                ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATcaca                                                         >  2:1417200/1‑145 (MQ=255)

                                                                                                ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTc                                                    >  1:738697/1‑150 (MQ=255)

                                                                                                ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTc                                                    >  1:59029/1‑150 (MQ=255)

                                                                                                 tCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCa                                                   <  2:2386930/150‑1 (MQ=255)

                                                                                                     ggCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGc                                               >  1:1657869/1‑150 (MQ=255)

                                                                                                       ccTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGGTTATAAAGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCACTCACACCGTCACGGCTc                                             >  1:1594143/1‑150 (MQ=255)

                                                                                                        cTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCg                                                                  >  1:652488/1‑128 (MQ=255)

                                                                                                           tgAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                                                 <  1:2221877/110‑1 (MQ=255)

                                                                                                           tgAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                                                 >  2:2221877/1‑110 (MQ=255)

                                                                                                           tgAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCAcc                                         >  2:2407128/1‑150 (MQ=255)

                                                                                                           ggAAAATCCGCCTTTGGACAGCAGGCGTACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCAcc                                         <  2:1269861/149‑1 (MQ=255)

                                                                                                            gAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCAccc                                        >  1:1032265/1‑150 (MQ=255)

                                                                                                                 tCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGCTCCAGGATGTCGCCCATCACATCGTCACGGCCCCACCCCatgg                                   >  2:643003/1‑149 (MQ=255)

                                                                                                                 tCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCatga                                   <  1:1151707/150‑1 (MQ=255)

                                                                                                                 tCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGACCGATGATGTCGCCAATCACATCGTCACGGCTCCACCCCatga                                   >  1:1842179/1‑150 (MQ=255)

                                                                                                                   cGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCatgatg                                 >  1:867012/1‑150 (MQ=255)

                                                                                                                      cTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCg                              >  1:675599/1‑150 (MQ=255)

                                                                                                                          ggACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGt                          >  2:1372000/1‑150 (MQ=255)

                                                                                                                           gACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTc                         <  1:1183595/150‑1 (MQ=255)

                                                                                                                               gcagGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTgg                     <  1:510569/150‑1 (MQ=255)

                                                                                                                               gcagGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTgg                     <  2:867012/150‑1 (MQ=255)

                                                                                                                                cagGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGc                    >  2:42135/1‑150 (MQ=255)

                                                                                                                                 agGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGcc                   >  1:1621349/1‑150 (MQ=255)

                                                                                                                                   gCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGcccc                 >  2:729005/1‑150 (MQ=255)

                                                                                                                                      aaCCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCtt              >  2:226516/1‑150 (MQ=255)

                                                                                                                                       aCCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTtc             >  1:307252/1‑150 (MQ=255)

                                                                                                                                          gccctgccctGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCg          >  1:1768067/1‑150 (MQ=255)

                                                                                                                                           ccctgccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATcc                                                                           <  2:783043/84‑1 (MQ=255)

                                                                                                                                           ccctgccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATcc                                                                           >  1:783043/1‑84 (MQ=255)

                                                                                                                                           ccctgccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGc         >  1:1620784/1‑150 (MQ=255)

                                                                                                                                             ctgccctgTCGCCGTCAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAg                                                                         >  1:122946/1‑84 (MQ=38)

                                                                                                                                             ctgccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAg       >  1:2014908/1‑150 (MQ=255)

                                                                                                                                             ctgccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAg       <  2:635721/150‑1 (MQ=255)

                                                                                                                                              tgccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTtctc           <  2:1828471/145‑1 (MQ=255)

                                                                                                                                              tgccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGa      >  2:2032753/1‑150 (MQ=255)

                                                                                                                                              tgccccgAGGTCGACAGGGTAGGACCGAATTCACTTCAGCCCATCACGTGCTAGTAATGCAGGTGGCGTTCGTACTGACCCAGGATGTCGGTAATTACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGa      <  2:461076/150‑1 (MQ=255)

                                                                                                                                               ggcctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGGAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGAt     <  2:786443/148‑1 (MQ=255)

                                                                                                                                                ccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTc                                                                                     >  2:636372/1‑69 (MQ=255)

                                                                                                                                                ccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTc                                                                                     <  1:636372/69‑1 (MQ=255)

                                                                                                                                                ccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATg    >  1:1493732/1‑150 (MQ=255)

CCCGAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCA  >  NC_002947/4463848‑4464143

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