Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,197,917 T→C 6.4% N223D (AAC→GAC)  PP_4578 ← MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,197,9170TC 6.4% 185.5 / 2.6 62N223D (AAC→GAC) PP_4578MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (25/33);  new base C (2/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GTTCGTTCTTCAGGAAGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAGGAGAACAGCAGCGCATAAAG  >  NC_002947/5197773‑5198064
                                                                                                                                                 |                                                                                                                                                   
gttcgttcTTCAGGAAGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTcc                                                                                                                                                 <  6:129434/149‑1 (MQ=255)
    gttcTTCAGGAAGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACt                                                                                                                                              <  3:195251/148‑1 (MQ=255)
     ttcttcAGGAAGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTgc                                                                                                                                            <  5:75194/149‑1 (MQ=255)
      tcttcAGGAAGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTgcc                                                                                                                                           >  7:209641/1‑149 (MQ=255)
          cAGGAAGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAGCCGCGAAGGTGCCTCGGCGTTTTCCACTgcct                                                                                                                                          <  1:34109/147‑1 (MQ=255)
          cAGGAAGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTgcctt                                                                                                                                         <  1:269281/147‑1 (MQ=255)
            ggAAGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCCAGCGCCCCGGTCGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTgccttggcc                                                                                                                                     >  1:217373/1‑149 (MQ=255)
               aGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTgccttggccttg                                                                                                                                  <  6:290039/149‑1 (MQ=255)
                 gTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGc                                                                                                                                 >  8:14165/1‑148 (MQ=255)
                     gCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTgcc                                                                                                                                           >  4:99999/1‑134 (MQ=255)
                     gCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTgcc                                                                                                                                           <  3:99999/134‑1 (MQ=255)
                     gCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGa                                                                                                                               >  6:97125/1‑146 (MQ=255)
                     gCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATa                                                                                                                             <  6:33944/148‑1 (MQ=255)
                       aGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGAc                                                                                                                          <  5:242385/149‑1 (MQ=255)
                           aCGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCa                                                                                                                      <  8:248974/149‑1 (MQ=255)
                             gAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAgg                                                                                                                    <  8:209641/149‑1 (MQ=255)
                             gAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAgg                                                                                                                    <  3:85653/149‑1 (MQ=255)
                             gAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAgg                                                                                                                    <  5:105681/149‑1 (MQ=255)
                                gTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTc                                                                                                                 >  6:794/1‑149 (MQ=255)
                                 tGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTcc                                                                                                                <  8:112374/149‑1 (MQ=255)
                                 tGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTcc                                                                                                                <  1:210877/149‑1 (MQ=255)
                                   aTGGCGTAGTAGCCGCGGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCGGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCtt                                                                                                              <  1:361669/149‑1 (MQ=255)
                                      gCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGAc                                                                                                           <  5:250805/149‑1 (MQ=255)
                                      gCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCAAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGa                                                                                                            <  4:308360/148‑1 (MQ=255)
                                             agCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAggcgg                                                                                                    >  6:147050/1‑149 (MQ=255)
                                              gccgccGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAggcggc                                                                                                   >  1:53074/1‑149 (MQ=255)
                                              gccgccGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAggcggc                                                                                                   <  4:19863/149‑1 (MQ=255)
                                                cgccgAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAgg                                                                                                                    <  8:63686/130‑1 (MQ=255)
                                                cgccgAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAgg                                                                                                                    >  7:63686/1‑130 (MQ=255)
                                                 gccgAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGcgaa                                                                                                >  4:309334/1‑149 (MQ=255)
                                                 gccgAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGcgaa                                                                                                >  5:245495/1‑149 (MQ=255)
                                                         ccccGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGAt                                                                                         >  1:127451/1‑148 (MQ=255)
                                                          cccGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCa                                                                                       >  8:159520/1‑149 (MQ=255)
                                                                 gtCAGCATGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGAt                                                                                <  5:185376/149‑1 (MQ=255)
                                                                    agcagGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTCTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGc                                                                             >  5:369345/1‑149 (MQ=255)
                                                                     gcagGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGcc                                                                            <  5:97125/149‑1 (MQ=255)
                                                                     gcagGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGcc                                                                            <  7:159520/149‑1 (MQ=255)
                                                                      cagGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCg                                                                           <  5:359380/149‑1 (MQ=255)
                                                                       agGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCgg                                                                          >  1:228950/1‑149 (MQ=255)
                                                                        gggATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTCTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCgg                                                                          >  8:225801/1‑148 (MQ=255)
                                                                            tGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCa                                                                     <  7:14165/149‑1 (MQ=255)
                                                                               ccGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGcc                                                                  >  2:91567/1‑149 (MQ=255)
                                                                                     aGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCcagca                                                            <  8:314678/149‑1 (MQ=255)
                                                                                          tCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGaa                                                        >  5:161664/1‑148 (MQ=255)
                                                                                          tCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGaa                                                        >  7:194152/1‑148 (MQ=255)
                                                                                            cAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTCTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGaacaa                                                     <  6:369345/149‑1 (MQ=255)
                                                                                                gCATGCCCGGGGCGAAGATTTAGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAactca                                                 <  6:312294/149‑1 (MQ=255)
                                                                                                    gCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAacgcacgc                                              >  5:216868/1‑148 (MQ=255)
                                                                                                            cGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGCCCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGc                                     >  8:68217/1‑149 (MQ=255)
                                                                                                                    tCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCgg                             >  2:338478/1‑149 (MQ=255)
                                                                                                                      gTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCgg                             >  3:241046/1‑147 (MQ=255)
                                                                                                                            g‑cgcGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAgg                     <  2:127451/149‑1 (MQ=255)
                                                                                                                               gcgAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAg                      >  4:131284/1‑146 (MQ=255)
                                                                                                                               gcgAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAGgag                   <  3:309334/149‑1 (MQ=255)
                                                                                                                                  aaGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCTATAGACCTCAGGGTCCTTGACCAGGCGCCGAACGAAGATCACGTAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAGgagatc                <  2:217373/149‑3 (MQ=255)
                                                                                                                                      tGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAGGAGAAcag              >  5:197442/1‑147 (MQ=255)
                                                                                                                                        ccTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAGGAGAACagcagc          <  1:91567/149‑1 (MQ=255)
                                                                                                                                        ccTCGGCGTCTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAGGAGAAcagcag           <  7:225801/148‑1 (MQ=255)
                                                                                                                                         cTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCa                       <  7:66309/135‑1 (MQ=255)
                                                                                                                                         cTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCa                       >  8:66309/1‑135 (MQ=255)
                                                                                                                                          tCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCCCGCCAGGCCTGCTCTGCCGGCAGCCAGGAGAACAGCAgcg         >  6:378879/1‑148 (MQ=255)
                                                                                                                                           cGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAGGAGAACAGCAGCGCa       <  6:245495/149‑1 (MQ=255)
                                                                                                                                                ttttCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAGGAGAACAGCAGCGCATAAAg  <  2:53074/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTTCGTTCTTCAGGAAGGTTGGCAGCCACGAGGTGATGGCGTAGTAGCCGCCGAGCGCCCCGGTGGTCAGCAGGGATGCCCGCACAGTGGTCCAGAGCATGCCCGGGGCGAAGATTTCGTAGAAG‑CGCGAAGGTGCCTCGGCGTTTTCCACTGCCTTGGCCTTGCGATAGACCTCAGGGTCCTTGACCAGGCGGCGAACGAAGATCACGAAGATTGCCGGTACCAGCCCCAGCAGGAACAACGCACGCCAGGCCTGCTCTGCCGGCAGCCAGGAGAACAGCAGCGCATAAAG  >  NC_002947/5197773‑5198064

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