Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,478,830 G→C 7.5% I33M (ATC→ATG)  PP_4816 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,478,8300GC 7.5% 288.4 / 10.6 93I33M (ATC→ATG) PP_4816hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/45);  new base C (2/5);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CAATGTCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTAGACCGAT  >  NC_002947/5478682‑5478977
                                                                                                                                                      |                                                                                                                                                   
cAATGTCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCgatga                                                                                                                                                    <  2:2323786/150‑1 (MQ=255)
  aTGTCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATc                                                                                                                                                  <  1:1637460/150‑1 (MQ=255)
   tGTCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCg                                                                                                                                                 <  2:1384822/150‑1 (MQ=255)
    gTCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATCATCgg                                                                                                                                                >  2:1096846/1‑150 (MQ=255)
     tCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGc                                                                                                                                               >  1:716624/1‑150 (MQ=255)
        ttGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATCATCGGcctg                                                                                                                                            <  1:1096846/150‑3 (MQ=255)
         tGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCgg                                                                                                                                           <  2:716624/150‑1 (MQ=255)
         tGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCgg                                                                                                                                           <  2:1830645/150‑1 (MQ=255)
          ggCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCggg                                                                                                                                          <  2:1554355/150‑1 (MQ=255)
          ggCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCggg                                                                                                                                          <  2:1046621/150‑1 (MQ=255)
          ggCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATCATCGGcctgtt                                                                                                                                          <  1:1687756/150‑5 (MQ=255)
          ggCATATACCAGCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCggg                                                                                                                                          <  1:1846632/150‑1 (MQ=255)
            cATATACCAGCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATCATCGGcctgttgt                                                                                                                                        <  2:1151816/150‑7 (MQ=255)
             atatACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACGGGCCGATCATCGGcctgttgtt                                                                                                                                       <  1:2047581/150‑8 (MQ=255)
                  ccATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACg                                                                                                                                  >  1:1607182/1‑150 (MQ=255)
                   cATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACgg                                                                                                                                 <  2:1458534/150‑1 (MQ=255)
                    aTCACCCCGACTGCCGGCGGCACCAGCAGCCGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGc                                                                                                                                >  2:2376518/1‑150 (MQ=255)
                    aTCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATc                                                                                                                                                  <  1:1067385/132‑1 (MQ=255)
                    aTCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATc                                                                                                                                                  >  2:1067385/1‑132 (MQ=255)
                    aTCACCCCGACTGCCGGCGGCACCAGCAGAAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCGCGGGGAAAAACGGc                                                                                                                                >  2:499846/1‑150 (MQ=255)
                      cACCCCGACGGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAc                                                                                                                              <  2:1128555/150‑1 (MQ=255)
                         cccGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGaa                                                                                                                           >  1:2076909/1‑150 (MQ=255)
                            gACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAg                                                                                                                        <  1:1325008/150‑1 (MQ=255)
                               tGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGtg                                                                                                                     >  1:1534009/1‑150 (MQ=255)
                                      ggcACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATc                                                                                                              <  2:2076909/150‑1 (MQ=255)
                                         aCCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCgg                                                                                                                                           >  1:793568/1‑118 (MQ=255)
                                         aCCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCgg                                                                                                                                           <  2:793568/118‑1 (MQ=255)
                                         aCCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGcc                                                                                                           <  1:1041119/150‑1 (MQ=255)
                                         aCCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTAATTGATTGct                                                                                                           <  2:645420/150‑2 (MQ=255)
                                           cagcagcagCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGGTCAGCGGCTTGCCGAACTAACCGCGGGACCAGCGCTGCACC‑‑CAGCAATGCAAGAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGGGGGGTCGCCGt                                                                                                         >  1:2350520/1‑150 (MQ=255)
                                                gcagcGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCca                                                                                                    >  2:1766218/1‑150 (MQ=255)
                                                   gcGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCcagcg                                                                                                 >  1:2271245/1‑149 (MQ=255)
                                                      ggAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCg                                                                                                                                            <  2:2243811/104‑1 (MQ=255)
                                                      ggAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCg                                                                                                                                            >  1:2243811/1‑104 (MQ=255)
                                                       gagaTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCa                                                                                                                               >  1:5416/1‑116 (MQ=255)
                                                       gagaTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCa                                                                                                                               <  2:5416/116‑1 (MQ=255)
                                                       gagaTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTg                                                                                             <  2:487877/150‑1 (MQ=255)
                                                          aTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGaa                                                                                          <  2:1534009/150‑1 (MQ=255)
                                                          aTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGaa                                                                                          <  2:1046282/150‑1 (MQ=255)
                                                          aTAGCGCAGCAAGCTGCTTGCGGCCTCGTCCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCCGCAA‑‑CAACCGGCCGATGATCGGCCCGGGTAAAACCGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGaa                                                                                          <  2:1897154/150‑1 (MQ=255)
                                                               gcagcaAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTCAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTa                                                                                     <  2:341958/150‑1 (MQ=255)
                                                               gcagcaAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATCATCGGcctgtt                                                                                                                                          >  1:432987/1‑93 (MQ=255)
                                                               gcagcaAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATCATCGGcctgtt                                                                                                                                          <  2:432987/97‑5 (MQ=255)
                                                                 agcaAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCcagca                                                                                                 <  2:1685734/136‑1 (MQ=255)
                                                                 agcaAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCcagca                                                                                                 >  1:1685734/1‑136 (MQ=255)
                                                                   caAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATc                                                                                                                                                  >  1:89488/1‑85 (MQ=255)
                                                                   caAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATc                                                                                                                                                  <  2:89488/85‑1 (MQ=255)
                                                                    aaGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCcagca                                                                                                 <  2:642682/133‑1 (MQ=255)
                                                                    aaGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCcagca                                                                                                 >  1:642682/1‑133 (MQ=255)
                                                                       ctgctTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTAccagccag                                                                             <  1:1566908/150‑1 (MQ=255)
                                                                        tgctTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGACGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATTGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAgg                                                                            <  1:100293/150‑1 (MQ=255)
                                                                           ttGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCa                                                                         >  1:196704/1‑150 (MQ=255)
                                                                             gCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCaaa                                                                       >  2:1764814/1‑150 (MQ=255)
                                                                               ggCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAAcc                                                                     >  2:1072169/1‑150 (MQ=255)
                                                                                 ccTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATc                                                                                                              <  1:315888/107‑1 (MQ=255)
                                                                                 ccTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATc                                                                                                              >  2:315888/1‑107 (MQ=255)
                                                                                     gTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCcagca                                                                                                 <  2:1994016/116‑1 (MQ=255)
                                                                                     gTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCcagca                                                                                                 >  1:1994016/1‑116 (MQ=255)
                                                                                        cAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAg                                                            <  2:2033616/150‑1 (MQ=255)
                                                                                         aGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTg                                                                                             >  2:432377/1‑116 (MQ=255)
                                                                                         aGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTg                                                                                             <  1:432377/116‑1 (MQ=255)
                                                                                         aGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGc                                                           <  1:1766218/150‑1 (MQ=255)
                                                                                         aGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGc                                                           <  2:2350520/150‑1 (MQ=255)
                                                                                             gCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATaa                                                       >  2:1282155/1‑150 (MQ=255)
                                                                                                 gCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCAt                                                   >  1:1275377/1‑150 (MQ=255)
                                                                                                  ccGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCAtc                                                  <  2:59207/150‑1 (MQ=255)
                                                                                                    gACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTAccagccag                                                                             >  1:1663223/1‑121 (MQ=255)
                                                                                                    gACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTAccagccag                                                                             <  2:1663223/121‑1 (MQ=255)
                                                                                                     aCCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGaa                                                                                          >  1:2022308/1‑107 (MQ=255)
                                                                                                     aCCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGaa                                                                                          <  2:2022308/107‑1 (MQ=255)
                                                                                                       cTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCt                                             >  1:2380990/1‑150 (MQ=255)
                                                                                                        tCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTc                                            >  2:1336961/1‑150 (MQ=255)
                                                                                                               gggCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGTACGAACAGGtgatt                                                                                                                  <  1:2284743/73‑1 (MQ=255)
                                                                                                               gggCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGTACGAACAGGtgatt                                                                                                                  >  2:2284743/1‑73 (MQ=255)
                                                                                                                 gCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGgc                                   >  2:2366561/1‑150 (MQ=255)
                                                                                                                    aTGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGGTCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGgccgc                                >  1:780687/1‑150 (MQ=255)
                                                                                                                    aTGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGgccgc                                <  2:2380990/150‑1 (MQ=255)
                                                                                                                     tGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTgg                                                                                            <  1:1284965/89‑1 (MQ=255)
                                                                                                                     tGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTgg                                                                                            >  2:1284965/1‑89 (MQ=255)
                                                                                                                        ggAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAAc                                                              >  1:1995198/1‑116 (MQ=255)
                                                                                                                        ggAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAAc                                                              <  2:1995198/116‑1 (MQ=255)
                                                                                                                         gAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCcagca                                                                                                 >  1:895063/1‑80 (MQ=255)
                                                                                                                         gAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCcagca                                                                                                 <  2:895063/80‑1 (MQ=255)
                                                                                                                         gAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATta                           <  1:554637/150‑1 (MQ=255)
                                                                                                                         gAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATta                           <  1:1336961/150‑1 (MQ=255)
                                                                                                                           agagCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTAccagccag                                                                             <  1:2256803/98‑1 (MQ=255)
                                                                                                                           agagCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTAccagccag                                                                             >  2:2256803/1‑98 (MQ=255)
                                                                                                                               cagcagcaA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCa                     >  2:1061863/1‑150 (MQ=255)
                                                                                                                                  cagcaA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCAGCCATTATAAGCACGc                  >  2:927744/1‑150 (MQ=255)
                                                                                                                                    gcaA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTc                >  1:692279/1‑150 (MQ=255)
                                                                                                                                    gcaA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTc                >  2:1937841/1‑150 (MQ=255)
                                                                                                                                     caa‑‑caaCAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCg               >  1:1512738/1‑150 (MQ=255)
                                                                                                                                             caGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTg                                                                                             <  1:1860943/66‑1 (MQ=255)
                                                                                                                                             caGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTg                                                                                             >  2:1860943/1‑66 (MQ=255)
                                                                                                                                                   gatgatCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTAGACCGa   >  2:1530384/1‑150 (MQ=255)
                                                                                                                                                    atgatCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCAt                                                   <  1:963875/101‑1 (MQ=255)
                                                                                                                                                    atgatCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCAt                                                   >  2:963875/1‑101 (MQ=255)
                                                                                                                                                    atgatCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTAGACCGat  <  1:1061863/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CAATGTCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAA‑‑CAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTAGACCGAT  >  NC_002947/5478682‑5478977

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