Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,767,350 C→A 8.5% M1174I (ATG→ATT) ‡ PP_5060 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,767,3500CA 8.5% 206.6 / 9.7 71M1174I (ATG→ATT) ‡PP_5060hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base A (2/4);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

TTCGGCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATTGATG  >  NC_002947/5767203‑5767497
                                                                                                                                                    |                                                                                                                                                   
ttcgGCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGccc                                                                                                                                                     <  4:214326/148‑1 (MQ=255)
ttcgGCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGccc                                                                                                                                                     >  3:214326/1‑148 (MQ=255)
    gCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGc                                                                                                                                                >  1:121547/1‑149 (MQ=255)
         aGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGa                                                                                                                                           <  8:200202/149‑1 (MQ=255)
         aGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggt                                                                                                                                           <  5:73457/149‑2 (MQ=255)
          ggCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggtc                                                                                                                                          <  5:274981/149‑3 (MQ=255)
          ggCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTgg                                                                                                                                            <  2:86003/147‑1 (MQ=255)
          ggCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACACTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTgg                                                                                                                                            >  1:86003/1‑147 (MQ=255)
          ggCAACTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggtc                                                                                                                                          <  4:245504/149‑3 (MQ=255)
           gCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGaa                                                                                                                                          >  2:250602/1‑148 (MQ=255)
                tCACAGAGCATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCAAGGCCTGGAAAAACCg                                                                                                                                    <  5:19773/149‑1 (MQ=255)
                 cacaGAACATCAGCAGTTGATGCGTTTGGGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTACGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAAcc                                                                                                                                     <  1:95190/148‑1 (MQ=255)
                  acaGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGa                                                                                                                                  >  8:245444/1‑149 (MQ=255)
                   caGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGAt                                                                                                                                 <  2:121547/149‑1 (MQ=255)
                    aGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGCCGCTGACAGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGCGTA‑TGCGTGCGTTTGTGCCACTGCAGGGCCAGCTCCAACACGCACAATGCCTGCAAAAAGCGCAtc                                                                                                                                >  1:304265/1‑148 (MQ=255)
                    aGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATg                                                                                                                                >  3:212411/1‑149 (MQ=255)
                    aGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATg                                                                                                                                >  6:1691/1‑149 (MQ=255)
                       aCATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGtt                                                                                                                              <  4:91920/148‑1 (MQ=255)
                        cATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGa                                                                                                                            <  7:245444/149‑1 (MQ=255)
                                 atgatgCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGgtcg                                                                                                                   >  7:86714/1‑149 (MQ=255)
                                        gTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACg                                                                                                            >  7:257441/1‑149 (MQ=255)
                                        gTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGACGTCTTACg                                                                                                            >  3:137192/1‑149 (MQ=255)
                                             tgCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTAc                                                                                                             >  2:209451/1‑143 (MQ=255)
                                             tgCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTAc                                                                                                             <  1:209451/143‑1 (MQ=255)
                                              gCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTcgc                                                                                                      >  5:268328/1‑149 (MQ=255)
                                               cGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTcgc                                                                                                      <  6:268328/148‑1 (MQ=255)
                                                      ccGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGc                                                                                              >  7:44998/1‑149 (MQ=255)
                                                        gCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTc                                                                                            <  8:257441/149‑1 (MQ=255)
                                                        gCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTc                                                                                            <  7:65532/149‑1 (MQ=255)
                                                          aGACGCTGACTGGCGAGCGCATGGCGGAAGATCGTGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGTAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCga                                                                                          <  1:243576/149‑1 (MQ=255)
                                                          aGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCga                                                                                          <  1:250602/149‑1 (MQ=255)
                                                          aGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCga                                                                                          <  6:188002/149‑1 (MQ=255)
                                                              gCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATa                                                                                                  >  6:292506/1‑137 (MQ=255)
                                                              gCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATa                                                                                                  <  5:292506/137‑1 (MQ=255)
                                                              gCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACg                                                                                      >  6:152388/1‑149 (MQ=255)
                                                               cTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGACGTCTTACGGCTCGCGATAATGCTCGAGACgg                                                                                     <  4:137192/149‑1 (MQ=255)
                                                                    tGGCTAGCGCATTGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCaga                                                                                <  7:292131/149‑1 (MQ=255)
                                                                            gcATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTAtt                                                                        <  4:320578/149‑1 (MQ=255)
                                                                            gcATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTAtt                                                                        >  3:196554/1‑149 (MQ=255)
                                                                             cATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTg                                                                       <  4:212411/149‑1 (MQ=255)
                                                                               tGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGAc                                                                     <  5:152388/149‑1 (MQ=255)
                                                                                   ggAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTAAGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGt                                                                   >  4:74/1‑147 (MQ=255)
                                                                                    gAAGATCATGTAGTGCTGGGTG‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAg                                                                >  5:335345/1‑149 (MQ=255)
                                                                                    gAAGATCATGTAGTGCTGGGTG‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAg                                                                >  5:330893/1‑149 (MQ=255)
                                                                                           aTGTTGTGGTGGGTA‑GTGGTGCGTTTGTGCCAGTGCAGGGCCAGCTCTAACAGGCCCATGGACTTGAAAAACCGGATGTTGATCGTGGTCGTGTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCAGtt                                                          <  2:304265/148‑1 (MQ=255)
                                                                                            tGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTg                                                         >  1:204323/1‑148 (MQ=255)
                                                                                                    tGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACtt                                                          <  3:195296/139‑1 (MQ=255)
                                                                                                    tGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACtt                                                          >  4:195296/1‑139 (MQ=255)
                                                                                                     gggTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgtcgagg                                               >  4:38146/1‑147 (MQ=255)
                                                                                                      ggTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgtcgatgt                                              >  5:244663/1‑149 (MQ=255)
                                                                                                      ggTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgtcgat                                                <  2:204323/147‑1 (MQ=255)
                                                                                                              gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACTGTCAGATTTGTATTGACGTAAg                                                                <  2:317446/124‑1 (MQ=255)
                                                                                                              gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTTGTATTGACGTAAg                                                                >  1:317446/1‑124 (MQ=255)
                                                                                                              gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGAt                                       >  2:33304/1‑149 (MQ=255)
                                                                                                               tGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCaga                                                                                <  8:159882/107‑1 (MQ=255)
                                                                                                               tGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCaga                                                                                >  7:159882/1‑107 (MQ=255)
                                                                                                               tGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATg                                      <  4:196554/149‑1 (MQ=255)
                                                                                                                 ccTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGc                                    >  5:173575/1‑149 (MQ=255)
                                                                                                                       tgCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACa                               <  3:74/148‑1 (MQ=255)
                                                                                                                       tgCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCACGATAA‑GCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAgg                             <  8:86714/149‑1 (MQ=255)
                                                                                                                           aCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGa                          >  7:226727/1‑149 (MQ=255)
                                                                                                                           aCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCAAACAGGCGa                          >  5:137567/1‑149 (MQ=255)
                                                                                                                           aCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGACCAACAGGCGa                          <  5:1691/149‑1 (MQ=255)
                                                                                                                           aCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGAGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGa                          >  1:87321/1‑149 (MQ=255)
                                                                                                                                  ggCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAgct                   <  2:87321/149‑1 (MQ=255)
                                                                                                                                     cAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCat                                                      <  3:56430/111‑1 (MQ=255)
                                                                                                                                     cAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCat                                                      >  4:56430/1‑111 (MQ=255)
                                                                                                                                     cAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAgctggc                >  5:200902/1‑149 (MQ=255)
                                                                                                                                       gtttttccAGGCCATGGGCCTGGAAAAATCGGATGTTGATCGTGGTCGTCTTACGGCTCGCAATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAgctggctg              >  3:287415/7‑149 (MQ=255)
                                                                                                                                         tCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTTGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGga            <  6:173575/149‑1 (MQ=255)
                                                                                                                                           gACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGga            <  8:226727/147‑1 (MQ=255)
                                                                                                                                            tccAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGgaga          >  4:53114/2‑148 (MQ=255)
                                                                                                                                              cAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACAt       >  5:226662/1‑149 (MQ=255)
                                                                                                                                                  cccATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtgat   <  6:335345/149‑1 (MQ=255)
                                                                                                                                                  cccATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtgat   <  6:330893/149‑1 (MQ=255)
                                                                                                                                                   ccATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCTATAATGCTCGAGACGGTCAGAGTGGTATTGAGTTTAGGCACTTGGATGTCGATGTCGGTGATGGCCTTCAGGCGATCGCGGTGGCTGGAGACATtgatg  >  6:86865/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCGGCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTA‑GGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATTGATG  >  NC_002947/5767203‑5767497

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