Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,843,873	C→A	6.1%	L241M (CTG→ATG) ‡	PP_4259 →	iron‑sulfur cluster‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,843,873	0	C	A	6.1%	209.3 / 2.0	65	L241M (CTG→ATG) ‡	PP_4259	iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base A (2/2);  total (34/31)
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.55e-01

CACGGGCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCC  >  NC_002947/4843725‑4844020                                                                                                                                                     |                                                                                                                                                       cACGGGCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGc                                                                                                                                                        >  6:260270/1‑149 (MQ=255)      gCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATc                                                                                                                                                   <  4:120798/149‑1 (MQ=255)         aGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCgtg                                                                                                                                                <  8:161957/149‑1 (MQ=255)           gCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTc                                                                                                                                              <  3:239192/149‑1 (MQ=255)            ccGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTcc                                                                                                                                             <  2:67941/149‑1 (MQ=255)             cGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGc                                                                                                                                                        >  8:298837/1‑137 (MQ=255)             cGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGc                                                                                                                                                        <  7:298837/137‑1 (MQ=255)             cGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCt                                                                                                                                            <  8:159509/149‑1 (MQ=255)              gATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                           <  7:246828/149‑2 (MQ=255)                 ggCTGGGCCTACTTCTGGGTCGGCTTGTTCAGCGTGGCCACTTACGGAAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCGCGCTTGCAGAGCGTGATGTTCGACAAGGAGACGATCAGCGTGTCCttgt                                                                                                                                         <  4:191606/148‑4 (MQ=25)                  gCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGa                                                                                                                                        <  1:299114/148‑1 (MQ=255)                   cTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGAcc                                                                                                                                      >  2:39195/1‑149 (MQ=255)                            cTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGt                                                                                                                              <  5:183432/148‑1 (MQ=255)                            cTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTg                                                                                                                             <  7:168078/149‑1 (MQ=255)                               cTGGGTCGGCTTCTTCACCCTCGCGACTTACGGCAACGCCTGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTgg                                                                                                                            <  5:322129/147‑1 (MQ=255)                               cTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGc                                                                                                                           <  6:187683/148‑1 (MQ=255)                                   gTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAAc                                                                                                                      <  1:139209/149‑1 (MQ=255)                                    tCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAAcc                                                                                                                     >  2:314534/1‑149 (MQ=255)                                       gCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCcgc                                                                                                                  >  5:67125/1‑149 (MQ=255)                                              cACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGccc                                                                                                             <  6:292368/147‑1 (MQ=255)                                              cACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCcgc                                                                                                           >  6:257387/1‑149 (MQ=255)                                                  cTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaa                                                                                                       <  4:295469/149‑1 (MQ=255)                                                   tCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCgtg                                                                                                                                                >  7:224286/1‑107 (MQ=255)                                                   tCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCgtg                                                                                                                                                <  8:224286/107‑1 (MQ=255)                                                      ccACTTACGGCAACTCCGGCTGGCTACGCTAACAGCTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaagaaa                                                                                                   <  4:167935/149‑1 (MQ=255)                                                      ccACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCCCGCCGTGGCGAAACCCGCGGCCCGCGCaaga                                                                                                     >  5:262948/1‑147 (MQ=255)                                                        aCTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTg                                                                                                                             <  8:41472/121‑1 (MQ=255)                                                        aCTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTg                                                                                                                             >  7:41472/1‑121 (MQ=255)                                                        aCTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaagaaag                                                                                                  >  8:94349/1‑148 (MQ=255)                                                         cTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaagaaaga                                                                                                 >  2:23721/1‑148 (MQ=255)                                                                       ggctggctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑Gc                                                                               >  5:64202/1‑149 (MQ=255)                                                                         ctggctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑Gcc                                                                              >  2:222365/1‑148 (MQ=255)                                                                            gctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCa                                                                            <  5:260270/147‑1 (MQ=255)                                                                            gctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAgg                                                                          <  6:265790/149‑1 (MQ=255)                                                                             ctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACGAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAggg                                                                         >  2:62830/1‑149 (MQ=255)                                                                                  cgAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTgg                                                                    <  7:296068/149‑1 (MQ=255)                                                                                     aCAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCTGGGCCttgtag                                                                 <  5:257387/149‑6 (MQ=255)                                                                                         gtgtGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGACTACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACt                                                              <  1:39195/149‑1 (MQ=255)                                                                                                tgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGAc                                                      >  4:88699/1‑149 (MQ=255)                                                                                                  tACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACt                                                     <  6:67125/148‑1 (MQ=255)                                                                                                       gtgCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCa                                               >  1:175355/1‑149 (MQ=255)                                                                                                        tgCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCAt                                              >  7:178964/1‑149 (MQ=255)                                                                                                        tgCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCAt                                              >  6:178804/1‑149 (MQ=255)                                                                                                         gCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGACGGACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCAt                                              >  1:108914/1‑149 (MQ=255)                                                                                                          ccccTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCAtgt                                            <  5:178804/149‑1 (MQ=255)                                                                                                           cccTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAGAGCCCAGGGCCTGGGTGACTGCATCGACTGCACCAtgt                                            <  1:222365/149‑1 (MQ=255)                                                                                                           cccTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCAtgtg                                           >  8:250540/1‑149 (MQ=255)                                                                                                            ccTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATgtgc                                          <  1:314534/149‑1 (MQ=255)                                                                                                                  cccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                    <  3:243499/149‑1 (MQ=255)                                                                                                                   ccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCCCCATGTGCGTGCAgg                                   >  6:218202/1‑149 (MQ=255)                                                                                                                   ccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAgg                                   >  6:11042/1‑149 (MQ=255)                                                                                                                   ccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAgg                                   >  5:74362/1‑149 (MQ=255)                                                                                                                           agagCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGccc                            >  1:68488/1‑148 (MQ=255)                                                                                                                                gTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCgg                      >  2:41806/1‑149 (MQ=255)                                                                                                                                      ttagACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATTTGCGTGCAGGTCTGCCCGACCGGCATCGa                >  2:295509/4‑149 (MQ=255)                                                                                                                                      ttCTACAAGGACACGCTATTCGTGTCCTACGACCCACGCCGTGTCGAAAGCCGCGGCCCGCGCAAGAATGA‑TAGCGAC‑TACTAG‑GCCCAGGGCCTGGGTGACTGCATCGACTTCACCATGTTCGTGCAGGTGTGCCCGACCGGCATCg                 <  2:317696/148‑1 (MQ=255)                                                                                                                                      ttCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGa                >  1:211724/1‑149 (MQ=255)                                                                                                                                         gACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGc                              <  8:119323/132‑1 (MQ=255)                                                                                                                                           cAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATcc           >  1:12706/1‑149 (MQ=255)                                                                                                                                            taGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCCGGTCTGCCCGACCGGCATCGACATCCg          >  6:176655/2‑149 (MQ=255)                                                                                                                                            aaGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                    <  7:126284/123‑1 (MQ=255)                                                                                                                                            aaGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                    >  8:126284/1‑123 (MQ=255)                                                                                                                                              ggACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGt         >  6:283763/1‑148 (MQ=255)                                                                                                                                               gACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGa       >  4:54451/1‑149 (MQ=255)                                                                                                                                                  acGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TCCCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACgg    >  4:243521/1‑149 (MQ=255)                                                                                                                                                    gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGcc  >  5:132361/1‑149 (MQ=255)                                                                                                                                                    gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGcc  >  4:256960/1‑149 (MQ=255)                                                                                                                                                    gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGcc  >  1:104456/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                       CACGGGCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGA‑TACCGAC‑TACAAG‑GCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCC  >  NC_002947/4843725‑4844020

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