Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,185,776:1 +G 100% intergenic (‑44/‑127) PP_1033 ← / → cumA sulfatase domain‑containing protein/multicopper oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,185,7761.G78.2% 816.2 / 249.5 294intergenic (‑44/‑127)PP_1033/cumAsulfatase domain‑containing protein/multicopper oxidase
Reads supporting (aligned to +/- strand):  ref base . (28/36);  new base G (127/102);  total (156/138)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 319 total aligned reads displayed.

ACATGGTGTCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTGCGCAGAAACGAGGTA  >  NC_002947/1185627‑1185920
                                                                                                                                                      |                                                                                                                                                
 caTGGTGTCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAg                                                                                                                                                 <  2:3268268/150‑1 (MQ=255)
  aTGGTGTCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAgg                                                                                                                                                >  1:3034112/1‑150 (MQ=255)
  aTGGTGTCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAgg                                                                                                                                                <  1:4753542/150‑1 (MQ=255)
   tGGTGTCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAggg                                                                                                                                               >  2:4855704/1‑150 (MQ=255)
   tGGTGTCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAggg                                                                                                                                               <  1:249386/150‑1 (MQ=255)
     gtgtCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAggggg                                                                                                                                             <  1:2906172/150‑1 (MQ=255)
     gtgtCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAggggg                                                                                                                                             <  1:4987127/150‑1 (MQ=255)
     gtgtCTACCGTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAggggg                                                                                                                                             <  1:5738575/150‑1 (MQ=255)
      tgtCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAgggggg                                                                                                                                            >  2:1811527/1‑150 (MQ=255)
        tCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAgggggggg                                                                                                                                          >  2:133394/1‑149 (MQ=255)
          tACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAggggggggcg                                                                                                                                        >  1:3869252/1‑147 (MQ=255)
             cTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGAt                                                                                                                                      >  2:3075969/1‑150 (MQ=255)
              ttCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGAtt                                                                                                                                     >  1:1419649/1‑150 (MQ=255)
              ttCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGAtt                                                                                                                                     <  1:4660301/150‑1 (MQ=255)
              ttCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGAtt                                                                                                                                     >  2:601218/1‑150 (MQ=255)
               tCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTg                                                                                                                                    >  2:5231180/1‑150 (MQ=255)
               tCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGc                                                                                                                                   <  2:1559524/150‑1 (MQ=255)
                 gCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGcc                                                                                                                                  <  2:877954/150‑1 (MQ=255)
                  cATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCa                                                                                                                                 >  2:5315889/1‑150 (MQ=255)
                     tCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTgg                                                                                                                             <  1:3122240/150‑1 (MQ=255)
                      cacaGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTgg                                                                                                                             <  1:2182529/150‑1 (MQ=255)
                       acaGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGc                                                                                                                            >  1:5822172/1‑150 (MQ=255)
                       acaGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGc                                                                                                                            >  2:6476180/1‑150 (MQ=255)
                           agCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATc                                                                                                                        <  2:2562094/150‑1 (MQ=255)
                              cGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTGTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATcgcg                                                                                                                     >  1:1202872/1‑150 (MQ=255)
                                gCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGc                                                                                                                   >  1:2024062/1‑150 (MQ=255)
                                     ttGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTa                                                                                                              >  2:6649269/1‑150 (MQ=255)
                                          tGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCCAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTACCCgg                                                                                                         >  1:1378054/1‑150 (MQ=255)
                                          tGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCgg                                                                                                         >  2:6624108/1‑150 (MQ=255)
                                            aaGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGcc                                                                                                       >  2:5197076/1‑150 (MQ=255)
                                              gCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGt                                                                                                     >  1:1945159/1‑150 (MQ=255)
                                                gATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGGTTAGCGGCCgcga                                                                                                   <  1:3985579/150‑4 (MQ=255)
                                                gATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGc                                                                                                   <  2:479687/150‑1 (MQ=255)
                                                gATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGc                                                                                                   <  2:1440959/150‑1 (MQ=255)
                                                    ggCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACaa                                                                                              <  2:3692781/150‑1 (MQ=255)
                                                    ggCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACaa                                                                                              >  1:2696135/1‑150 (MQ=255)
                                                         gCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACt                                                                                          >  1:3737353/1‑150 (MQ=255)
                                                          cAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACtt                                                                                         >  2:2190172/1‑150 (MQ=255)
                                                           aaaaGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCACGTCCGCCCCGGGGCGGTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACCAGACTTc                                                                                        >  1:649671/1‑150 (MQ=255)
                                                           aaaaGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTc                                                                                        >  1:1492084/1‑150 (MQ=255)
                                                           aaaaGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTc                                                                                        <  2:6222806/150‑1 (MQ=255)
                                                           aaaaGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGCACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTc                                                                                        >  2:3576106/1‑150 (MQ=255)
                                                             aaGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGTCTTCAc                                                                                      <  1:928335/150‑1 (MQ=255)
                                                             aaGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGCACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCAc                                                                                      <  1:3576106/150‑1 (MQ=255)
                                                                 ccAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTaa                                                                                  <  1:179559/150‑1 (MQ=255)
                                                                  cAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAat                                                                                 <  2:734979/150‑1 (MQ=255)
                                                                     cACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATAtt                                                                              >  1:562187/1‑150 (MQ=255)
                                                                     cACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATAtt                                                                              <  1:4571996/150‑1 (MQ=255)
                                                                       ccTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTc                                                                            <  2:4105798/150‑1 (MQ=255)
                                                                         tGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACATAATATTTCGc                                                                          >  1:6127361/1‑150 (MQ=255)
                                                                            gcACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCaaa                                                                       >  2:3383311/1‑150 (MQ=255)
                                                                                cccGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCTCTTAATATTTCGCAAATGCa                                                                   >  2:323013/1‑150 (MQ=255)
                                                                                cccGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCa                                                                   >  1:6724123/1‑150 (MQ=255)
                                                                                cccGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCa                                                                   >  1:1654659/1‑150 (MQ=255)
                                                                                cccGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCa                                                                   >  2:2948608/1‑150 (MQ=255)
                                                                                cccGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCa                                                                   >  2:5439747/1‑150 (MQ=255)
                                                                                cccGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTAGCAAATGCa                                                                   >  1:6386381/1‑150 (MQ=255)
                                                                                cccGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAg                                                                  <  2:3825269/150‑1 (MQ=255)
                                                                                cccGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAg                                                                  <  2:1911882/150‑1 (MQ=255)
                                                                                 ccGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAg                                                                  <  1:5489887/150‑1 (MQ=255)
                                                                                   gTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGgt                                                                >  2:5678228/1‑150 (MQ=255)
                                                                                   gTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGgt                                                                <  2:2015435/150‑1 (MQ=255)
                                                                                   gTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGgt                                                                >  2:6310562/1‑150 (MQ=255)
                                                                                       ttAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGGAt                                                            >  2:2769899/1‑150 (MQ=255)
                                                                                       ttAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATc                                                           <  2:1963566/150‑1 (MQ=255)
                                                                                             ttAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGGCTTCACTTAATATTTCGCAAATGCAGGTGTATCGTAtt                                                      >  2:4296437/1‑150 (MQ=255)
                                                                                             ttAAGCAGCCCACCCTGAATGCGGCACGTCAGCCCCGGGGCGTTGGGTGGGTCCCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTAtg                                                      >  2:1054233/1‑149 (MQ=255)
                                                                                              tAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCg                                                    <  1:1739168/150‑1 (MQ=255)
                                                                                               aaGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCg                                                    <  2:654894/150‑1 (MQ=255)
                                                                                                agcagcCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCACAAATGCAGGTGTATCGTATTAGc                                                   >  2:6783293/1‑150 (MQ=255)
                                                                                                 gcagcCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCCGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGcc                                                  >  1:5299606/1‑150 (MQ=255)
                                                                                                  cagcCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGccc                                                 <  2:2591146/150‑1 (MQ=255)
                                                                                                     cccACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTc                                              >  2:2484627/1‑150 (MQ=255)
                                                                                                      ccACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCAt                                            >  1:709897/1‑150 (MQ=255)
                                                                                                       cACCCTGAATTCGGCAAGTCTGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTTCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATc                                           <  2:1086127/150‑1 (MQ=255)
                                                                                                           cTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTCACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGcc                                        >  2:4127743/1‑150 (MQ=255)
                                                                                                            tGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCa                                       <  2:1673938/150‑1 (MQ=255)
                                                                                                            tGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCa                                       <  2:2107011/150‑1 (MQ=255)
                                                                                                            tGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCa                                       <  1:4021937/150‑1 (MQ=255)
                                                                                                            tGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTGCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCa                                       >  1:4819757/1‑150 (MQ=255)
                                                                                                               aTTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGcc                                    >  1:3438828/1‑150 (MQ=255)
                                                                                                                    gCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCctgc                               >  2:3872702/1‑150 (MQ=255)
                                                                                                                    gCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCctgc                               >  2:5073255/1‑150 (MQ=255)
                                                                                                                    gCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCctgct                              >  2:3453442/1‑150 (MQ=255)
                                                                                                                      aaGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCctgctg                             >  1:2785987/1‑150 (MQ=255)
                                                                                                                      aaGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCctgctg                             >  1:1080458/1‑150 (MQ=255)
                                                                                                                       aGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTgg                            <  2:5765453/150‑1 (MQ=255)
                                                                                                                       aGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGc                           <  2:5506725/150‑1 (MQ=255)
                                                                                                                                cGGGGCGTTGGTTGGGTACCGTGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTGAATATATCGCAAATGCAGGTGTATCGTATTCGCCCGTCATCGCCAGCCCCTGCTGGCGGACGGAtt                  >  1:105402/1‑150 (MQ=255)
                                                                                                                                cGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGAt                   >  2:4368775/1‑150 (MQ=255)
                                                                                                                                cGGGGCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGAt                   >  1:3943889/1‑150 (MQ=255)
                                                                                                                                 ggggCGTTGGTTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGAtt                  <  1:5231180/150‑1 (MQ=255)
                                                                                                                                     cGTTGGGTGGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTgcgc              <  1:4117797/150‑1 (MQ=255)
                                                                                                                                      gttggttgGGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTGCGCa             >  2:2902203/1‑150 (MQ=255)
                                                                                                                                         ggttggGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTGCGCAGaa          >  2:4300632/1‑150 (MQ=255)
                                                                                                                                         ggttggGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTGCGCAGaa          >  2:3986187/1‑150 (MQ=255)
                                                                                                                                          gttggGTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTGCGCAGaaa         >  1:2600123/1‑150 (MQ=255)
                                                                                                                                             gggTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTGCGCAGAAACGa      <  1:6310562/150‑1 (MQ=255)
                                                                                                                                               gTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAGTATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTGCGCAGAAACGAgg    >  1:3222621/1‑150 (MQ=255)
                                                                                                                                               gTACCGAGGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTGCGCAGAAACGAgg    >  2:4610724/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
ACATGGTGTCTACCTTCGCATTCACAGAGCCGGCGTCTTGTCTGAAGCGATAGGCGGGCAAAAGGCCAGCACCTGCGCACCCCGTTTTTAGACTTAAGCAGCCCACCCTGAATTCGGCAAGTCAGCCCCGGGGCGTTGGTTGGGTACCGAGGGGGGGCGGATTGCCACTGGCGATCGCGGCCGCTAACCGGCCGTGCGACAAGACTTCACTTAATATTTCGCAAATGCAGGTGTATCGTATTCGCCCTTCATCGCCAGCCCCTGCTGGCAGACTGATTGCGCAGAAACGAGGTA  >  NC_002947/1185627‑1185920

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