Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 134,397 A→C 5.0% D98A (GAT→GCT) ‡ PP_0129 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947134,3970AC 5.0% 611.2 / 13.6 178D98A (GAT→GCT) ‡PP_0129GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (86/83);  new base C (5/4);  total (91/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 178 total aligned reads displayed.

AAGGTGGTAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCGGCT  >  NC_002947/134249‑134543
                                                                                                                                                    |                                                                                                                                                  
aaGGTGGTAGACGAGTGCATGAAGGAACTGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGAGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGAt                                                                                                                                                   <  2:2965432/150‑1 (MQ=255)
aaGGTGGTAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGAt                                                                                                                                                   >  1:3324826/1‑150 (MQ=255)
  ggtggtAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATgc                                                                                                                                                 >  1:748448/1‑150 (MQ=255)
     ggtAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCt                                                                                                                                              <  1:4551230/150‑1 (MQ=255)
     ggtAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTGGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCt                                                                                                                                              >  1:2711608/1‑150 (MQ=255)
      gtAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTg                                                                                                                                             >  2:3516118/1‑150 (MQ=255)
        aGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGAc                                                                                                                                           >  2:2149185/1‑150 (MQ=255)
           cGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGCATCGCTGAcctg                                                                                                                                        <  2:4344601/150‑3 (MQ=255)
             aGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCAc                                                                                                                                      <  1:1240220/150‑1 (MQ=255)
                gCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTgg                                                                                                                                   >  2:4229233/1‑150 (MQ=255)
                gCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTgg                                                                                                                                   <  1:4714796/150‑1 (MQ=255)
                gCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTgg                                                                                                                                   <  2:949262/150‑1 (MQ=255)
                  aTGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGtt                                                                                                                                 >  2:1581055/1‑150 (MQ=255)
                  aTGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGtt                                                                                                                                 <  1:3738325/150‑1 (MQ=255)
                   tGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTg                                                                                                                                >  2:3915395/1‑150 (MQ=255)
                   tGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTg                                                                                                                                >  2:1368180/1‑150 (MQ=255)
                   tGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTg                                                                                                                                <  1:4220151/150‑1 (MQ=255)
                    gtaGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCTCTGGTTgg                                                                                                                               <  2:1036521/148‑1 (MQ=255)
                    gAATGAAATGCGCGATGTCATTCGCAGCGACAACATGGTCGCCGGCCTCGCTGGCTTGATTCCGCTCCTTTATATAGCGGTACTGGATTCCGAACAGCGTCGCGAATCCCGCATGATCCAGGTCATCGATGCGCTGACCGCACTGGTTgg                                                                                                                               <  1:4004675/150‑1 (MQ=255)
                              cgcgATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCaa                                                                                                                     <  1:1100975/150‑1 (MQ=255)
                                     tCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgc                                                                                                              <  1:2515186/150‑1 (MQ=255)
                                     tCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgc                                                                                                              >  1:380579/1‑150 (MQ=255)
                                      cATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgcc                                                                                                             <  1:1740424/150‑1 (MQ=255)
                                        ttCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgcccc                                                                                                           >  2:318924/1‑150 (MQ=255)
                                           tcAGAGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACAGGATGGCCAACTGCAAGGcctgcccctgc                                                                                                        >  2:3795573/2‑150 (MQ=255)
                                           gCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgcccctgc                                                                                                        >  2:512558/1‑150 (MQ=255)
                                           gCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgcccctgc                                                                                                        >  2:4064863/1‑150 (MQ=255)
                                           gCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgcccctgc                                                                                                        <  2:1395531/150‑1 (MQ=255)
                                                gACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGc                                                                                                   >  1:4355459/1‑150 (MQ=255)
                                                gACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGc                                                                                                   <  1:2286264/150‑1 (MQ=255)
                                                   aacaTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGat                                                                                                >  1:4104709/1‑150 (MQ=255)
                                                   aacaTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGat                                                                                                >  1:3439612/1‑150 (MQ=255)
                                                      aTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATAtc                                                                                             >  2:1567113/1‑150 (MQ=255)
                                                       tGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATAtct                                                                                            >  2:2512254/1‑150 (MQ=255)
                                                       tGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATAtct                                                                                            >  1:4628802/1‑150 (MQ=255)
                                                          aCGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTcgc                                                                                         <  2:5047596/150‑1 (MQ=255)
                                                            gCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCgg                                                                                       <  2:770803/150‑1 (MQ=255)
                                                               ggCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCa                                                                                    >  1:2498552/1‑150 (MQ=255)
                                                                   tCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGa                                                                                >  1:298998/1‑150 (MQ=255)
                                                                   tCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGa                                                                                <  2:873045/150‑1 (MQ=255)
                                                                   tCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGa                                                                                >  1:219781/1‑150 (MQ=255)
                                                                       tGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGaaaaa                                                                            >  1:1332689/1‑150 (MQ=255)
                                                                       tGGCTTGATTCCGCGCCTGGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGaaaaa                                                                            >  1:4807329/1‑150 (MQ=255)
                                                                        ggCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGaaaaaa                                                                           <  2:1080870/150‑1 (MQ=255)
                                                                           ttGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTg                                                                        >  1:471874/1‑150 (MQ=255)
                                                                           ttGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTg                                                                        >  1:1418793/1‑150 (MQ=255)
                                                                           ttGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTg                                                                        >  2:3973262/1‑150 (MQ=255)
                                                                           ttGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTg                                                                        >  1:1901818/1‑150 (MQ=255)
                                                                           ttGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTg                                                                        >  2:1480089/1‑150 (MQ=255)
                                                                           ttGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTg                                                                        >  1:4220628/1‑150 (MQ=255)
                                                                           ttGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTg                                                                        >  2:881086/1‑150 (MQ=255)
                                                                             gATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGc                                                                      <  1:555917/150‑1 (MQ=255)
                                                                               ttCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGctc                                                                    >  2:49574/1‑148 (MQ=255)
                                                                               ttCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCga                                                                    <  1:154975/150‑1 (MQ=255)
                                                                                  cgcgCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCgaaga                                                                 >  1:2532838/1‑150 (MQ=255)
                                                                                   gcgcCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCgaagaa                                                                >  1:507031/1‑150 (MQ=255)
                                                                                     gcCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGc                                                              >  2:437951/1‑150 (MQ=255)
                                                                                       cTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTc                                                            <  1:1189740/150‑1 (MQ=255)
                                                                                       cTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTc                                                            <  1:2352090/150‑1 (MQ=255)
                                                                                       cTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTc                                                            <  1:2856537/150‑1 (MQ=255)
                                                                                         tGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGa                                                          >  1:4367957/1‑150 (MQ=255)
                                                                                           agaAAGCGGTACTGGATTCCGAACATCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGAACTGGTTGGCCAACTGCTATGCCTGCCCCTGACCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAgg                                                        <  1:3517865/150‑1 (MQ=255)
                                                                                             aaaGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCgcgtga                                                      >  2:683207/1‑145 (MQ=255)
                                                                                               aGCGGTACTGTATTCCGAACTGCGTCGCGAAACCCGCATGAACCCGGTCAGCGATGCGCTGACGGCACTGGTTGGCCAACGGAAAGGCCGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCggg                                                   <  2:3866078/150‑1 (MQ=255)
                                                                                               aGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCgg                                                    >  2:4985088/1‑150 (MQ=255)
                                                                                                gCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCggg                                                   <  1:2410587/150‑1 (MQ=255)
                                                                                                gCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCggg                                                   <  1:2149185/150‑1 (MQ=255)
                                                                                                 cGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCgggg                                                  <  2:5010312/150‑1 (MQ=255)
                                                                                                 cGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCgggg                                                  <  1:314061/150‑1 (MQ=255)
                                                                                                 cGGTACTGGATTCAGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCgggg                                                  <  2:3282056/150‑1 (MQ=255)
                                                                                                    tACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCCCTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTg                                               >  2:1425793/1‑150 (MQ=255)
                                                                                                    tACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTg                                               <  1:1282796/150‑1 (MQ=255)
                                                                                                     aCTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGc                                              >  2:4763662/1‑150 (MQ=255)
                                                                                                      cTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGcc                                             <  1:4503399/150‑1 (MQ=255)
                                                                                                      cTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGcc                                             <  2:2749317/150‑1 (MQ=255)
                                                                                                        ggATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCa                                           <  2:2393369/150‑1 (MQ=255)
                                                                                                        ggATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCa                                           <  2:3191684/150‑1 (MQ=255)
                                                                                                        ggATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCa                                           >  1:1981444/1‑150 (MQ=255)
                                                                                                        ggATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGCTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCa                                           <  2:1034418/150‑1 (MQ=255)
                                                                                                         gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATGGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAg                                          >  2:2292626/1‑150 (MQ=255)
                                                                                                             ccGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAc                                      <  2:1711810/150‑1 (MQ=255)
                                                                                                             ccGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAc                                      <  2:4373134/150‑1 (MQ=255)
                                                                                                               gAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAcgc                                    >  2:192921/1‑150 (MQ=255)
                                                                                                                 aCAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGa                                  <  2:3554286/150‑1 (MQ=255)
                                                                                                                    gCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCt                               <  2:4547349/150‑1 (MQ=255)
                                                                                                                       tCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTgcc                            <  2:2195460/150‑1 (MQ=255)
                                                                                                                           gAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGc                                                              >  1:722931/1‑112 (MQ=255)
                                                                                                                           gAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGc                                                              <  2:722931/112‑1 (MQ=255)
                                                                                                                               cccGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGt                    <  1:4435850/150‑1 (MQ=255)
                                                                                                                                ccGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGtt                   >  2:2475348/1‑150 (MQ=255)
                                                                                                                                 cGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTg                  <  1:2019210/150‑1 (MQ=255)
                                                                                                                                  gCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTgg                 <  2:4451130/150‑1 (MQ=255)
                                                                                                                                      gAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGa             <  2:223254/150‑1 (MQ=255)
                                                                                                                                         ccAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGtt          <  2:4361872/150‑1 (MQ=255)
                                                                                                                                           cggTCAGCGCATCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGa        >  2:4815109/2‑150 (MQ=255)
                                                                                                                                            ggTCAGCGCATCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAg       >  2:3548365/1‑150 (MQ=255)
                                                                                                                                            ggTCAGCGCATCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAg       >  2:3089777/1‑150 (MQ=255)
                                                                                                                                             gTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGc      <  2:4831546/150‑1 (MQ=255)
                                                                                                                                             gTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGc      <  1:102793/150‑1 (MQ=255)
                                                                                                                                                aGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCGGc   <  1:4957858/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAGGTGGTAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCGGCT  >  NC_002947/134249‑134543

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