Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,697,314 C→A 6.9% 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 6.9% 271.8 / 7.3 86intergenic (+24/‑53)hslV/hslUATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base A (2/4);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CCCTGCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAAC  >  NC_002947/5697166‑5697462
                                                                                                                                                    |                                                                                                                                                    
cccTGCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTcc                                                                                                                                                      >  2:315191/1‑149 (MQ=255)
  ctgctgAACAAGACCGACCTCTCGGCCTGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCg                                                                                                                                                    >  1:11542/1‑149 (MQ=255)
     ctgAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCgg                                                                                                                                                   <  1:170321/147‑1 (MQ=255)
        aaCAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCAGTGaa                                                                                                                                               <  2:8398/148‑1 (MQ=255)
         aCAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                             <  6:61402/149‑1 (MQ=255)
          cAAGACCTACCTCTCGGCCAGCGAAATCGCCGAGACCTCCCTGAACAACGCCGGCGACATCATCCTCATCACCAACCACAACCTGACCATCGAGGAGCACTACCAGGCCGACTGAATCAGCTGTTCTGGCGCCCATTCCAGGTGAAggg                                                                                                                                            >  2:25168/1‑149 (MQ=255)
          cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAggg                                                                                                                                            >  6:375986/1‑149 (MQ=255)
          cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCCGCGACATCTGCGTCTTCACCACCCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCACCTGGACTGGCGCCCCTTCCCGGTGAAggg                                                                                                                                            >  7:107933/1‑149 (MQ=255)
           aaGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgggg                                                                                                                                           <  5:292663/149‑1 (MQ=255)
           aaGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAgggg                                                                                                                                           <  3:368608/149‑1 (MQ=255)
            aGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAGGGGc                                                                                                                                          <  5:375986/149‑1 (MQ=255)
              accgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTTATTCCGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCtt                                                                                                                                        <  5:29681/149‑1 (MQ=255)
               ccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCAGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAGGGGcgcc                                                                                                                                       <  8:349131/149‑4 (MQ=255)
                 gaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGTGCTTTTc                                                                                                                                     <  3:103039/149‑1 (MQ=255)
                 gaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTc                                                                                                                                     <  3:54257/149‑1 (MQ=255)
                 gaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAGGGGcgccag                                                                                                                                     <  3:320388/149‑6 (MQ=255)
                  accTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTcc                                                                                                                                    <  7:201793/149‑1 (MQ=255)
                      ctcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTGACCAACCACAAGCTGACGATCGAGGAGAAGGACCTGGCCGACTGATTGAGCTGTTTTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCGACGc                                                                                                                                >  3:284979/1‑149 (MQ=255)
                      ctcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGCCCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTGCCGGTGAAGGGGCTTTTCCACGc                                                                                                                                >  3:289986/1‑149 (MQ=255)
                         ggCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCAGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa                                                                                                                             >  3:259892/1‑149 (MQ=255)
                             cgcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCa                                                                                                                         >  2:243810/1‑149 (MQ=255)
                             cgcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCa                                                                                                                         >  1:336897/1‑149 (MQ=255)
                              gcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCAc                                                                                                                        <  8:93061/149‑1 (MQ=255)
                              gcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCAc                                                                                                                        >  7:323239/1‑149 (MQ=255)
                                gAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCt                                                                                                                      <  5:87179/149‑1 (MQ=255)
                                 aaaTCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTg                                                                                                                     <  8:299526/149‑1 (MQ=255)
                                  aaTCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGc                                                                                                                    <  6:269828/149‑1 (MQ=255)
                                      gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGCCCATCGAGGAGCAGGACCTGGCCGACTGATACCGCTGTTCTGGCGCCCCTGCCCGGTGAAGGGGCTTTTCCAGGCTGACTCACCTGccc                                                                                                                  >  4:44761/1‑147 (MQ=255)
                                         gagaCCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGa                                                                                                             <  8:323239/149‑1 (MQ=255)
                                              cGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTa                                                                                                        >  3:193453/1‑149 (MQ=255)
                                               gCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTAt                                                                                                       >  5:341535/1‑149 (MQ=255)
                                                cccTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTAtt                                                                                                      >  5:96016/1‑149 (MQ=255)
                                                 ccTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTAtt                                                                                                      >  6:188154/1‑148 (MQ=255)
                                                 ccTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTAtt                                                                                                      >  1:139600/1‑148 (MQ=255)
                                                    gAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCGTCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCt                                                                                                                      <  6:275622/129‑1 (MQ=255)
                                                    gAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCGTCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCt                                                                                                                      >  5:275622/1‑129 (MQ=255)
                                                    gAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATc                                                                                                  <  4:243543/149‑1 (MQ=255)
                                                      aCATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGt                                                                                                         <  6:164112/140‑1 (MQ=255)
                                                      aCATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGt                                                                                                         >  5:164112/1‑140 (MQ=255)
                                                      aCATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCAt                                                                                                <  6:96016/149‑1 (MQ=255)
                                                       cATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATg                                                                                               >  3:351864/1‑149 (MQ=255)
                                                        aTCGCCGGCGATATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGt                                                                                              <  1:315191/149‑1 (MQ=255)
                                                          cGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTcc                                                                                            <  4:193453/149‑1 (MQ=255)
                                                          cGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTc                                                                                             >  4:345650/1‑148 (MQ=255)
                                                            ccGTCGACATATGCGTGTTCACCAACCACAACCTCAACATCGAGGATCAGGACCTGGCCGACTGATTAAGAAGTACTGGCGCCCCTTCACGCTGAATGGTATTATCCACGCTCCCTCACCTCCACAAGGACCGTATTGATAAAGTcc                                                                                            >  2:130015/1‑147 (MQ=255)
                                                                  aCATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGTACCTGGCCGACTGATTCAGCTCTTCTGGCGCCCCTTCCCGGTGCAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGAcccc                                                                                    <  8:317195/149‑1 (MQ=255)
                                                                        gCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGGGATGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                                >  8:318837/1‑147 (MQ=255)
                                                                         cGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                                <  2:139600/146‑1 (MQ=255)
                                                                           tgtTGACCAACCACAACCTTACCATCGAGGAGCTGGACCTGGCCGTCTGCTTCAGCTGTTCTGGCGCCGCTTGCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTTTTGATCATGTCCATGACCCCCCGCGAGAt                                                                           <  8:107933/149‑1 (MQ=255)
                                                                            gtTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAACGACCGTATTGATCATGTCCATGACCCCCCGCGAGATc                                                                          <  2:11542/149‑1 (MQ=255)
                                                                             ttCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATc                                                                          <  4:351864/148‑1 (MQ=255)
                                                                               cACCAACCACTACCTGACCATCGAGGAGCAGGACCTTGCCGACTGATTCAGCTGTTCTGGCGCCCCTTGCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGt                                                                        <  1:25168/148‑1 (MQ=255)
                                                                                accaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc                                                                       >  6:191901/1‑148 (MQ=255)
                                                                                accaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc                                                                       >  2:144515/1‑148 (MQ=255)
                                                                                 ccaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCACCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCa                                                                     >  7:346292/1‑149 (MQ=255)
                                                                                  caaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCAc                                                                    >  2:86217/1‑149 (MQ=255)
                                                                                     ccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGaa                                                                 <  5:188154/149‑1 (MQ=255)
                                                                                            cTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCCCGAGATCGTCCACGAACTCAAcc                                                          >  3:39727/1‑149 (MQ=255)
                                                                                                 cATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCca                                                      >  4:257811/1‑148 (MQ=255)
                                                                                                  aTCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCcaca                                                    >  2:300139/1‑149 (MQ=255)
                                                                                                     gaggagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatc                                                  >  6:353766/1‑148 (MQ=255)
                                                                                                      aggagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCCcatcat                                                <  7:318837/149‑1 (MQ=255)
                                                                                                        gagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCg                                              <  1:300139/149‑1 (MQ=255)
                                                                                                             ggACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGTCTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAg                                         <  2:116908/149‑1 (MQ=255)
                                                                                                              gACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATAGGCCAg                                         >  1:127354/1‑148 (MQ=255)
                                                                                                                    gCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGgacgac                                   <  6:341535/148‑1 (MQ=255)
                                                                                                                    gCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGacgacg                                  >  5:364220/1‑149 (MQ=255)
                                                                                                                        aCTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCaa                              <  6:364220/149‑1 (MQ=255)
                                                                                                                            aTTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcgc                          >  5:177132/1‑149 (MQ=255)
                                                                                                                                aGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGc                                 <  7:340331/138‑1 (MQ=255)
                                                                                                                                aGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCg                      >  3:111231/1‑149 (MQ=255)
                                                                                                                                aGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCAACGAACTCAACCGCCACATCATCGGCCAGGACGACGc                                 >  8:340331/1‑138 (MQ=255)
                                                                                                                                 gCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCaa                              <  7:348067/140‑1 (MQ=255)
                                                                                                                                 gCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCaa                              >  8:348067/1‑140 (MQ=255)
                                                                                                                                   tGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                                <  4:226847/88‑1 (MQ=255)
                                                                                                                                   tGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                                >  3:226847/1‑88 (MQ=255)
                                                                                                                                    gTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGc                  >  4:359636/1‑149 (MQ=255)
                                                                                                                                     ttCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCCCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGc                  >  8:224324/1‑148 (MQ=255)
                                                                                                                                       cTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCa                <  1:173664/148‑1 (MQ=255)
                                                                                                                                            gcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGccc          >  5:24371/1‑149 (MQ=255)
                                                                                                                                             ccccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCt         <  7:224324/149‑1 (MQ=255)
                                                                                                                                               ccatCCCGGTGACGGGGGTTTTCCACGCTGACTCACCTGCCCAAGGTCCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATAGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgc       <  4:284979/146‑1 (MQ=255)
                                                                                                                                               ccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCGACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgc       <  1:86217/149‑1 (MQ=255)
                                                                                                                                               ccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgc       <  5:353766/149‑1 (MQ=255)
                                                                                                                                               ccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgc       >  5:269975/1‑149 (MQ=255)
                                                                                                                                               ccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgc       >  3:243315/1‑149 (MQ=255)
                                                                                                                                               ccTGCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCTATATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgc       <  4:289986/149‑1 (MQ=255)
                                                                                                                                                    ccGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAAc  <  8:346292/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCTGCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAAC  >  NC_002947/5697166‑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: ‑