Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,921,589 C→G 5.2% A666A (GCC→GCG)  PP_4330 → toxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,921,5890CG 5.2% 252.6 / 3.3 76A666A (GCC→GCG) PP_4330toxin‑like protein
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base G (2/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

GCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACC  >  NC_002947/4921441‑4921737
                                                                                                                                                      |                                                                                                                                                    
gCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGcc                                                                                                                                                      >  1:197765/1‑149 (MQ=255)
  gACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTCGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCATACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGccc                                                                                                                                                     >  1:318453/1‑149 (MQ=255)
     aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAgg                                                                                                                                                 >  4:349103/1‑149 (MQ=255)
           aTTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTGcc                                                                                                                                           <  6:77606/149‑2 (MQ=255)
            ttAGCGAGCCATTGCTGGCCAACCTTGCCGCCACCGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAg                                                                                                                                          <  4:331118/149‑1 (MQ=255)
                cGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTAt                                                                                                                                      <  7:289562/149‑1 (MQ=255)
                      aTTCCTGGCCACCCTTGCCTCCACCGCCGCGCTCCT‑GAGCGACAGTGACTACAGCATGCGCCTGTTGCTGTCGAAGT‑CGTCGCCTACCGCTAACGCACAGACCGCAAGGCTATTGGCAGCCCGGGCCCCCGGGCGGCAGGTATTGCCGc                                                                                                                                <  6:222433/149‑1 (MQ=255)
                            ggCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCCACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCCGCCCCCAGGGCTGCACGTATTGCCCCTTg                                                                                                                             >  2:322600/1‑146 (MQ=255)
                             gCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGaa                                                                                                                           >  1:195597/1‑147 (MQ=255)
                             gCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGaa                                                                                                                           <  2:195597/147‑1 (MQ=255)
                             gCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAag                                                                                                                          <  5:269999/148‑1 (MQ=255)
                             gCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAACATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGaa                                                                                                                           <  1:321960/147‑1 (MQ=255)
                               cAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGc                                                                                                                       <  7:193185/149‑1 (MQ=255)
                                 aCCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCaa                                                                                                                     <  2:197765/149‑1 (MQ=255)
                                  ccTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAg                                                                                                                    <  5:249617/149‑1 (MQ=255)
                                   cTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCt                                                                                                                               <  4:405714/137‑1 (MQ=255)
                                   cTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCt                                                                                                                               >  3:405714/1‑137 (MQ=255)
                                    ttGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGc                                                                                                                  <  4:221451/149‑1 (MQ=255)
                                        cgccAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGGAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttct                                                                                                              >  6:118449/1‑149 (MQ=255)
                                        cgccAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttct                                                                                                              <  3:349103/149‑1 (MQ=255)
                                        cgccAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttct                                                                                                              >  3:111936/1‑149 (MQ=255)
                                          ccAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCAAGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCttcttc                                                                                                            <  3:114938/149‑1 (MQ=255)
                                            aaCGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCg                                                                                                          <  1:260706/149‑1 (MQ=255)
                                            aaCGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCg                                                                                                          <  8:80753/149‑1 (MQ=255)
                                               gccgcGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGc                                                                                                       <  7:150973/149‑1 (MQ=255)
                                                    gcgcCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCt                                                                                                  <  2:262449/149‑1 (MQ=255)
                                                    gcgcCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCt                                                                                                  <  2:209489/149‑1 (MQ=255)
                                                        cT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGc                                                                                                >  1:259721/1‑147 (MQ=255)
                                                              cGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGcc                                                                                         >  8:136889/1‑149 (MQ=255)
                                                              cGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGcc                                                                                         >  3:134722/1‑149 (MQ=255)
                                                                    tGAATACAGCATGCGCCTGTTTCTGTCGAAGGACTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTATTCCGGGCCTTCTGACAGAGCCGCTAt                                                                                    >  1:256079/1‑149 (MQ=255)
                                                                      aaTACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTAc                                                                                 >  8:47586/1‑149 (MQ=255)
                                                                           aGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttct                                                                                                              >  5:144242/1‑115 (MQ=255)
                                                                           aGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttct                                                                                                              <  6:144242/115‑1 (MQ=255)
                                                                           aGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTAcgcg                                                                              >  8:338168/1‑147 (MQ=255)
                                                                                 cgcCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGa                                                                      <  7:338168/149‑1 (MQ=255)
                                                                                  gcCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAggg                                                                                                                                                >  3:268629/1‑74 (MQ=255)
                                                                                  gcCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAggg                                                                                                                                                <  4:268629/74‑1 (MQ=255)
                                                                                  gcCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                     >  4:366637/1‑149 (MQ=255)
                                                                                   ccTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAgg                                                                    >  1:155192/1‑149 (MQ=255)
                                                                                          cTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCt                                                             >  8:294138/1‑149 (MQ=255)
                                                                                            gTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTgg                                                           >  7:134100/1‑149 (MQ=255)
                                                                                               gAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCCGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATg                                                        >  5:237813/1‑149 (MQ=255)
                                                                                               gAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGa                                                                      <  1:266574/135‑1 (MQ=255)
                                                                                               gAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGa                                                                      >  2:266574/1‑135 (MQ=255)
                                                                                               gAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGAt                                                         <  7:47586/148‑1 (MQ=255)
                                                                                               gAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGAt                                                         <  7:83594/148‑1 (MQ=255)
                                                                                               gAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATg                                                        >  7:341105/1‑149 (MQ=255)
                                                                                                aaGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATgg                                                       >  1:272394/1‑149 (MQ=255)
                                                                                                     cTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAAc                                                   >  8:87328/1‑149 (MQ=255)
                                                                                                         gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTg                                               <  2:256079/149‑1 (MQ=255)
                                                                                                              ccGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATc                                          <  2:155192/149‑1 (MQ=255)
                                                                                                                    aCGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGccac                                     <  2:272394/148‑1 (MQ=255)
                                                                                                                      gAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCaa                                  >  2:389047/1‑149 (MQ=255)
                                                                                                                      gAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCaa                                  >  1:47305/1‑149 (MQ=255)
                                                                                                                       aaCAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAAt                                 >  3:396488/1‑149 (MQ=255)
                                                                                                                        aCAGAACGCAAGGCTATTGGCAGCCCGGGCGCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAAt                                 >  1:35649/1‑148 (MQ=255)
                                                                                                                          aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATg                                                        >  5:213017/1‑123 (MQ=255)
                                                                                                                          aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATg                                                        <  6:213017/123‑1 (MQ=255)
                                                                                                                          aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTg                                               <  4:169344/132‑1 (MQ=255)
                                                                                                                          aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTg                                               >  3:169344/1‑132 (MQ=255)
                                                                                                                            aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                            >  8:176622/1‑149 (MQ=255)
                                                                                                                            aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                            <  4:111936/149‑1 (MQ=255)
                                                                                                                               gCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTc                         >  1:374113/1‑149 (MQ=255)
                                                                                                                               gCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTc                         >  5:111011/1‑149 (MQ=255)
                                                                                                                                cAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                        <  7:294138/149‑1 (MQ=255)
                                                                                                                                 aaGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAgcgc                                                              <  6:261705/110‑1 (MQ=255)
                                                                                                                                 aaGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAgcgc                                                              >  5:261705/1‑110 (MQ=255)
                                                                                                                                 aaGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                                >  4:404004/1‑124 (MQ=255)
                                                                                                                                 aaGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                                <  3:404004/124‑1 (MQ=255)
                                                                                                                                  aGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGc                      <  7:87328/149‑1 (MQ=255)
                                                                                                                                      tATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGAcc                  >  8:29644/1‑149 (MQ=255)
                                                                                                                                        ttGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCg                 <  7:136889/148‑1 (MQ=255)
                                                                                                                                         ctgCAGCCCTGGGGCCCGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATc                                          <  7:232984/120‑1 (MQ=255)
                                                                                                                                         ctgCAGCCCTGGGGCCCGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATc                                          >  8:232984/3‑122 (MQ=255)
                                                                                                                                                ccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgctgc        >  6:197923/1‑149 (MQ=255)
                                                                                                                                                 cGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTgcg       <  5:118449/149‑1 (MQ=255)
                                                                                                                                                    gCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTa    <  2:374113/149‑1 (MQ=255)
                                                                                                                                                      cccAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTAcc  >  1:97009/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCT‑GAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGG‑CTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACC  >  NC_002947/4921441‑4921737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑