Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 220,561 A→T 5.8% intergenic (+18/‑110) PP_0168 → / → PP_5431 putative surface adhesion protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947220,5610AT 5.8% 266.6 / 6.0 86intergenic (+18/‑110)PP_0168/PP_5431putative surface adhesion protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/43);  new base T (2/3);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

AACGCCGAAGGTGGTATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTGCAACGCAAC  >  NC_002947/220412‑220709
                                                                                                                                                     |                                                                                                                                                    
aaCGCCGAAGGTGGTATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCAGGGCCGca                                                                                                                                                      >  1:1780921/1‑150 (MQ=255)
    ccGAAGGTGGTATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGc                                                                                                                                                  >  2:1889586/1‑150 (MQ=255)
           tggtATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCg                                                                                                                                           >  1:792401/1‑150 (MQ=255)
             gtATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGc                                                                                                                                         >  1:2095149/1‑150 (MQ=255)
              tATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCTGTGCGGCCCCggtg                                                                                                                                        <  2:400610/150‑3 (MQ=255)
                  gCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCggtggtct                                                                                                                                    <  2:2231943/150‑7 (MQ=255)
                   ccAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAAACAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTc                                                                                                                                   >  1:441580/1‑150 (MQ=255)
                      aCGTCGACGTGACCAACAAGCTGGAAGGCGCCACCGTGTTTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGTACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTg                                                                                                                                <  2:802669/150‑1 (MQ=255)
                       cGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACATCGGCCCCGGCTCTTTCCTGc                                                                                                                               <  2:498455/150‑1 (MQ=255)
                          cGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTAc                                                                                                                            <  1:1619533/150‑1 (MQ=255)
                              gTGACCATCAAGCTGGAAGGCGTCAAATGGTCTAACAAGACTATCAACTCGTTGATCAGCGGCGCCGATCCGACGATCATCATCCAAAACAAGGACAGCTGACAGCCAACCGGGGCAGCACAGCGGCCCCGGCTAATTCCGGATACCCgg                                                                                                                        >  1:1120977/1‑149 (MQ=255)
                                   cATCAAGCTGGAAGGCGTCTACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACa                                                                                                                   <  2:1590570/150‑1 (MQ=255)
                                   cATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACa                                                                                                                   >  2:1331915/1‑150 (MQ=255)
                                     tCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAgc                                                                                                                 <  1:1331915/150‑1 (MQ=255)
                                       aaGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAgcgc                                                                                                               <  2:331975/150‑1 (MQ=255)
                                           tGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATAc                                                                                                           <  1:698698/150‑1 (MQ=255)
                                             gAAGGCGACAACTGGTCTAACCAGACGCTCAACTCGTTGATCAGCGGGGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGTACAGCGGCCCCGGCTCTTGCCTGCTACCCGCCAACAGCGCATActc                                                                                                         <  2:1021235/150‑1 (MQ=255)
                                               aGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTg                                                                                                       <  1:1973116/150‑1 (MQ=255)
                                               aGGCGTCAACTGGTCTAAACAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTg                                                                                                       <  2:441580/150‑1 (MQ=255)
                                                ggCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCTGTGCGGCCCCggt                                                                                                                                         >  2:1103525/1‑114 (MQ=255)
                                                ggCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCTGTGCGGCCCCggt                                                                                                                                         <  1:1103525/115‑2 (MQ=255)
                                                   gTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACa                                                                                                                   <  2:989437/134‑1 (MQ=255)
                                                   gTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACa                                                                                                                   >  1:989437/1‑134 (MQ=255)
                                                    tCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCaa                                                                                                                     >  2:1914782/1‑131 (MQ=255)
                                                    tCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCaa                                                                                                                     <  1:1914782/131‑1 (MQ=255)
                                                      aaCTGGTCTAACCAGACGATCAACTCGTTCCTAAGGGGATCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCgg                                                                                                <  2:2000053/150‑1 (MQ=255)
                                                       aCTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCggg                                                                                               >  1:1313267/1‑150 (MQ=255)
                                                        cTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGc                                                                                              >  2:775865/1‑150 (MQ=255)
                                                        cTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGc                                                                                              <  2:2341182/150‑1 (MQ=255)
                                                                aCCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgcgcgc                                                                                      >  2:406518/1‑150 (MQ=255)
                                                                aCCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgcgcgc                                                                                      >  1:938083/1‑150 (MQ=255)
                                                                  cAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGccc                                                                                    >  2:733054/1‑150 (MQ=255)
                                                                  cAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGccc                                                                                    >  1:6952/1‑150 (MQ=255)
                                                                     aCGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGc                                                                                 <  2:1272294/150‑1 (MQ=255)
                                                                          cAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCt                                                                            <  1:1889586/150‑1 (MQ=255)
                                                                            aCTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCTGTGCGGCCCCggtg                                                                                                                                        >  1:2212757/1‑86 (MQ=38)
                                                                            aCTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCTGTGCGGCCCCggtg                                                                                                                                        <  2:2212757/88‑3 (MQ=38)
                                                                              tCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgc                                                                                          <  1:2060702/132‑1 (MQ=255)
                                                                              tCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgc                                                                                          >  2:2060702/1‑132 (MQ=255)
                                                                              tCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                        <  2:1780921/150‑1 (MQ=255)
                                                                              tCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                        <  2:107866/150‑1 (MQ=255)
                                                                               gtcggATCGGCGCCGCTGATCCGACCATCATCATCCCCAACAAGGACAGCTGACAGACCACCGGGGCCGCACATCGGCCCCGGCTCTGTCCTGCCACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCCATCTTTGcc                                                                       >  2:407434/5‑150 (MQ=255)
                                                                                 ttGATCAGCGGCGCCGATCCGACCATCATCATCCCCAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTa                                                                     <  1:524253/150‑1 (MQ=255)
                                                                                 ttGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTa                                                                     >  2:1837925/1‑150 (MQ=255)
                                                                                        gcggcgCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATgctgct                                                              >  1:1215205/1‑150 (MQ=255)
                                                                                        gcggcgCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATgctgct                                                              >  1:2337368/1‑150 (MQ=255)
                                                                                          ggcgCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGt                                                            <  1:775865/150‑1 (MQ=255)
                                                                                             gcCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTa                                                         >  1:1237859/1‑150 (MQ=255)
                                                                                             gcCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTa                                                         >  2:850758/1‑150 (MQ=255)
                                                                                              ccGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTAc                                                        <  1:406518/150‑1 (MQ=255)
                                                                                                    cGACCCTCATCATCCCCAACACGAACAGCGGACACACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCTCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACcaacaa                                                  <  2:463938/150‑1 (MQ=255)
                                                                                                    cGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACcaacaa                                                  <  1:733054/150‑1 (MQ=255)
                                                                                                    cGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACcaacaa                                                  <  1:1837925/150‑1 (MQ=255)
                                                                                                     gACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGgcg                                                                                                    <  2:683954/99‑1 (MQ=255)
                                                                                                     gACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGgcg                                                                                                    >  1:683954/1‑99 (MQ=255)
                                                                                                     gACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACcaacaa                                                  <  1:1634409/149‑1 (MQ=255)
                                                                                                     gACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACcaacaa                                                  >  2:1634409/1‑149 (MQ=255)
                                                                                                     gACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAg                                                 <  2:1215205/150‑1 (MQ=255)
                                                                                                         atcatcatcCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGaaa                                             >  2:1532745/1‑150 (MQ=255)
                                                                                                           catcatcCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGATGTCTAACAACAAAAAaaca                                           >  2:1118749/1‑150 (MQ=255)
                                                                                                            atcatcCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAaacaa                                          <  1:1118749/150‑1 (MQ=255)
                                                                                                             tcatcCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAaacaac                                         >  1:2175848/1‑150 (MQ=255)
                                                                                                               atcCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGt                                       <  2:1237859/150‑1 (MQ=255)
                                                                                                                  cacaACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGt                                       >  2:291924/1‑147 (MQ=255)
                                                                                                                  cacaACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGt                                       <  1:291924/147‑1 (MQ=255)
                                                                                                                     aacaaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgc                                                                                          <  2:507647/93‑1 (MQ=255)
                                                                                                                     aacaaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                        <  2:1440455/111‑1 (MQ=255)
                                                                                                                     aacaaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                        >  1:1440455/1‑111 (MQ=255)
                                                                                                                     aacaaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAACGCATACTCTGGCGCCGGGCCTgc                                                                                          >  1:507647/1‑93 (MQ=255)
                                                                                                                       caaGCACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAggg                               >  2:2203731/1‑150 (MQ=255)
                                                                                                                        aaGGACAGCTGACATACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGa                              <  1:470909/150‑1 (MQ=255)
                                                                                                                        aaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGa                              >  2:2175346/1‑150 (MQ=255)
                                                                                                                          ggACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTAcc                                                       <  1:1116141/123‑1 (MQ=255)
                                                                                                                          ggACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTAcc                                                       >  2:1116141/1‑123 (MQ=255)
                                                                                                                              aGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgcgcgc                                                                                      >  2:1576892/1‑88 (MQ=255)
                                                                                                                              aGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgcgcgc                                                                                      <  1:1576892/88‑1 (MQ=255)
                                                                                                                                 tGACAGACCACCTGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCAt                     <  1:2175346/150‑1 (MQ=255)
                                                                                                                                  gACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCgg                                                                                                <  1:1196184/74‑1 (MQ=255)
                                                                                                                                  gACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCgg                                                                                                >  2:1196184/1‑74 (MQ=255)
                                                                                                                                    cAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                        >  2:1399953/1‑96 (MQ=255)
                                                                                                                                    cAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                        <  1:1399953/96‑1 (MQ=255)
                                                                                                                                      gACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCt                <  2:792401/150‑1 (MQ=255)
                                                                                                                                       accaccGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTa               >  2:1100917/1‑150 (MQ=255)
                                                                                                                                       accaccGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTa               >  2:54331/1‑150 (MQ=255)
                                                                                                                                         caccGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACg             >  1:1036088/1‑150 (MQ=255)
                                                                                                                                          cccGGGGCCACACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCAATGCTGCTCTCTACCAACCCGCAAACAACGTTACGAGGGACATCGCCATGTTCTACGt            <  1:54331/149‑1 (MQ=255)
                                                                                                                                                gCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTgcaacg      <  2:6952/150‑1 (MQ=255)
                                                                                                                                                    cacaGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTgcaacgcaac  <  1:850758/150‑1 (MQ=255)
                                                                                                                                                    cacaGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTgcaacgcaac  <  2:1036088/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACGCCGAAGGTGGTATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTGCAACGCAAC  >  NC_002947/220412‑220709

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