Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 681,339 C→T 9.1% R646H (CGC→CAC)  PP_0584 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947681,3390CT 9.1% 178.4 / 7.3 66R646H (CGC→CAC) PP_0584methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (35/25);  new base T (3/3);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

TCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGATGTCCTCAGCCACCTGAT  >  NC_002947/681195‑681486
                                                                                                                                                |                                                                                                                                                   
tCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAAc                                                                                                                                                <  2:138665/150‑1 (MQ=255)
tCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAAc                                                                                                                                                >  2:2058011/1‑150 (MQ=255)
   aCCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCtc                                                                                                                                             >  2:206134/1‑149 (MQ=255)
                  gCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATcc                                                                                                                              <  1:1273121/150‑1 (MQ=255)
                           aaGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCAcc                                                                                                                                       >  2:1335496/1‑132 (MQ=255)
                           aaGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCAcc                                                                                                                                       <  1:1335496/132‑1 (MQ=255)
                           aaGCAGCTATTTGGCGCCACACCGGGGCAATTGAAGCTCTGATTTTCGGTTTCCTGGGCCGGCCCTTTCGCGGGGTTACCCGCAAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTg                                                                                                                     <  1:901988/150‑1 (MQ=255)
                            agcagcTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTgg                                                                                                                    >  2:1022890/1‑150 (MQ=255)
                                        gcgcGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCa                                                                                                        >  2:2226494/1‑150 (MQ=255)
                                        gccctacacCGGCGGAATTGAAGCTCTGATTTTCGTTTTCCTGTGCCGCCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGTACAGCCCGGCCa                                                                                                        <  1:947494/145‑1 (MQ=255)
                                                  gggggAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGa                                                                                              <  2:2054575/150‑1 (MQ=255)
                                                      gAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGTGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTcc                                                                                          >  1:2136239/1‑150 (MQ=255)
                                                         ttGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTgg                                                                                       <  1:2058011/150‑1 (MQ=255)
                                                                tctGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGc                                                                                                                                    >  1:2291982/1‑98 (MQ=255)
                                                                tctGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGc                                                                                                                                    <  2:2291982/98‑1 (MQ=255)
                                                                tctGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCa                                                                                                                             <  2:381964/105‑1 (MQ=255)
                                                                tctGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCa                                                                                                                             >  1:381964/1‑105 (MQ=255)
                                                                       ttCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGc                                                                         <  1:2226494/150‑1 (MQ=255)
                                                                                       ggCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCTGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAAc                                                         >  1:1087512/1‑150 (MQ=255)
                                                                                        gCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCtctc                                                                                                 <  2:294014/109‑1 (MQ=255)
                                                                                        gCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCtctc                                                                                                 >  1:294014/1‑109 (MQ=255)
                                                                                           cTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcg                                                     >  1:1905296/1‑150 (MQ=255)
                                                                                           cTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcg                                                     >  2:1628625/1‑150 (MQ=255)
                                                                                           cTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcg                                                     >  1:2057081/1‑150 (MQ=255)
                                                                                            tttCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGTGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcgc                                                    >  2:1621979/1‑150 (MQ=255)
                                                                                            tttCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcgc                                                    >  1:1553170/1‑150 (MQ=255)
                                                                                            tttCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcgc                                                    >  2:1755400/1‑150 (MQ=255)
                                                                                                gcgGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGcca                                                              <  1:218070/136‑1 (MQ=255)
                                                                                                gcgGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGcca                                                              >  2:218070/1‑136 (MQ=255)
                                                                                                gcgGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATcc                                                >  1:2176935/1‑150 (MQ=255)
                                                                                                   gggTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCg                                             >  1:332601/1‑150 (MQ=255)
                                                                                                    ggTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGt                                            <  2:130637/150‑1 (MQ=255)
                                                                                                         tCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATc                                       >  1:1769417/1‑150 (MQ=255)
                                                                                                         tCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATc                                       >  1:1212418/1‑150 (MQ=255)
                                                                                                          ccGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCTCTCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCt                                      <  1:1022890/150‑1 (MQ=255)
                                                                                                            gcgAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTgg                                    <  2:1212418/150‑1 (MQ=255)
                                                                                                              gAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCt                                  >  2:215668/1‑150 (MQ=255)
                                                                                                                 aGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGAt                               <  1:1755400/150‑1 (MQ=255)
                                                                                                                   ggacGGCAAAGGCCATGAATCATCAGGTGCGGAACCTGGGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGAtgt                             >  2:874897/4‑150 (MQ=255)
                                                                                                                    gCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAg                                                                                             <  2:2105474/85‑1 (MQ=255)
                                                                                                                    gCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAg                                                                                             >  1:2105474/1‑85 (MQ=255)
                                                                                                                    gCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATGTTGCTCAtgtg                            >  1:498655/1‑150 (MQ=255)
                                                                                                                     ccGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGc                           >  1:2294384/1‑150 (MQ=255)
                                                                                                                       ggCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCg                         <  1:1628625/150‑1 (MQ=255)
                                                                                                                        gCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTTGCTGATGTGCCgg                        >  1:309169/1‑150 (MQ=255)
                                                                                                                          acaGGCCATGCATCATCAGGTGTGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGtt                      <  2:2136239/150‑1 (MQ=255)
                                                                                                                          acaGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGtt                      <  1:2372592/150‑1 (MQ=255)
                                                                                                                           caGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTg                     >  1:848630/1‑150 (MQ=255)
                                                                                                                            aGGCCATGCCTCTTCCGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGtctg                    <  1:461933/150‑4 (MQ=255)
                                                                                                                            aGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGa                    <  2:2057081/150‑1 (MQ=255)
                                                                                                                             ggCCATGCATCATCAGGTGTGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGAt                   <  1:1621979/150‑1 (MQ=255)
                                                                                                                             ggCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGc                                                                         <  1:829545/96‑1 (MQ=255)
                                                                                                                             ggCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGc                                                                         >  2:829545/1‑96 (MQ=255)
                                                                                                                             ggCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGAt                   >  2:461933/1‑150 (MQ=255)
                                                                                                                             ggCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCATGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGAt                   <  1:166158/150‑1 (MQ=255)
                                                                                                                                     atcatcAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGATGTCCTCAg           >  1:2309063/1‑150 (MQ=255)
                                                                                                                                      tcatcaGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGATGTCCTCAGc          >  1:2197468/1‑150 (MQ=255)
                                                                                                                                        atcaGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGATGTCCTCAGCCa        >  2:2094013/1‑150 (MQ=255)
                                                                                                                                         tcaGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACt                                                        >  2:2060946/1‑101 (MQ=255)
                                                                                                                                         tcaGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACt                                                        <  1:2060946/101‑1 (MQ=255)
                                                                                                                                           aGGTGTGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGATGTc                <  1:1168106/139‑1 (MQ=255)
                                                                                                                                           aGGTGTGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGATGTc                >  2:1168106/1‑139 (MQ=255)
                                                                                                                                             gTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGATGTCCTCAGCCACCTGa   >  2:2369522/1‑150 (MQ=255)
                                                                                                                                             gTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTATCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGt                       >  1:2016327/1‑130 (MQ=255)
                                                                                                                                             gTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTATCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGt                       <  2:2016327/130‑1 (MQ=255)
                                                                                                                                              tGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGAAAGGCTTTGGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGATGTCCTCAGCCACCTGAt  >  1:1241447/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCCGTGAATCTGGCTGATGTGCCGGTTGATGTCCTCAGCCACCTGAT  >  NC_002947/681195‑681486

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