Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 682,317 A→C 7.8% V319G (GTG→GGG)  PP_0584 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947682,3170AC 7.8% 230.1 / 9.7 77V319G (GTG→GGG) PP_0584methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base C (3/3);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

AGGCTGCGAATGGCACCGAGGAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGATGAAGGTCCAGCCCAGC  >  NC_002947/682170‑682460
                                                                                                                                                   |                                                                                                                                               
aGGCGGCGAAGGGCACCGTGAAACTGGTTCAACCAGTTATCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCa                                                                                                                                                 <  1:123522/148‑1 (MQ=255)
  gCTGCGAATGGCACCGAGGAACTGGTTGAACGAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGa                                                                                                                                              <  1:68239/149‑1 (MQ=255)
  gCTGCGAATGGCACCGAGGAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGa                                                                                                                                              >  5:94514/1‑149 (MQ=255)
        aaTGGCACCGAGGAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGTGCACCCGATCaaca                                                                                                                                        <  8:224606/149‑4 (MQ=255)
              aCCGAGGAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCg                                                                                                                                   >  2:213497/1‑148 (MQ=255)
              aCCGAGGAACTGGTTGAACCAGGTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCg                                                                                                                                   >  7:233565/1‑148 (MQ=255)
               ccGAGGAACTTGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGaa                                                                                                                                 <  7:18324/149‑1 (MQ=255)
                cGAGGAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAg                                                                                                                                >  7:88523/1‑149 (MQ=255)
                  aGGAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGc                                                                                                                              <  8:348587/149‑1 (MQ=255)
                   ccaaCTGGTTCAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGcc                                                                                                                             >  8:124076/3‑149 (MQ=255)
                    gAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGCCGCTGTTCATCGGGCGCACGATCAGTTTGGCGCAGGCCg                                                                                                                            >  8:172258/1‑149 (MQ=255)
                    gAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAgg                                                                                                                               >  5:202350/1‑146 (MQ=255)
                    gAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAgg                                                                                                                               <  6:202350/146‑1 (MQ=255)
                               aCCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGc                                                                                                                           >  3:95903/1‑139 (MQ=255)
                               aCCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGc                                                                                                                           <  4:95903/139‑1 (MQ=255)
                                    ttGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGaa                                                                                                            >  4:20048/1‑149 (MQ=255)
                                    ttGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGaa                                                                                                            <  1:213497/149‑1 (MQ=255)
                                       gCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAg                                                                                                         >  3:254087/1‑149 (MQ=255)
                                         cAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAgcg                                                                                                       >  7:90409/1‑149 (MQ=255)
                                         cAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAgcg                                                                                                       >  2:12756/1‑149 (MQ=255)
                                          aGTTGCGCGGTTTCGTCGCGGCAGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCAAGGCAATTGGTGACTCTGTTGATCTGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCTCGAACAgcgc                                                                                                      >  4:49592/1‑149 (MQ=255)
                                           gctGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGcgcg                                                                                                     <  8:233565/147‑1 (MQ=255)
                                                gcgGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCaa                                                                                                >  1:191064/1‑149 (MQ=255)
                                                  ggTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAaca                                                                                              <  1:12756/149‑1 (MQ=255)
                                                      tcgtcgCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAAcacca                                                                                           <  2:191064/148‑1 (MQ=255)
                                                      tcgtcgCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAAcacca                                                                                           <  8:88523/148‑1 (MQ=255)
                                                        gtggctgCCGCTGATTTCCAGGTTTCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTTACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACcaca                                                                                         <  5:308741/142‑1 (MQ=255)
                                                        gtcgCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGTGCGCCCGATCaaca                                                                                                                                        <  1:212838/101‑4 (MQ=255)
                                                        gtcgCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGTGCGCCCGATCaaca                                                                                                                                        >  2:212838/1‑98 (MQ=255)
                                                            cgGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGAt                                                                                    <  5:340400/149‑1 (MQ=255)
                                                              gCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGa                                                                                  <  3:20048/149‑1 (MQ=255)
                                                                cgcgGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGaca                                                                                >  2:100319/1‑149 (MQ=255)
                                                                   ggATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCaa                                                                                                <  8:346675/130‑1 (MQ=255)
                                                                   ggATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCaa                                                                                                >  7:346675/1‑130 (MQ=255)
                                                                      tttCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATc                                                                          <  6:94514/149‑1 (MQ=255)
                                                                       ttCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGAt                                                                           >  6:114138/1‑147 (MQ=255)
                                                                         ccAGGTTGCGGGTCAGGTCACCTTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGc                                                                       >  5:77538/1‑149 (MQ=255)
                                                                         ccAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGc                                                                       >  8:224425/1‑149 (MQ=255)
                                                                          cAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGAGACCGATCAGc                                                                       <  8:90409/148‑1 (MQ=255)
                                                                           aGGTTGCGGGTCAGGTCACCCTCACCTTTGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGTCGAAGGCCGCACCGACCACCGCTAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCa                                                                     <  3:49592/149‑1 (MQ=255)
                                                                           aGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAgcgc                                                                                                      >  8:310919/1‑116 (MQ=255)
                                                                           aGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAgcgc                                                                                                      <  7:310919/116‑1 (MQ=255)
                                                                            ggTTGCGGGTCAGGTCACCATCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAg                                                                    <  6:77538/149‑1 (MQ=255)
                                                                              ttGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGt                                                                  >  2:347805/1‑149 (MQ=255)
                                                                                      caggtcaCCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCggg                                                          >  5:301854/1‑149 (MQ=255)
                                                                                       aggtcaCCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGc                                                                                                   <  5:74421/107‑1 (MQ=255)
                                                                                       aggtcaCCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGc                                                                                                   >  6:74421/1‑107 (MQ=255)
                                                                                       aggtcaCCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCggg                                                          <  7:124076/148‑1 (MQ=255)
                                                                                              ccTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTc                                                  >  4:332483/1‑149 (MQ=255)
                                                                                                tCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCACCACCACAGCGATGACACCGATCAGCCAGGTGAGCCTGGTGGTGGTCTg                                                >  8:314893/1‑149 (MQ=255)
                                                                                                  aCCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGc                                               <  5:114138/148‑1 (MQ=255)
                                                                                                        ggCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACtt                                        <  1:100319/149‑1 (MQ=255)
                                                                                                         gCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTc                                       >  5:69601/1‑149 (MQ=255)
                                                                                                         gCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTc                                       >  7:257821/1‑149 (MQ=255)
                                                                                                          cAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGCCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCg                                      <  7:224425/149‑1 (MQ=255)
                                                                                                             tGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCtt                                   <  6:301854/149‑1 (MQ=255)
                                                                                                               tCTTCGAGGCCATTGGTGACGCTGTTGATCGGGCGAACGATCAGTTTGGCGAAGGCAGCACCGACCACAGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTg                                 <  2:111197/149‑1 (MQ=255)
                                                                                                               tCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTg                                 <  7:172258/149‑1 (MQ=255)
                                                                                                                         cATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGa                            <  2:259938/144‑1 (MQ=255)
                                                                                                                         cATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGa                            >  1:259938/1‑144 (MQ=255)
                                                                                                                          aTTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGc                      <  4:254087/149‑1 (MQ=255)
                                                                                                                              gTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACCCCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCAGTTCGCTGTGTTCGATCAGGCCGAt                  >  4:40414/1‑149 (MQ=255)
                                                                                                                              gTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGAt                           >  3:320891/1‑140 (MQ=255)
                                                                                                                              gTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGAt                           <  4:320891/140‑1 (MQ=255)
                                                                                                                              gTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGAt                  >  2:160678/1‑149 (MQ=255)
                                                                                                                              gTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGAt                  >  6:253220/1‑149 (MQ=255)
                                                                                                                                  ccaaactGATCGTGCGCCCGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGATg                 >  1:133300/7‑146 (MQ=255)
                                                                                                                                   caaactGATCGTGCGCCCGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGGGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGATGAAgg             >  5:334669/6‑149 (MQ=255)
                                                                                                                                         gATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCa                                              >  3:13080/1‑110 (MQ=255)
                                                                                                                                         gATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCa                                              <  4:13080/110‑1 (MQ=255)
                                                                                                                                         gATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGtt                               >  5:276521/1‑125 (MQ=255)
                                                                                                                                         gATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGtt                               <  6:276521/125‑1 (MQ=255)
                                                                                                                                         gATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGATGAAGGTccagc       >  3:86409/1‑149 (MQ=255)
                                                                                                                                          aTCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGATGAAGGTccagcc      >  4:227936/1‑149 (MQ=255)
                                                                                                                                          aTCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGATGAAGGTccagc       <  8:59325/148‑1 (MQ=255)
                                                                                                                                           tCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGtt                               <  8:331221/123‑1 (MQ=255)
                                                                                                                                           tCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGtt                               >  7:331221/1‑123 (MQ=255)
                                                                                                                                              tgcgcCCGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGATGAAGGTccagcccagc  <  2:133300/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGGCTGCGAATGGCACCGAGGAACTGGTTGAACCAGTTGGCCAGTTGCGCGGTTTCGTCGCGGCCGCGGATTTCCAGGTTGCGGGTCAGGTCACCCTCACCTTGGGCAATGTCTTCCAGGCCATTGGTGACGCTGTTGATCGGGCGCACGATCAGTTTGGCGAAGGCCGCACCGACCACCGCGAACAGCGCGGCCAACACCACAGCGATGACACCGATCAGCCAGGTGAGCCGGGTGGTGGTCTGCATCACTTCGCTTTGTTCGATCAGGCCGATGAAGGTCCAGCCCAGC  >  NC_002947/682170‑682460

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