Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,692,531:1 +A 5.8% coding (628/1623 nt) PP_1488 → methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,692,5311.A 5.8% 282.1 / 8.7 69coding (628/1623 nt)PP_1488methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base . (33/32);  new base A (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

AAGTGATCGAGAATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCTTCAACGAAATGGCCGAGGCACTGG  >  NC_002947/1692394‑1692678
                                                                                                                                          |                                                                                                                                                   
aaGTGATCGAGAATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTAt                                                                                                                                           <  7:304713/148‑1 (MQ=255)
  gTGATCGAGAATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATcac                                                                                                                                        <  3:6921/149‑1 (MQ=255)
  gTGATCGAGAATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATca                                                                                                                                         >  7:100900/1‑148 (MQ=255)
   tGATCGAGAATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCATGGGCGCtgtga                                                                                                                                         >  2:334640/1‑144 (MQ=255)
    gATCGAGAATAATCGGGCCTCTGCCGATGCCGCCACCAACGATGTCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCATGGGCGCtgtga                                                                                                                                         >  2:367471/1‑143 (MQ=255)
    gATCGAGAATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATcaca                                                                                                                                       <  1:35306/148‑1 (MQ=255)
      tCGAGAATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAgcg                                                                                                                                    <  8:243958/149‑1 (MQ=255)
        tagaATAATCGGGGCTCTGCCGAGGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGTGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCTGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGcc                                                                                                                                  <  8:222492/148‑1 (MQ=255)
         agaATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGccc                                                                                                                                 <  6:158586/149‑1 (MQ=255)
          gaATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCa                                                                                                                                >  5:33275/1‑149 (MQ=255)
            ataatCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATg                                                                                                                              <  1:200429/149‑1 (MQ=255)
             taatCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATg                                                                                                                              >  7:63917/1‑148 (MQ=255)
              aatCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATg                                                                                                                              >  3:363184/1‑147 (MQ=255)
               atCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAg                                                                                                                           <  8:100900/149‑1 (MQ=255)
                tCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGc                                                                                                                          >  3:198815/1‑149 (MQ=255)
                tCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGc                                                                                                                          <  5:108040/149‑1 (MQ=255)
                   ggCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTa                                                                                                                       <  3:107344/149‑1 (MQ=255)
                    gCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCATGGGCGCTGTGATAGCGCGCagcagcaaca                                                                                                                      <  1:334640/149‑5 (MQ=38)
                     ccTCTGCCGATGCCGCCACCAACGATGTCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCATGGGCGCTGTGATAGCGCGCagcagcaacag                                                                                                                     <  1:367471/149‑6 (MQ=37)
                      ctctGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCg                                                                                                                    >  6:211203/1‑149 (MQ=255)
                         tGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGt                                                                                                                   <  5:119066/147‑1 (MQ=255)
                            cGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATg                                                                                                              >  7:10935/1‑149 (MQ=255)
                              aTGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGcc                                                                                                            <  2:285637/149‑1 (MQ=255)
                                   gccACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGa                                                                                                       >  4:71088/1‑149 (MQ=255)
                                   gccACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGa                                                                                                       >  1:65647/1‑149 (MQ=255)
                                        caACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGc                                                                                                  >  2:315505/1‑149 (MQ=255)
                                        caACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGc                                                                                                  <  3:71088/149‑1 (MQ=255)
                                         aaCGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCt                                                                                                 >  3:354769/1‑149 (MQ=255)
                                         aaCGATATCGTCAACACAGACAACACGGCAAAGAGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGACGCCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCt                                                                                                 >  8:54604/1‑149 (MQ=255)
                                          aCGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTg                                                                                                >  5:296546/1‑149 (MQ=255)
                                             atatCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCAt                                                                                                                               <  7:105723/115‑1 (MQ=255)
                                             atatCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCAt                                                                                                                               >  8:105723/1‑115 (MQ=255)
                                               atCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCg                                                                                                                         >  4:365753/1‑119 (MQ=255)
                                               atCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCg                                                                                                                         <  3:365753/119‑1 (MQ=255)
                                                  gTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCgg                                                                                        >  5:305206/1‑149 (MQ=255)
                                                  gTCAACGCAGGCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCg                                                                                         <  4:198815/148‑1 (MQ=255)
                                                      aCGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTgcc                                                                                    >  5:271845/1‑147 (MQ=255)
                                                        gCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCt                                                                                  <  1:181539/149‑1 (MQ=255)
                                                        gCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCt                                                                                  <  8:10935/149‑1 (MQ=255)
                                                         cAGTCAACACGGCAAAGTGCAGCATGGTGGTCTCACTGCTGCTTGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTg                                                                                 <  7:264083/149‑1 (MQ=255)
                                                              aaCACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGATCTGAGc                                                                              <  6:33275/147‑1 (MQ=255)
                                                              aaCACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCAt                                                                            >  1:226174/1‑149 (MQ=255)
                                                              aaCACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCAt                                                                            >  1:56830/1‑149 (MQ=255)
                                                               acacGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATg                                                                           >  2:221360/1‑149 (MQ=255)
                                                                cacGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATgc                                                                          >  7:186290/1‑149 (MQ=255)
                                                                      aaaGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCTTGa                                                                    >  1:157180/1‑149 (MQ=255)
                                                                      aaaGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGa                                                                    <  5:211203/149‑1 (MQ=255)
                                                                          ggCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCt                                                                <  4:354769/149‑1 (MQ=255)
                                                                           gcagcaTGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTg                                                               <  6:296546/149‑1 (MQ=255)
                                                                                   gtggtCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGc                                                         >  3:94719/1‑147 (MQ=255)
                                                                                   gtggtCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCaa                                                       <  2:65647/149‑1 (MQ=255)
                                                                                                  ctgGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGt                                          <  8:63917/147‑1 (MQ=255)
                                                                                                   tgGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGc                                        >  6:82869/1‑148 (MQ=255)
                                                                                                       aaTCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGa                                                    >  5:131964/1‑132 (MQ=255)
                                                                                                       aaTCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGa                                                    <  6:131964/132‑1 (MQ=255)
                                                                                                       aaTCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGt                                          >  3:344414/1‑142 (MQ=255)
                                                                                                       aaTCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGt                                          <  4:344414/142‑1 (MQ=255)
                                                                                                       aaTCGTCGCTGCCGCCATCTGCTGCATGTTGCTGATGCGCGCTATCACAGCGCACATGCAGCGTATCGTGCATGACCTCGACAAGCTGCGTTCCGGTGACCTCAGCATCCGCCTGAGCCTGGATCGAAAGGACGAATACGGTGCCATc                                    >  2:352398/1‑148 (MQ=255)
                                                                                                               cTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGc                           >  2:47680/1‑149 (MQ=255)
                                                                                                               cTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGc                           >  8:176484/1‑149 (MQ=255)
                                                                                                                  ccGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCtt                         <  5:362870/148‑1 (MQ=255)
                                                                                                                  ccGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCtt                         >  6:171558/1‑148 (MQ=255)
                                                                                                                    ggCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCTTCaa                      <  2:226174/149‑1 (MQ=255)
                                                                                                                     gCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCTTCAAc                     >  2:105703/1‑149 (MQ=255)
                                                                                                                            gggCCTGTTGCTGCTGCGCGCTTTCACAGCGCCCATGCAGCGTATCGTGCATGCCGTCGAGAAGCTGGGTTCCGGTGACCTGAGCATGCGCGTGAGCCTGGATCGCTAGGACGTATTCGGTGCCATCGAAAGCGGCTTCAACGaaa                 <  7:54604/145‑1 (MQ=255)
                                                                                                                             ggCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCTTCAACGAAATGGc             >  6:200350/1‑149 (MQ=255)
                                                                                                                              gCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCTTCAACGAAATGGcc            <  4:94719/149‑1 (MQ=255)
                                                                                                                                cTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGTCATCGAAAGCGGCTTCAACGAAATGGCCGa          <  1:221360/149‑1 (MQ=255)
                                                                                                                                   ttGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCTTCAACGAAATGGCCGAgg        >  6:160168/1‑148 (MQ=255)
                                                                                                                                       tgctgcGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGa                                                    <  7:1573/100‑1 (MQ=255)
                                                                                                                                       tgctgcGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGa                                                    >  8:1573/1‑100 (MQ=255)
                                                                                                                                        gctgcGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAACTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCTTCAACGAAATGGCCGAGGCACTgg  >  5:33394/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
AAGTGATCGAGAATAATCGGGCCTCTGCCGATGCCGCCACCAACGATATCGTCAACGCAGTCAACACGGCAAAGGGCAGCATGGTGGTCTCACTGCTGCTGGCAATCGTCGCTGCCGGCATCTGCGGCCTGTTGCTGCTGCGCGCTATCACAGCGCCCATGCAGCGTATCGTGCATGCCCTCGACAAGCTGCGTTCCGGTGACCTGAGCATGCGCCTGAGCCTGGATCGCAAGGACGAATTCGGTGCCATCGAAAGCGGCTTCAACGAAATGGCCGAGGCACTGG  >  NC_002947/1692394‑1692678

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