Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,697,316 G→C 7.4% intergenic (+26/‑51) 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,3160GC 7.4% 255.4 / 8.7 81intergenic (+26/‑51)hslV/hslUATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base G (36/39);  new base C (3/3);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

CTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCG  >  NC_002947/5697171‑5697464
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                                                              <  8:213179/149‑1 (MQ=255)
     cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAggg                                                                                                                                              >  7:185243/1‑149 (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                                                                                                                         >  7:22601/1‑140 (MQ=255)
                                   cGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCAcc                                                                                                                         <  8:22601/140‑1 (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                                                                  >  3:46657/1‑148 (MQ=255)
                                                                                  acaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAAc                                                                  <  4:46657/148‑1 (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                                                  >  6:74962/1‑149 (MQ=255)
                                                                                                 aggagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat                                                  >  2:72651/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             <  7:130539/149‑1 (MQ=255)
                                                                                                                                      cgcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGcc             <  6:348268/149‑1 (MQ=255)
                                                                                                                                       gcCCCTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGccc            >  7:355928/1‑149 (MQ=255)
                                                                                                                                       gcCCCTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGccc            >  2:59593/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)
                                                                                                                                                 ggTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc                                                                         <  5:241331/78‑1 (MQ=255)
                                                                                                                                                 ggTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc                                                                         >  6:241331/1‑78 (MQ=255)
                                                                                                                                                 ggTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCg  <  3:209371/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCG  >  NC_002947/5697171‑5697464

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: ‑