Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,103,979 A→C 5.7% E271A (GAG→GCG)  PP_5353 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,103,9790AC 5.7% 276.9 / 6.8 88E271A (GAG→GCG) PP_5353hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base C (3/2);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

GATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGT  >  NC_002947/6103831‑6104120
                                                                                                                                                     |                                                                                                                                             
gATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCga                                                                                                                                               <  4:59204/149‑1 (MQ=255)
gATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGCTCGCGCCCAAGGCga                                                                                                                                               <  7:4900/149‑1 (MQ=255)
     ccTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAAc                                                                                                                                           >  4:59887/1‑148 (MQ=255)
      cTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAAc                                                                                                                                           >  5:243608/1‑147 (MQ=255)
             tGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCTCGCCTTGGGCGCGAACCatgtgcca                                                                                                                                  <  8:336343/149‑8 (MQ=255)
             tGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTgg                                                                                                                                  <  4:241805/149‑1 (MQ=255)
              gCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCTCGCCTTGGGCGCGAACCatgtgcca                                                                                                                                  <  4:307558/148‑8 (MQ=255)
               cAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGc                                                                                                                                 <  4:354260/148‑1 (MQ=255)
               cAAATGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTTTCGAAGATTTCGCGCTTCGACGAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATTGTTCGCGGCCAAGGCGAGAACCAGCCTTGGc                                                                                                                                 <  3:122498/148‑1 (MQ=255)
                aaCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTGCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTg                                                                                                                               <  3:34012/149‑1 (MQ=255)
                 aCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              <  6:249035/149‑1 (MQ=255)
                 aCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              >  5:247633/1‑149 (MQ=255)
                       cGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCTCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCGCAAGTCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        <  8:149030/149‑1 (MQ=255)
                       cGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        >  3:181409/1‑149 (MQ=255)
                        ggTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCgg                                                                                                                       <  5:251697/149‑1 (MQ=255)
                            tCGATGACCGCGCCGAGCGGACCGACGCTACCCTGACCTTCCGGGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCTAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCTCCCAAGGCGAGAACCAGCGTTGGCTGGTAGCCGGTAt                                                                                                                    <  4:324506/148‑1 (MQ=255)
                              gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              >  1:301422/1‑136 (MQ=255)
                              gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              <  2:301422/136‑1 (MQ=255)
                              gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCg                                                                                                                 >  4:185334/1‑149 (MQ=255)
                                   ccGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCACGGCGAA‑GTGTTCAGCGAAACCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCAGTCaa                                                                                                             >  6:174235/1‑148 (MQ=255)
                                       gcCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAg                                                                                                        <  3:132764/149‑1 (MQ=255)
                                       gcCCACCGCACCTAACGTACCCTGTCCTTACGTGGGGTTTCGAATAATTGGGGGTTTGCTCATGGGGAATGTCTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACa                                                                                                         <  1:117811/149‑1 (MQ=255)
                                               cACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGc                                                                                                >  3:296732/1‑149 (MQ=255)
                                               cACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGc                                                                                                >  7:127943/1‑149 (MQ=255)
                                                  cGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                             >  5:328896/1‑149 (MQ=255)
                                                  cGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                             >  5:133575/1‑149 (MQ=255)
                                                     cGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTg                                                                                          >  5:70200/1‑149 (MQ=255)
                                                          cccTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaa                                                                                      <  2:346478/148‑1 (MQ=255)
                                                          cccTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaa                                                                                      <  4:287399/148‑1 (MQ=255)
                                                           ccTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaaa                                                                                     >  7:345296/1‑148 (MQ=255)
                                                           ccTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaaa                                                                                     >  2:166267/1‑148 (MQ=255)
                                                             tGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCGTTGTCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                             <  3:351565/138‑1 (MQ=255)
                                                             tGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACAAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTCGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCgacaacgg                                                                                     >  4:351565/1‑139 (MQ=255)
                                                                ccTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAg                                                                                >  7:346079/1‑148 (MQ=255)
                                                                ccTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAg                                                                                >  5:342965/1‑148 (MQ=255)
                                                                  ttCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAc                                                                             >  4:276051/1‑149 (MQ=255)
                                                                  ttCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAc                                                                             >  3:63542/1‑149 (MQ=255)
                                                                     cGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAcccc                                                                          <  6:247633/149‑1 (MQ=255)
                                                                         gCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAATCAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                       <  6:243608/148‑1 (MQ=255)
                                                                         gCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                       <  2:254878/148‑1 (MQ=255)
                                                                         gCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                       <  8:132991/148‑1 (MQ=255)
                                                                         gCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTCGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCACGGGc                                                                      >  6:264426/1‑149 (MQ=255)
                                                                              tCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        <  4:95013/94‑1 (MQ=255)
                                                                              tCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        >  3:95013/1‑94 (MQ=255)
                                                                              tCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGcc                                                                 >  1:259659/1‑149 (MQ=255)
                                                                                gaagaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCg                                                               <  4:102833/149‑1 (MQ=255)
                                                                                 aagaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                              >  6:219508/1‑149 (MQ=255)
                                                                                 aagaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                              >  6:190046/1‑149 (MQ=255)
                                                                                 aagaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCAGGTATCCGTCAAAACAGCTAATCGCAGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                              >  6:346944/1‑149 (MQ=255)
                                                                                 aagaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCACCGAAAGCTGGCACATGGTAAGCGCCCAAGGGGAGAACCAGCCTTGGCTGGTAGCCGCTATCCGTCAAAACAGCTAATCGCTGCCTGCTAAACACAGACCCACGGGCTTGCCCgg                                                              >  4:31648/1‑149 (MQ=255)
                                                                                  agaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                             >  1:144031/1‑149 (MQ=255)
                                                                                  agaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                             <  8:346079/149‑1 (MQ=255)
                                                                                  agaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                             <  6:342965/149‑1 (MQ=255)
                                                                                  agaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCggg                                                             <  4:296732/149‑1 (MQ=255)
                                                                                   gaaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgggg                                                            >  2:95508/1‑149 (MQ=255)
                                                                                    aaTTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGt                                                           <  6:70200/149‑1 (MQ=255)
                                                                                         gcgcTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCGGGGttctgt                                                      <  1:145511/149‑5 (MQ=255)
                                                                                            cTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctt                                                   >  5:161083/1‑149 (MQ=255)
                                                                                               cGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttt                                                  <  5:190046/147‑1 (MQ=255)
                                                                                                gACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttttg                                                <  4:181409/148‑1 (MQ=255)
                                                                                                        cGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTATTCTTTTGCCAAAGTgg                                       >  7:7162/1‑149 (MQ=255)
                                                                                                          aa‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAac                                                                                    >  5:94509/1‑102 (MQ=255)
                                                                                                          aa‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAac                                                                                    <  6:94509/102‑1 (MQ=255)
                                                                                                               gtTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                  >  7:17414/1‑148 (MQ=255)
                                                                                                               gtTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                  >  5:213418/1‑148 (MQ=255)
                                                                                                                ttCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                  >  6:14685/1‑147 (MQ=255)
                                                                                                                   aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAggg                             <  3:276051/149‑1 (MQ=255)
                                                                                                                   aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  6:180329/148‑1 (MQ=255)
                                                                                                                   aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  8:345296/148‑1 (MQ=255)
                                                                                                                   aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  2:259659/148‑1 (MQ=255)
                                                                                                                   aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  5:219508/148‑1 (MQ=255)
                                                                                                                       aaaGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              >  2:34280/1‑144 (MQ=255)
                                                                                                                       aaaGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGCTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  1:34280/144‑1 (MQ=255)
                                                                                                                         aGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAg                       <  5:310025/149‑1 (MQ=255)
                                                                                                                          gCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAgg                      <  3:237052/149‑1 (MQ=255)
                                                                                                                           cTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgc                     >  8:160733/1‑149 (MQ=255)
                                                                                                                              aaggctGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgcg                    >  2:372894/6‑147 (MQ=255)
                                                                                                                               aggctGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACACCGGGCTTGGCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAg                 >  5:242158/5‑149 (MQ=255)
                                                                                                                                acaTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGc                <  8:17414/149‑1 (MQ=255)
                                                                                                                                acaTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGc                >  2:127263/1‑149 (MQ=255)
                                                                                                                                   tGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt             >  8:99709/1‑149 (MQ=255)
                                                                                                                                     gTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt           <  5:224881/149‑1 (MQ=255)
                                                                                                                                     gTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt           >  6:71724/1‑149 (MQ=255)
                                                                                                                                      ttCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt           >  7:128062/1‑148 (MQ=255)
                                                                                                                                       tCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCTGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTc          <  6:328896/148‑1 (MQ=255)
                                                                                                                                         gcgcCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATc       <  6:213418/149‑1 (MQ=255)
                                                                                                                                            cccAATGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAg    <  3:59887/149‑1 (MQ=255)
                                                                                                                                             ccAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATACGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCAGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAgg   >  5:298397/1‑149 (MQ=255)
                                                                                                                                              cAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGt  <  8:128062/149‑1 (MQ=255)
                                                                                                                                              cAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGt  <  1:166267/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGT  >  NC_002947/6103831‑6104120

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