Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,084,406 C→A 5.1% S183* (TCG→TAG)  PP_3594 → amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,084,4060CA 5.1% 262.6 / 2.9 82S183* (TCG→TAG) PP_3594amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base A (2/2);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

CTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTACATCGAACGGCGCAGCAACCG  >  NC_002947/4084261‑4084553
                                                                                                                                                 |                                                                                                                                                   
cTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCt                                                                                                                                                  >  8:135750/1‑149 (MQ=255)
   gAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTgg                                                                                                                                                >  5:232138/1‑148 (MQ=255)
   gAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTg                                                                                                                                                 <  3:326116/147‑1 (MQ=255)
   gAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGt                                                                                                                                               <  6:232138/149‑1 (MQ=255)
    aaGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTc                                                                                                                                              >  5:342808/1‑149 (MQ=255)
      gCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCAGCGAGATGGTCGa                                                                                                                                            >  7:183431/1‑149 (MQ=255)
       ccTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCAGCGAGATGGTCGAc                                                                                                                                           <  8:183431/149‑1 (MQ=255)
        cTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGa                                                                                                                                            >  6:139986/1‑147 (MQ=255)
        cTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGa                                                                                                                                            <  5:139986/147‑1 (MQ=255)
              ggCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATgc                                                                                                                                    >  1:71186/1‑149 (MQ=255)
                      ggCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGAc                                                                                                                             >  2:32986/1‑148 (MQ=255)
                       gCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGAcc                                                                                                                            >  7:75929/1‑148 (MQ=255)
                       gCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGAcc                                                                                                                            >  2:244684/1‑148 (MQ=255)
                              aaGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGc                                                                                                                    >  8:284727/1‑149 (MQ=255)
                               aGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCg                                                                                                                   <  7:135750/149‑1 (MQ=255)
                                  gcACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGCCCATCTCGCTGGTCGACATGATGCGCCAGACCAACCTGGCGGt                                                                                                                 >  1:65355/1‑148 (MQ=255)
                                  gcACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCGTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGt                                                                                                                 >  1:113391/1‑148 (MQ=255)
                                   cACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGAt                                                                                                                                      <  3:130848/126‑1 (MQ=255)
                                   cACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGAt                                                                                                                                      >  4:130848/1‑126 (MQ=255)
                                    aCCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCTCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCTCCAGACCACCCTGGCGGTCAg                                                                                                              >  8:375382/1‑149 (MQ=255)
                                    aCCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAg                                                                                                              >  5:229487/1‑149 (MQ=255)
                                        gTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCAGCGAGATGGTCGAcacca                                                                                                                                       >  8:300039/1‑116 (MQ=255)
                                        gTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCAGCGAGATGGTCGAcacca                                                                                                                                       <  7:300039/120‑5 (MQ=255)
                                        gTTTCTTCAAGGTGATGCTGCCACAACTGACCCGCATCGCGCTGCACGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTa                                                                                                          <  2:65355/149‑1 (MQ=255)
                                           tcttcAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAAc                                                                                                        <  8:75929/148‑1 (MQ=255)
                                           tcttcAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGAAGCGCCAGACCAACCTGGCGGTCAGCGTAAc                                                                                                        <  1:32986/148‑1 (MQ=255)
                                               cAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCa                                                                                                               <  6:172974/137‑1 (MQ=255)
                                               cAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCa                                                                                                               >  5:172974/1‑137 (MQ=255)
                                                 aGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTAGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGaa                                                                                                 >  2:214493/1‑149 (MQ=255)
                                                      ctgctgcCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGa                                                                                             <  4:243546/148‑1 (MQ=255)
                                                      ctgctgcCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGAt                                                                                            <  2:276485/149‑1 (MQ=255)
                                                      ctgctgcCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGAt                                                                                            <  1:214493/149‑1 (MQ=255)
                                                       tgctgcCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATg                                                                                           >  3:265181/1‑149 (MQ=255)
                                                            ccACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGttc                                                                                        <  2:71186/147‑1 (MQ=255)
                                                              acaACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCt                                                                                                                                                  <  7:145014/87‑1 (MQ=255)
                                                              acaACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCt                                                                                                                                                  >  8:145014/1‑87 (MQ=255)
                                                              acaACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGttcttc                                                                                     >  4:158302/1‑148 (MQ=255)
                                                                  cTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACgg                                                                                >  1:134025/1‑149 (MQ=255)
                                                                   tGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGc                                                                                                                                                   >  8:188996/1‑81 (MQ=255)
                                                                   tGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGc                                                                                                                                                   <  7:188996/81‑1 (MQ=255)
                                                                       ccGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCg                                                                           >  6:116789/1‑149 (MQ=255)
                                                                            tCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCgcctgc                                                                      <  3:158302/149‑1 (MQ=255)
                                                                               cccTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCt                                                                                                                                                  <  8:383514/70‑1 (MQ=255)
                                                                               cccTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCt                                                                                                                                                  >  7:383514/1‑70 (MQ=255)
                                                                               cccTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGcctgcctg                                                                   >  3:403885/1‑149 (MQ=255)
                                                                                ccTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTgg                                                                  <  5:116789/149‑1 (MQ=255)
                                                                                 cTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTggg                                                                 <  1:244684/149‑1 (MQ=255)
                                                                                  tGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCgcctgc                                                                      >  7:92415/1‑143 (MQ=255)
                                                                                  tGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCgcctgc                                                                      <  8:92415/143‑1 (MQ=255)
                                                                                  tGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCACACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTATACGGCGTCGCCTTCATGGGc                                                                >  6:133615/1‑149 (MQ=255)
                                                                                    cccGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCgcctgcc                                                                     <  3:161599/142‑1 (MQ=255)
                                                                                    cccGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCgcctgcc                                                                     >  4:161599/1‑142 (MQ=255)
                                                                                     ccGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTAc                                                             <  5:133615/149‑1 (MQ=255)
                                                                                        gCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACctg                                                          <  2:355524/149‑1 (MQ=255)
                                                                                            gTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGt                                                      >  7:211594/1‑149 (MQ=255)
                                                                                               caacaaCTGGCTGTCGCTGCTCCAGGACACCTCGCTGGTGTCGACCATCTCGCTGGGCGACCTGATGCGCCAGCCCAACCTGGCGGTCAGCGTAACCAAgg                                                                                                   <  5:90198/101‑1 (MQ=255)
                                                                                               caacaaCTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTaa                                                                                                         >  5:160389/1‑95 (MQ=255)
                                                                                               caacaaCTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTaa                                                                                                         <  6:160389/95‑1 (MQ=255)
                                                                                               caacaaCTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTGTTCTACGGCGTCGcctgcctg                                                                   <  8:284473/133‑1 (MQ=255)
                                                                                               caacaaCTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACg                                                                                 <  4:5074/119‑1 (MQ=255)
                                                                                               caacaaCTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACg                                                                                 >  3:5074/1‑119 (MQ=255)
                                                                                               caacaaCTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTcg                                                   >  7:100121/1‑149 (MQ=255)
                                                                                               caacaaCTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTATTCTACGGCGTCGcctgcctg                                                                   >  7:284473/1‑133 (MQ=255)
                                                                                                 acaacTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTcgcg                                                 <  6:229487/149‑1 (MQ=255)
                                                                                                       ggCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTg                                           >  7:164390/1‑149 (MQ=255)
                                                                                                         cTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTc                                         <  2:134025/149‑1 (MQ=255)
                                                                                                           gTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTGGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGt                                          <  2:173222/146‑1 (MQ=255)
                                                                                                           gTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTGGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGt                                          >  1:173222/1‑146 (MQ=255)
                                                                                                            tCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCgg                                      >  2:350307/1‑149 (MQ=255)
                                                                                                                  ctCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGgcgc                                  <  8:211594/147‑1 (MQ=255)
                                                                                                                    cAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACctgct                                                        <  8:334047/123‑1 (MQ=255)
                                                                                                                    cAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACctgct                                                        >  7:334047/1‑123 (MQ=255)
                                                                                                                     aaGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGtt                             >  3:221542/1‑149 (MQ=255)
                                                                                                                          cacCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCt                        <  4:403885/149‑1 (MQ=255)
                                                                                                                            ccTCGCTGGTGTCGACCATCTCGCTGTTCGAGCTGATGCGCCAGACCAAGCTGGCGTTCAGCGTATCCAAGGAAGCGTTGTTCTTCTAGGGCGGCGCCTGGCTGGGCTACCTGCTGTTCGCGTGAGTGTCGGGGCGCGTGTTCGCCTAc                      <  8:164390/149‑1 (MQ=255)
                                                                                                                            ccTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTa                       >  5:265408/1‑148 (MQ=255)
                                                                                                                              tCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCa                                              >  2:104662/1‑123 (MQ=255)
                                                                                                                              tCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCa                                              <  1:104662/123‑1 (MQ=255)
                                                                                                                                gCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTACAt                    <  1:350307/147‑1 (MQ=255)
                                                                                                                                gCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTACATCg                  <  6:265408/149‑1 (MQ=255)
                                                                                                                                  tGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTACATCGaa                <  8:100121/149‑1 (MQ=255)
                                                                                                                                         gACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTACATCGAACGGCgca         <  6:342808/149‑1 (MQ=255)
                                                                                                                                              tctcGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGgcgc                                  <  5:267702/119‑1 (MQ=255)
                                                                                                                                              tctcGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTACATCGAACGGCGCAGCAAc    >  3:220825/1‑149 (MQ=255)
                                                                                                                                              tctcCCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGgcgc                                  >  6:267702/1‑119 (MQ=255)
                                                                                                                                                tcGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTACATCGAACGGCGCAGCAACCg  >  6:421510/1‑149 (MQ=255)
                                                                                                                                                tcGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTACATCGAACGGCGCAGCAACCg  >  5:367668/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTTCGCCTACATCGAACGGCGCAGCAACCG  >  NC_002947/4084261‑4084553

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: