Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 675,594 A→T 8.4% A36A (GCT→GCA)  PP_0577 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947675,5940AT 8.4% 249.9 / 11.7 85A36A (GCT→GCA) PP_0577hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/42);  new base T (3/4);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

GCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTCCAT  >  NC_002947/675447‑675741
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                          <  1:1025096/150‑1 (MQ=255)
                         tAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAgg                                                                                                                          <  2:1986099/150‑1 (MQ=255)
                              cccATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCacaa                                                                                                                     <  1:1162849/150‑1 (MQ=255)
                                cATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCacaaca                                                                                                                   >  1:392358/1‑150 (MQ=255)
                                cATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCacaaca                                                                                                                   >  1:2286708/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:1338850/150‑1 (MQ=255)
                                       aCGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTa                                                                                                            >  2:216145/1‑150 (MQ=255)
                                         gAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCTGGCTGCCAGGGCag                                                                                                                                           >  2:707456/1‑116 (MQ=255)
                                         gAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCTGGCTGCCAGGGCag                                                                                                                                           <  1:707456/117‑2 (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                                                                                   <  1:448186/148‑1 (MQ=255)
                                                                  gTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCcggcgcggc                                                                                   >  2:448186/1‑148 (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)
                                                                                                                                                   |                                                                                                                                                   
GCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTCCAT  >  NC_002947/675447‑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: ‑