Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,137,497 T→C 5.1% intergenic (+41/‑49) PP_5733 → / → czcC metal‑binding protein/CzcC family metal RND transporter outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,137,4970TC 5.1% 288.5 / 2.9 95intergenic (+41/‑49)PP_5733/czcCmetal‑binding protein/CzcC family metal RND transporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (46/44);  new base C (2/3);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACTCGCTCACTC  >  NC_002947/6137353‑6137639
                                                                                                                                                 |                                                                                                                                              
tCGAGCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAAc                                                                                                                                            <  2:303667/149‑1 (MQ=255)
   aGCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGAt                                                                                                                                         >  1:77308/1‑149 (MQ=255)
    gCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGAt                                                                                                                                         >  8:292930/1‑148 (MQ=255)
     cccGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGa                                                                                                                                       >  2:85368/1‑149 (MQ=255)
      ccGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGAt                                                                                                                                         >  1:381440/1‑145 (MQ=255)
      ccGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGAt                                                                                                                                         <  2:381440/145‑1 (MQ=255)
        gCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGacac                                                                                                                                    <  6:320002/149‑1 (MQ=255)
          aCTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACAcgc                                                                                                                                  >  5:352530/1‑149 (MQ=255)
           cTCACCTTCGCCACTACACCTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACAcgc                                                                                                                                  <  1:85368/149‑1 (MQ=255)
            tCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACgcgc                                                                                                                                <  6:352530/149‑1 (MQ=255)
              aCCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCt                                                                                                                                                <  1:173627/131‑1 (MQ=255)
              aCCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCt                                                                                                                                                >  2:173627/1‑131 (MQ=255)
               ccTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGAGGACACGCGCCtg                                                                                                                             >  6:252887/1‑149 (MQ=255)
                 ttCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCtgtg                                                                                                                           >  2:275613/1‑149 (MQ=255)
                    gCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTgcgc                                                                                                                        >  3:279476/1‑149 (MQ=255)
                    gCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTgcgc                                                                                                                        >  3:284733/1‑149 (MQ=255)
                         tACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGtt                                                                                                                   >  8:184053/1‑149 (MQ=255)
                            acGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGttttt                                                                                                                <  7:292930/149‑1 (MQ=255)
                              gTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTc                                                                                                              <  4:318326/149‑1 (MQ=255)
                              gTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTc                                                                                                              >  4:115341/1‑149 (MQ=255)
                              gTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTc                                                                                                              <  2:84607/149‑1 (MQ=255)
                                ccAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGa                                                                                                            <  8:2073/149‑1 (MQ=255)
                                ccAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCCGAACGATGACACGCGCCTGTGCGCTGGTTTTTTc                                                                                                              >  3:212007/1‑147 (MQ=255)
                                 cAGTCGATACCTACGGTCTGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGa                                                                                                            <  4:320972/148‑1 (MQ=255)
                                 cAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGa                                                                                                            <  5:222277/148‑1 (MQ=255)
                                 cAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAg                                                                                                           <  5:113971/149‑1 (MQ=255)
                                  aGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAg                                                                                                           >  6:316197/1‑148 (MQ=255)
                                   ggcGATACCTACGCGCAGAGCCTTTTCACCAGGGCGCCCCCCGATCCCTGGCGCCCTCCTCCCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCGCGGCGGCCGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAggc                                                                                                         <  4:339710/147‑2 (MQ=255)
                                    tCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTgcgc                                                                                                                        >  2:387025/1‑133 (MQ=255)
                                    tCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTgcgc                                                                                                                        <  1:387025/133‑1 (MQ=255)
                                    tCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGaa                                                                                                        <  2:181619/149‑1 (MQ=255)
                                         aCCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAatcatc                                                                                                   <  3:286561/149‑1 (MQ=255)
                                           cTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTc                                                                                                              <  2:63655/136‑1 (MQ=255)
                                           cTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTc                                                                                                              >  1:63655/1‑136 (MQ=255)
                                           cTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCAtctc                                                                                                 >  6:121694/1‑149 (MQ=255)
                                              cGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCa                                                                                               <  2:77308/148‑1 (MQ=255)
                                                gTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGa                                                                                            <  3:115341/149‑1 (MQ=255)
                                                   agagCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTc                                                                                         <  4:284733/149‑1 (MQ=255)
                                                   agagCCGTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTc                                                                                         <  4:279476/149‑1 (MQ=255)
                                                       ccTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCt                                                                                     <  1:275613/149‑1 (MQ=255)
                                                       ccTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCt                                                                                     <  8:260229/149‑1 (MQ=255)
                                                        cTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCCGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTa                                                                                    >  1:13162/1‑149 (MQ=255)
                                                         ttttAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTAt                                                                                   >  1:133765/1‑149 (MQ=255)
                                                              aCCAGGGCTCCCCCCTATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTTTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTcc                                                                              <  5:1241/149‑1 (MQ=255)
                                                               ccAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTcc                                                                              >  7:117514/1‑148 (MQ=255)
                                                                cAGGGCTCCCACCGATCACTGGCGCCCTACTCGCGACTGATTACTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCACCATGCCTCCCGGATGCACCCa                                                                            >  4:365416/1‑149 (MQ=255)
                                                                 aGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCAt                                                                           <  8:144024/149‑1 (MQ=255)
                                                                      tCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGtttt                                                                                                                 <  8:79632/106‑1 (MQ=255)
                                                                      tCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGtttt                                                                                                                 >  7:79632/1‑106 (MQ=255)
                                                                        cccccGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTAt                                                                                   <  8:42792/134‑1 (MQ=255)
                                                                        cccccGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTAt                                                                                   >  7:42792/1‑134 (MQ=255)
                                                                        cccccGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGc                                                                    >  4:1558/1‑149 (MQ=255)
                                                                         ccccGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAatc                                                                                                      >  4:174755/1‑114 (MQ=255)
                                                                                  cTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGAtgcttgc                                                          <  5:316197/149‑1 (MQ=255)
                                                                                   tGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCAtt                                                                          >  6:372902/1‑132 (MQ=255)
                                                                                     gcgcCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTccc                                                                             <  7:319076/127‑1 (MQ=255)
                                                                                     gcgcCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTccc                                                                             >  8:319076/1‑127 (MQ=255)
                                                                                     gcgcCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTtgctgc                                                       >  8:32132/1‑149 (MQ=255)
                                                                                     gcgcCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTtgctg                                                        <  2:133765/148‑1 (MQ=255)
                                                                                       gcCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGtt                                                                                                                   <  8:336285/87‑1 (MQ=255)
                                                                                       gcCCTCCTCGCGCCAGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGtt                                                                                                                   >  7:336285/1‑87 (MQ=255)
                                                                                          ctcctcGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCt                                                  <  5:28314/149‑1 (MQ=255)
                                                                                           tcctcGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCtt                                                 >  5:104748/1‑149 (MQ=255)
                                                                                            cctcGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTg                                                <  4:102702/149‑1 (MQ=255)
                                                                                             ctcGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTg                                                <  7:184053/148‑1 (MQ=255)
                                                                                              tcGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTgg                                               >  2:352078/1‑148 (MQ=255)
                                                                                                gcgcCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTAt                                                                                   <  8:253878/110‑1 (MQ=255)
                                                                                                gcgcCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTAt                                                                                   >  7:253878/1‑110 (MQ=255)
                                                                                                  gcCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCCGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTc                                          <  2:13162/149‑1 (MQ=255)
                                                                                                    cTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCct                                        >  5:141940/1‑149 (MQ=255)
                                                                                                       aTTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGAtgcttg                                                           <  2:131121/127‑1 (MQ=255)
                                                                                                       aTTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGAtgcttg                                                           >  1:131121/1‑127 (MQ=255)
                                                                                                       aTTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCctg                                       <  3:10042/147‑1 (MQ=255)
                                                                                                        ttCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCctgctg                                    >  6:2842/1‑149 (MQ=255)
                                                                                                        ttCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCctgctg                                    >  5:180770/1‑149 (MQ=255)
                                                                                                        ttCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCCGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCctgct                                     <  4:212007/148‑1 (MQ=255)
                                                                                                          ccTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCtgctgc                                   <  3:4210/148‑1 (MQ=255)
                                                                                                            ttCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCtgt                                                                                                                            <  4:304907/57‑1 (MQ=255)
                                                                                                            ttCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCtgt                                                                                                                            >  3:304907/1‑57 (MQ=255)
                                                                                                            ttCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCgg                                >  4:41383/1‑149 (MQ=255)
                                                                                                            ttCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCgg                                >  5:187806/1‑149 (MQ=255)
                                                                                                            ttCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCgg                                >  6:114906/1‑149 (MQ=255)
                                                                                                               gTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGgcgc                             <  4:192654/149‑1 (MQ=255)
                                                                                                                       tCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGgc                      >  6:48876/1‑149 (MQ=255)
                                                                                                                          cacCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCa                   <  3:1558/149‑1 (MQ=255)
                                                                                                                             cATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCGCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGc                >  5:305854/1‑149 (MQ=255)
                                                                                                                             cATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGgcg                              >  7:74218/1‑135 (MQ=255)
                                                                                                                             cATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGgcg                              <  8:74218/135‑1 (MQ=255)
                                                                                                                             cATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGc                <  6:305854/149‑1 (MQ=255)
                                                                                                                             cATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGc                >  2:393365/1‑149 (MQ=255)
                                                                                                                                  ttCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACt           >  5:245536/1‑149 (MQ=255)
                                                                                                                                       ccTCGGCGGCTGAATGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACTCGCt       >  4:249302/1‑148 (MQ=255)
                                                                                                                                       ccTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACTCg         <  5:2842/146‑1 (MQ=255)
                                                                                                                                       ccTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACTCGCt       >  4:174113/1‑148 (MQ=255)
                                                                                                                                         tCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCg                        >  7:141596/1‑129 (MQ=255)
                                                                                                                                         tCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCg                        <  8:141596/129‑1 (MQ=255)
                                                                                                                                           ggcggcTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACTCGCTCActc  >  2:82898/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TCGAGCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACTCGCTCACTC  >  NC_002947/6137353‑6137639

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: