Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,511,734	C→T	6.1%	C134C (TGC→TGT)	PP_3108 →	rhs family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,511,734	0	C	T	6.1%	192.9 / 2.4	65	C134C (TGC→TGT)	PP_3108	rhs family protein
Reads supporting (aligned to +/- strand):  ref base C (26/35);  new base T (2/2);  total (28/37)
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.06e-01

GTTCATCAACAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGCACCACCGTGGATGTGCTGATGGCCATCGCCGGTGCCGCAGGCGGTA  >  NC_002947/3511587‑3511876                                                                                                                                                    |                                                                                                                                               gtTCATCAACAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTgcg                                                                                                                                               >  1:361021/1‑149 (MQ=255) gtTCATCAACAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTgcg                                                                                                                                               >  6:138734/1‑149 (MQ=255)  ttCATCAACAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTgcgg                                                                                                                                              <  2:72043/149‑1 (MQ=255)  ttCATCAACAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTgcgg                                                                                                                                              >  5:125373/1‑149 (MQ=255)     atcaACAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCATGGGCGACAACCTGCCCTGCGgcgc                                                                                                                                           >  7:161581/1‑149 (MQ=255)     atcaACAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGTGACAAGCTGACCTCCGgcgc                                                                                                                                           >  4:213051/1‑149 (MQ=255)        aaCAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGccac                                                                                                                                        >  4:90863/1‑149 (MQ=255)        aaCAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGcca                                                                                                                                         <  6:125373/148‑1 (MQ=255)              cGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTc                                                                                                                                  >  6:308977/1‑149 (MQ=255)              cGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTc                                                                                                                                  <  5:329052/149‑1 (MQ=255)              cGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCACGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTc                                                                                                                                  <  8:335433/149‑1 (MQ=255)               gCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTcc                                                                                                                                 >  4:130624/1‑149 (MQ=255)                 aaGGCCGCCCGCGCCGAAAAAAGCATCGGTGCCAGCGACAAGCAACCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCgg                                                                                                                               <  8:88513/149‑1 (MQ=255)                  aGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCTGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGTGGCGCCACCATTTCCggt                                                                                                                              >  8:37894/1‑149 (MQ=255)                  aGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCTGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGTGGCGCCACCATTTCCggt                                                                                                                              >  8:32160/1‑149 (MQ=255)                   ggCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCggtg                                                                                                                             >  5:318532/1‑149 (MQ=255)                      cccccgcGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACCTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGAGGtt                                                                                                                          >  5:364319/3‑149 (MQ=255)                       gcccgcGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTc                                                                                                                         <  4:30325/149‑1 (MQ=255)                        cccgcGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCg                                                                                                                        <  3:254589/149‑1 (MQ=255)                              tCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTTCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGTGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAAc                                                                                                                  <  3:213051/149‑1 (MQ=255)                               cGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACg                                                                                                                 >  8:165959/1‑149 (MQ=255)                                 aaaaaaGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGt                                                                                                                >  8:78777/1‑148 (MQ=255)                                    aaaGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGtcatc                                                                                                            <  5:308977/149‑1 (MQ=255)                                      aGCATCGGTGCCTTCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTcatcat                                                                                                          <  3:367905/149‑1 (MQ=255)                                              tGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATcggcggcg                                                                                                  <  1:214393/149‑1 (MQ=255)                                                 cTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCa                                                                                               >  7:20560/1‑149 (MQ=255)                                                  tGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCAc                                                                                              >  8:102142/1‑149 (MQ=255)                                                   gCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCAcc                                                                                             >  6:295282/1‑149 (MQ=255)                                                    cGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCAcc                                                                                             <  7:165959/148‑1 (MQ=255)                                                    cGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCAcc                                                                                             <  4:96261/148‑1 (MQ=255)                                                    cGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCAcc                                                                                             <  7:102142/148‑1 (MQ=255)                                                      aCAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTAc                                                                                          <  7:197112/149‑1 (MQ=255)                                                       cAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTAcc                                                                                         <  5:50736/149‑1 (MQ=255)                                                            aCCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTAc                                                                                    <  3:130624/149‑1 (MQ=255)                                                            aCCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTAc                                                                                    <  8:20560/149‑1 (MQ=255)                                                             cccAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGtcatc                                                                                                            <  5:91908/124‑1 (MQ=255)                                                             cccAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGtcatc                                                                                                            >  6:91908/1‑124 (MQ=255)                                                              ccAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGTTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCt                                                                                  >  2:90608/1‑149 (MQ=255)                                                                aGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTGCCGTTACCTGc                                                                                <  4:248032/149‑1 (MQ=255)                                                                     cGGTGCTGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGTGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTg                                                                           <  7:37894/149‑1 (MQ=255)                                                                     cGGTGCTGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGTGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTg                                                                           <  7:32160/149‑1 (MQ=255)                                                                          cAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGa                                                                      <  3:90863/149‑1 (MQ=255)                                                                          cAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCAACATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGa                                                                      <  5:362009/149‑1 (MQ=255)                                                                            gATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGAt                                                                     >  4:228943/1‑148 (MQ=255)                                                                            gATTGCCGAAGGCTCGACCAACGGCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATg                                                                    <  3:201275/149‑1 (MQ=255)                                                                              ttGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGaa                                                                  <  2:328820/149‑1 (MQ=255)                                                                              ttGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGaa                                                                  <  5:138734/149‑1 (MQ=255)                                                                              ttGCCGAAGGCTCGACCAACGTAGTCATCAACAGCGTAGCCGCCGCGCGCAAGGGCGACAAGCTGACATGCGGCGCCACCATTTACGGTGGTTCGGACACCGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGAtt                                                                    >  8:10921/1‑146 (MQ=255)                                                                                  cGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTcc                                                              >  1:11107/1‑149 (MQ=255)                                                                                   gAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTccc                                                             >  1:250331/1‑149 (MQ=255)                                                                                      ggCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGa                                                          <  7:78777/149‑1 (MQ=255)                                                                                      ggCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGa                                                          <  3:141565/149‑1 (MQ=255)                                                                                        cTCGACCTACGTCTTCATCAACAGCGTCGCCGGCGCGCGGAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTGATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGt                                                        <  6:364319/149‑1 (MQ=255)                                                                                         tCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGTAGTg                                                       >  8:50845/1‑149 (MQ=255)                                                                                           gACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGTTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGc                                                     <  1:90608/149‑1 (MQ=255)                                                                                              cAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGAACTGCGGCGCCACCATTTCCGGTCGTTCGGACAACGTCATCATCCGCGGCGGCACCTACCGTTACCTGCCCGTGGACCATGAAGTCCCGGAGTGGCTgc                                                  >  4:359931/1‑149 (MQ=255)                                                                                               aaCGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAAAGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTgcg                                                 >  8:266398/1‑149 (MQ=255)                                                                                                aCGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTgcgc                                                <  3:228943/149‑1 (MQ=255)                                                                                                  gTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGcac                                              <  6:318532/149‑1 (MQ=255)                                                                                                  gTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGCCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGcac                                              <  2:361021/149‑1 (MQ=255)                                                                                                    cTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGcacca                                            >  5:31799/1‑149 (MQ=255)                                                                                                     ttCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGcaccac                                           >  3:268865/1‑149 (MQ=255)                                                                                                       catcaACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGCACCACCg                                         >  8:220933/1‑149 (MQ=255)                                                                                                                cggcgccgccgcGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGGGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTTGACGATGAAGTCCCGGAGTGGCTGCGCACCACCGTGGATGTGc                                <  7:56931/146‑1 (MQ=255)                                                                                                                    gccgccgcGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGCACCACCGTGGATGTGCTGAt                            <  5:24203/149‑1 (MQ=255)                                                                                                                      cgccgcGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGCACCACCGTGGATGTGCTGATgg                          >  3:198428/1‑149 (MQ=255)                                                                                                                      cgccgcGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGCACCACCGTGGATGTGCTGATgg                          >  7:81651/1‑149 (MQ=255)                                                                                                                             cgcAAGGGCGAAAAGCTGACCTGCGGCGACACCATTTCCTGTGGTTCGTACAACGTCATCATCGGCGGCGGAACCTACCGTTACCTTCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGCACCACCGTGGATGTGCTAATGGCCATCGc                   <  7:10921/149‑1 (MQ=255)                                                                                                                                               aCCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGCACCACCGTGGATGTGCTGATGGCCATCGCCGGTGCCGCAGGCGGTa  <  2:350527/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               GTTCATCAACAGCCGCAAGGCCGCCCGCGTCGAAAAAAGCATCGGTGCCTGCGACAAGCACCCAGGGCCGGTGCAGATTGCCGAAGGCTCGACCAACGTCTTCATCAACAGCGTCGCCGCCGCGCGCAAGGGCGACAAGCTGACCTGCGGCGCCACCATTTCCGGTGGTTCGGACAACGTCATCATCGGCGGCGGCACCTACCGTTACCTGCCCGTGGACGATGAAGTCCCGGAGTGGCTGCGCACCACCGTGGATGTGCTGATGGCCATCGCCGGTGCCGCAGGCGGTA  >  NC_002947/3511587‑3511876

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