Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,697,314 C→A 7.6% intergenic (+24/‑53) hslV → / → hslU ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,697,3140CA 7.6% 236.5 / 6.7 79intergenic (+24/‑53)hslV/hslUATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base A (3/3);  total (38/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 = 8.75e-01

CTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAAC  >  NC_002947/5697171‑5697462
                                                                                                                                               |                                                                                                                                                    
ctgAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGt                                                                                                                                                  <  1:374723/148‑1 (MQ=255)
   aaCAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAg                                                                                                                                              <  3:261482/149‑1 (MQ=255)
   aaCAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAg                                                                                                                                              >  3:102278/1‑149 (MQ=255)
   aaCAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAg                                                                                                                                              <  2:76643/149‑1 (MQ=255)
    aCAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAgg                                                                                                                                             >  2:138722/1‑149 (MQ=255)
     cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAggg                                                                                                                                            >  7:185243/1‑149 (MQ=255)
     cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAggg                                                                                                                                            <  8:213179/149‑1 (MQ=255)
      aaGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgggg                                                                                                                                           <  1:174539/149‑1 (MQ=255)
      aaGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgggg                                                                                                                                           >  1:107909/1‑149 (MQ=255)
        gaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAGGGGcg                                                                                                                                         <  1:138722/149‑2 (MQ=255)
         accgaccTCTCTGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAGGGGcgc                                                                                                                                        <  2:238877/149‑3 (MQ=255)
         accgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAGGGGcgc                                                                                                                                        <  5:292684/149‑3 (MQ=255)
                 ctcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGc                                                                                                                                          >  4:176774/1‑139 (MQ=255)
                 ctcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGc                                                                                                                                          <  3:176774/139‑1 (MQ=255)
                  tcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCACCTTCCCGGTGAAGGGGCTTTTCCACGCt                                                                                                                               >  6:313545/1‑149 (MQ=255)
                    ggCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCt                                                                                                                               <  5:335364/147‑1 (MQ=255)
                    ggCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa                                                                                                                             <  8:185243/149‑1 (MQ=255)
                       ccGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                             >  5:260373/1‑130 (MQ=255)
                       ccGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                             <  6:260373/130‑1 (MQ=255)
                        cgcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCa                                                                                                                         >  2:205692/1‑149 (MQ=255)
                        cgcgAAATCGCCGAGACCGCCCTGAACATCGACGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTTACTCa                                                                                                                         >  5:335724/1‑149 (MQ=255)
                         gcgAACTCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCAc                                                                                                                        <  8:237815/149‑1 (MQ=255)
                               tCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGccc                                                                                                                  <  2:107909/149‑1 (MQ=255)
                                cGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGccc                                                                                                                  >  4:172671/1‑148 (MQ=255)
                                   cGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCAcc                                                                                                                       <  8:22601/140‑1 (MQ=255)
                                   cGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCAcc                                                                                                                       >  7:22601/1‑140 (MQ=255)
                                    gagaCCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCttt                                                                                                                                       >  3:284280/1‑123 (MQ=255)
                                    gagaCCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGCCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCttt                                                                                                                                       <  4:284280/123‑1 (MQ=255)
                                      gaCCGCCCTGAACCTCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTACCACGCTGACTCACCTGCCCAAGGAcc                                                                                                           <  8:65624/149‑1 (MQ=255)
                                      gaCCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGAcc                                                                                                           <  3:160979/149‑1 (MQ=255)
                                            ccTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCATGCt                                                                                                                               <  7:345112/123‑1 (MQ=255)
                                            ccTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGATGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCAGGCt                                                                                                                               >  8:345112/1‑123 (MQ=255)
                                                aaCATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATc                                                                                                  <  5:313545/148‑1 (MQ=255)
                                                    tCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGAt                                                                                                   >  3:229970/1‑143 (MQ=255)
                                                    tCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGAt                                                                                                   <  4:229970/143‑1 (MQ=255)
                                                    tCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTc                                                                                             >  6:194977/1‑149 (MQ=255)
                                                      gCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTcc                                                                                            <  6:45135/148‑1 (MQ=255)
                                                      gCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCa                                                                                           <  1:43113/149‑1 (MQ=255)
                                                       ccGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAt                                                                                          <  3:172671/149‑1 (MQ=255)
                                                        cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAt                                                                                          >  5:246813/1‑148 (MQ=255)
                                                                    cGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGAcccccccc                                                                                >  3:121596/1‑144 (MQ=255)
                                                                    cGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCAGGCCGAATGATTCAGCTGTTATGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                                >  7:9113/1‑146 (MQ=255)
                                                                       gtTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTc                                                                                                                          <  4:377482/101‑1 (MQ=255)
                                                                       gtTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTc                                                                                                                          >  3:377482/1‑101 (MQ=255)
                                                                         tCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGt                                                                        >  3:257691/1‑149 (MQ=255)
                                                                           accaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCAAATGACCGTATTGATCATGTCGATGACCCCCCGCGAGATCGTc                                                                       <  6:335724/148‑1 (MQ=255)
                                                                             caaccaCAACCTGACCATCGTGGAGCAGTCCCTGGCCGACTGATTCAGCTGCTCTGGCGCCCCTTCCCGGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCAc                                                                    <  8:339561/149‑1 (MQ=255)
                                                                              aaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATAATGTCCATGACCCCCCGCGAGATCGTCCAc                                                                    >  3:13845/1‑148 (MQ=255)
                                                                                 cacaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCCCGAGATCGTCCACGAAc                                                                <  4:121596/149‑1 (MQ=255)
                                                                                  acaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAAc                                                                <  4:46657/148‑1 (MQ=255)
                                                                                  acaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAAc                                                                >  3:46657/1‑148 (MQ=255)
                                                                                  acaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACt                                                               <  1:205692/149‑1 (MQ=255)
                                                                                     aCCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCaa                                                            <  5:194977/149‑1 (MQ=255)
                                                                                          aCCATCGAGGAGCAGGACCTGGCCGACTTATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATTTCCATGACCCCCCGCGTTATCGTCCACGAACTCAACCGc                                                        <  8:235924/148‑1 (MQ=255)
                                                                                            cATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCca                                                      >  8:130539/1‑148 (MQ=255)
                                                                                            cATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCca                                                      >  4:209371/1‑148 (MQ=255)
                                                                                             aTCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCcaca                                                    >  2:350149/1‑149 (MQ=255)
                                                                                              tCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcat                                                   >  5:348268/1‑149 (MQ=255)
                                                                                                 aggagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat                                                >  2:72651/1‑149 (MQ=255)
                                                                                                 aggagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat                                                >  6:74962/1‑149 (MQ=255)
                                                                                                            cTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGgac                                      >  2:259953/1‑148 (MQ=255)
                                                                                                                ccGACTGATTCAGCTGTTGTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTTTTGATCATGTCCATGACCCCCCGCGAGATCGTGCACGAACTCAACCGCCACATCATCTGCCAGGACGACGc                                 <  8:9113/149‑1 (MQ=255)
                                                                                                                   aCTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCCACCGCCACATCATCGCCCACGCCGACGCCac                              >  4:211299/1‑148 (MQ=255)
                                                                                                                         tCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcgcgc                        >  8:153835/1‑149 (MQ=255)
                                                                                                                           aGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTGTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCg                      <  8:32531/149‑1 (MQ=255)
                                                                                                                            gCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGt                     >  4:196817/1‑149 (MQ=255)
                                                                                                                             cTGTTCTGGCGCCCCTTGCCGGGTAAGGGGCTTTTCCACGCTGTCTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTg                    <  3:211299/149‑1 (MQ=255)
                                                                                                                             cTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAgg                                        <  7:285742/129‑1 (MQ=255)
                                                                                                                             cTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAgg                                        >  8:285742/1‑129 (MQ=255)
                                                                                                                             cTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTg                    >  5:228060/1‑149 (MQ=255)
                                                                                                                              tGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACatcatc                                               >  1:5500/1‑121 (MQ=255)
                                                                                                                              tGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACatcatc                                               <  2:5500/121‑1 (MQ=255)
                                                                                                                                     gcgcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCAt               <  1:259953/146‑1 (MQ=255)
                                                                                                                                      cgcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGcc           <  6:348268/149‑1 (MQ=255)
                                                                                                                                      cgcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGcc           <  7:130539/149‑1 (MQ=255)
                                                                                                                                       gcCCCTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGccc          >  2:59593/1‑149 (MQ=255)
                                                                                                                                       gcCCCTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGccc          >  7:355928/1‑149 (MQ=255)
                                                                                                                                          ccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgc       <  1:350149/149‑1 (MQ=255)
                                                                                                                                           ctacccGGTGAAGGGGCTTTTCCACGCTGACTCACCTCCCCAAGGACCGAATTCATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgcg      >  3:173059/4‑149 (MQ=255)
                                                                                                                                            ttCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCCAGGACCGTATTGATCATGTCCATGCCCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgcg      >  4:275626/1‑148 (MQ=255)
                                                                                                                                               ccGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAAc  <  4:257691/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAAC  >  NC_002947/5697171‑5697462

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