Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 675,592 C→G 9.1% G37A (GGC→GCC)  PP_0577 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947675,5920CG 9.1% 257.3 / 12.1 85G37A (GGC→GCC) PP_0577hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/41);  new base G (3/5);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTCCAT  >  NC_002947/675443‑675741
                                                                                                                                                     |                                                                                                                                                     
aCTGGCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGcagc                                                                                                                                                       <  2:1465416/150‑1 (MQ=255)
  tGGCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCAGGGcagcaa                                                                                                                                                     <  1:179988/150‑3 (MQ=255)
  tGGCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCAGGGcagcaa                                                                                                                                                     <  1:857084/150‑3 (MQ=255)
    gCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGGCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGcagccagc                                                                                                                                                   >  1:2154728/1‑150 (MQ=255)
    gCCGTTCCGCCGGCTCCAGCTGCGGGAGTGCCCATAATCACGAGCGAGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGGCTGGGCGGGGGCCGGTGGTAACGCGCGGCTGCTGGCCACATACGGGTTGCTGCCCTTGcatccagc                                                                                                                                                   >  2:1955158/1‑150 (MQ=255)
       gTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGccagccag                                                                                                                                                <  2:2362659/150‑1 (MQ=255)
        ttGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCTCTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAgg                                                                                                                                               <  2:327026/150‑1 (MQ=255)
        ttGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAgg                                                                                                                                               >  1:329267/1‑150 (MQ=255)
          gcgcCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGc                                                                                                                                             <  2:1937917/150‑1 (MQ=255)
           cgcCAGCTCCAGCTGCGGTAGTGGCCGCGACCACGAGCGACCGCAAGGTAAGCACTGGCGTTAAAGGTCTGCTCCTCCTGGGC‑GGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCTGGCTGCCAGGGCag                                                                                                                                           <  1:44411/150‑2 (MQ=255)
           cgcCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCTGGCTGCCAGGGCa                                                                                                                                            <  2:238479/150‑1 (MQ=255)
            gcCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCGGGGTAAGCACTGGCGTAGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGGCGCTGCCCTGGCGGCCGGCCAGGGCaa                                                                                                                                           >  1:859352/1‑150 (MQ=255)
            gcCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCaa                                                                                                                                           >  2:1025096/1‑150 (MQ=255)
             ccagctccagctGCGGTAGTGCCCATAATCGCGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAAt                                                                                                                                          <  2:202221/150‑1 (MQ=255)
             ccagctccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAAt                                                                                                                                          >  2:525664/1‑150 (MQ=255)
             ccagctccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGCTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAAt                                                                                                                                          >  2:2349838/1‑150 (MQ=255)
             ccagctccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGACCGGTGGTAACGGGCGGCTGCTGGCCCCATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAAt                                                                                                                                          >  2:317152/1‑150 (MQ=255)
             ccagctccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGCGGTGCCAGCC‑AGGAGGCCGAAGCTCGTACCGGTCTGCTGCTGGCCATTGACGTGGTTCTGCTCTGGCAGCCAGCCAGGGCAATg                                                                                                                                         <  2:2080845/150‑1 (MQ=255)
              cagctccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCGCTGTTGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATg                                                                                                                                         <  1:1955158/150‑1 (MQ=255)
              cagctccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATg                                                                                                                                         >  2:2407932/1‑150 (MQ=255)
                  tccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACCCGAGTGGTGCTGCCCTGGCAGCCAGCCAGGGCAATGGgca                                                                                                                                     >  1:1330995/1‑150 (MQ=255)
                     agctGCGGTAGTGCCCATGATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGgcagca                                                                                                                                  <  2:1211436/150‑1 (MQ=255)
                      gctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTTGGGGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGcagcag                                                                                                                                 >  1:2352268/1‑150 (MQ=255)
                      gctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGcagcag                                                                                                                                 >  2:1105628/1‑150 (MQ=255)
                       ctGCGGTAGTGCCCATAATCAGGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCAGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTCGCCAAATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGAAgg                                                                                                                                <  2:342467/150‑1 (MQ=255)
                       ctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAgg                                                                                                                                <  1:528689/150‑1 (MQ=255)
                             tAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCCgg                                                                                                                          <  1:2407932/150‑1 (MQ=255)
                             tAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAgg                                                                                                                          <  2:1986099/150‑1 (MQ=255)
                             tAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAgg                                                                                                                          <  1:1025096/150‑1 (MQ=255)
                                  cccATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCacaa                                                                                                                     <  1:1162849/150‑1 (MQ=255)
                                    cATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCacaaca                                                                                                                   >  1:2286708/1‑150 (MQ=255)
                                    cATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCacaaca                                                                                                                   >  1:392358/1‑150 (MQ=255)
                                    cATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCacaaca                                                                                                                   <  1:639339/150‑1 (MQ=255)
                                     aTAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAg                                                                                                                  >  1:146590/1‑150 (MQ=255)
                                         tCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACgg                                                                                                              <  2:2286708/150‑1 (MQ=255)
                                           aCGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGCCgccg                                                                                                            <  1:1289505/150‑4 (MQ=255)
                                           aCGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTa                                                                                                            >  2:216145/1‑150 (MQ=255)
                                           aCGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTa                                                                                                            <  2:1338850/150‑1 (MQ=255)
                                             gAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCTGGCTGCCAGGGCag                                                                                                                                           <  1:707456/117‑2 (MQ=255)
                                             gAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCTGGCTGCCAGGGCag                                                                                                                                           >  2:707456/1‑116 (MQ=255)
                                               gCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAAAAGACGGTACGGc                                                                                                        <  2:1535313/150‑1 (MQ=255)
                                                cGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCa                                                                                                       <  1:1665213/150‑1 (MQ=255)
                                                      cAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTAtt                                                                                                 <  2:141540/150‑1 (MQ=255)
                                                         ggTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTtct                                                                                              >  1:1316125/1‑150 (MQ=255)
                                                           tAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTcc                                                                                            <  2:146590/150‑1 (MQ=255)
                                                           tAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTcc                                                                                            >  1:2136384/1‑150 (MQ=255)
                                                           tAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGAAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTcc                                                                                            >  1:1600560/1‑150 (MQ=255)
                                                                  tGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCcggcgcg                                                                                     >  2:1465969/1‑150 (MQ=255)
                                                                  tGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCcggcgcg                                                                                     >  2:629397/1‑150 (MQ=255)
                                                                      gTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCcggcgcggc                                                                                   >  2:448186/1‑148 (MQ=255)
                                                                      gTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCcggcgcggc                                                                                   <  1:448186/148‑1 (MQ=255)
                                                                         gAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAgg                                                                                                                          >  2:118216/1‑106 (MQ=255)
                                                                         gAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAgg                                                                                                                          <  1:118216/106‑1 (MQ=255)
                                                                          aaGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGCGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTAGGGCATGTATTTCTCCCGGCGCGGCGGATTa                                                                             <  2:1828930/150‑1 (MQ=255)
                                                                            ggtggtGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGc                                                                           <  2:629711/150‑1 (MQ=255)
                                                                              tggtggCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCgg                                                                         >  2:673973/1‑150 (MQ=255)
                                                                               ggtggCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCggg                                                                        <  1:216145/150‑1 (MQ=255)
                                                                                 tggCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCggggg                                                                      >  2:219290/1‑150 (MQ=255)
                                                                                     cgccTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGggcggc                                                                  <  1:1105628/150‑1 (MQ=255)
                                                                                      gccTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTa                                                                                                            >  2:992179/1‑107 (MQ=255)
                                                                                      gccTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTa                                                                                                            <  1:992179/107‑1 (MQ=255)
                                                                                            gCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCa                                                           <  2:1005795/150‑1 (MQ=255)
                                                                                              gggggCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGt                                                         <  2:2154728/150‑1 (MQ=255)
                                                                                                  gCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAg                                                     <  1:932401/150‑1 (MQ=255)
                                                                                                       tggtAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAgg                                                                                                                                               >  2:492790/1‑55 (MQ=255)
                                                                                                       tggtAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAgg                                                                                                                                               <  1:492790/55‑1 (MQ=255)
                                                                                                        ggtAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCg                                               <  2:1316125/150‑1 (MQ=255)
                                                                                                         gtAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCgg                                              >  1:651777/1‑150 (MQ=255)
                                                                                                                cgaggccctgGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGCTTAGGGGGGGCGGGAAAGCCAGTGCAGGTAGCGGCCAAGcc                                       >  2:570992/8‑150 (MQ=255)
                                                                                                                  ggcTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCg                                     <  1:673973/150‑1 (MQ=255)
                                                                                                                   gctgctGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCgg                                    <  1:317152/150‑1 (MQ=255)
                                                                                                                    ctgctgGCCACATACGGGTTGATGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCACGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGc                                   >  1:1520596/1‑150 (MQ=255)
                                                                                                                     tgctgGCCACATACGGGTTGCTGCGCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCg                                  >  2:794173/1‑150 (MQ=255)
                                                                                                                       ctgGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCATTGCAGGTAGCGGCCAAGCCCGGCGaa                                <  1:570992/150‑1 (MQ=255)
                                                                                                                        tgGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAg                               <  2:464684/150‑1 (MQ=255)
                                                                                                                          gcgaactACGGGTTCTGGGCCCGGCAGGCAGCCGGGGCATTGGGGAGCAGGTGCTGGCACAGCAGAGGGTACGGCATGTATTTCCCGCGGCGCGGCGGATTAGCGGGGGCGGCTAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggt                             <  2:859352/144‑1 (MQ=255)
                                                                                                                            cacaTACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtgg                           >  1:218202/1‑150 (MQ=255)
                                                                                                                             acaTACGGGTTTTTGCCTTGGCAGCCAGCCAGGGCAATGGGCACCAGGGCCAGGCAGAACAGCCGGTACGGCATGTATTCCTCCCGGGGCGGTTGATTAGCGGGGGCGCCGATCCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtggg                          <  2:329267/150‑1 (MQ=255)
                                                                                                                                  cGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGc                     <  2:2352268/150‑1 (MQ=255)
                                                                                                                                   gggTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCg                    >  1:2231284/1‑150 (MQ=255)
                                                                                                                                        attgcCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGggcgg                                                                   <  2:1274888/96‑1 (MQ=255)
                                                                                                                                        attgcCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGggcgg                                                                   >  1:1274888/3‑98 (MQ=255)
                                                                                                                                         ttgcCCTGGCTGGCTGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCa                                                      >  1:2165129/2‑110 (MQ=255)
                                                                                                                                         ttgcCCTGGCTGGCTGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCa                                                      <  2:2165129/109‑1 (MQ=255)
                                                                                                                                         ttgcCCTGGCTGGCTGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTg              >  2:810351/2‑150 (MQ=255)
                                                                                                                                          tgcCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGTCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTgg             <  1:1625008/150‑1 (MQ=255)
                                                                                                                                             cttGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGcc          >  2:602917/3‑150 (MQ=255)
                                                                                                                                             ccTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGcc          >  1:1630849/1‑150 (MQ=255)
                                                                                                                                              cTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtgggg                         <  1:1765653/134‑1 (MQ=255)
                                                                                                                                              cTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtgggg                         >  2:1765653/1‑134 (MQ=255)
                                                                                                                                               tGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAg        >  2:1959907/1‑150 (MQ=255)
                                                                                                                                                   agccagccAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAATCCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTcc    <  1:794173/150‑1 (MQ=255)
                                                                                                                                                    gccagccaGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTCCa   >  2:1363401/1‑150 (MQ=255)
                                                                                                                                                     ccagccagGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTCCAt  <  1:629397/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTGGCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTCCAT  >  NC_002947/675443‑675741

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