Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,685,967 G→C 6.9% R328G (CGC→GGC) ‡ PP_4988 ← signal transduction histidine kinase CheA/CheY‑like receiver

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,685,9670GC 6.9% 222.5 / 7.6 72R328G (CGC→GGC) ‡PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base C (3/2);  total (37/35)
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.66e-01

TCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTCCAAAGGGCCGA  >  NC_002947/5685819‑5686116
                                                                                                                                                    |                                                                                                                                                     
tcatcaATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGt                                                                                                                                                      <  1:362629/150‑1 (MQ=255)
  atcaATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTg                                                                                                                                                    >  2:294275/1‑150 (MQ=255)
     aaTCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCg                                                                                                                                                 <  2:729956/150‑1 (MQ=255)
     aaTCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCg                                                                                                                                                 >  2:1230668/1‑150 (MQ=255)
     aaTCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCg                                                                                                                                                 <  1:837742/150‑1 (MQ=255)
         ggcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGc                                                                                                                                             <  2:1438806/150‑1 (MQ=255)
          gcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGcc                                                                                                                                            >  2:65656/1‑150 (MQ=255)
              ctctTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCGGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTc                                                                                                                                        <  2:1139459/150‑1 (MQ=255)
               tctTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCg                                                                                                                                       <  1:531323/150‑1 (MQ=255)
               tctTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCg                                                                                                                                       >  2:1175173/1‑150 (MQ=255)
                 ttCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGa                                                                                                                                     <  1:546357/150‑1 (MQ=255)
                  tcctccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGAc                                                                                                                                    <  1:65656/150‑1 (MQ=255)
                      ccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCa                                                                                                                                <  2:163507/150‑1 (MQ=255)
                       cAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                               <  1:1424231/150‑1 (MQ=255)
                       cAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                               <  2:1229989/150‑1 (MQ=255)
                       cAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                               <  2:556754/150‑1 (MQ=255)
                       cAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                               >  2:1028434/1‑150 (MQ=255)
                         ggTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGt                                                                                                                             >  2:1238256/1‑150 (MQ=255)
                              tgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCg                                                                                                                        >  2:74744/1‑150 (MQ=255)
                                    caacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATc                                                                                                                  >  2:752372/1‑150 (MQ=255)
                                          tGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATcc                                                                                                            >  1:1418560/1‑150 (MQ=255)
                                             acaTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCATGTAATCGCCGATCAGATCCTTg                                                                                                         >  1:1016489/1‑150 (MQ=255)
                                             acaTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTg                                                                                                         <  2:1552590/150‑1 (MQ=255)
                                               aTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCTGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTCGc                                                                                                       >  2:1435394/1‑150 (MQ=255)
                                               aTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGc                                                                                                       >  2:1386625/1‑150 (MQ=255)
                                               aTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGc                                                                                                       >  1:647638/1‑150 (MQ=255)
                                               aTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGc                                                                                                       <  1:1230668/150‑1 (MQ=255)
                                                  gCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgg                                                                                                                                                <  2:1506925/106‑1 (MQ=255)
                                                  gCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgg                                                                                                                                                >  1:1506925/1‑106 (MQ=255)
                                                  gCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATcc                                                                                                            >  2:1395498/1‑142 (MQ=255)
                                                  gCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATcc                                                                                                            <  1:1395498/142‑1 (MQ=255)
                                                   cTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCctgc                                                                                                   >  1:712995/1‑150 (MQ=255)
                                                         cAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAac                                                                                             <  2:712995/150‑1 (MQ=255)
                                                              gACAAACAGCTCGGCCACGGCCGGCAGCATGAACATGGCTAATGCCCCGCGGATCTGCTGCAACGGCCCGGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGGAACACGtt                                                                                        >  2:45561/1‑150 (MQ=255)
                                                              gACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGtt                                                                                        >  2:517255/1‑150 (MQ=255)
                                                                cAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTc                                                                                                                                        <  2:1538776/100‑1 (MQ=255)
                                                                cAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTc                                                                                                                                        >  1:1538776/1‑100 (MQ=255)
                                                                 aaaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAg                                                                                     <  1:157333/150‑1 (MQ=255)
                                                                        tCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                                    >  1:1315322/1‑144 (MQ=255)
                                                                        tCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                                    <  2:1315322/144‑1 (MQ=255)
                                                                          ggCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATc                                                                                                                  <  1:1068682/112‑1 (MQ=255)
                                                                          ggCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATc                                                                                                                  >  2:1068682/1‑112 (MQ=255)
                                                                          ggCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGattcatt                                                                            <  1:1238256/150‑1 (MQ=255)
                                                                            ccACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCCAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGa                                                                                  <  2:55544/142‑1 (MQ=255)
                                                                            ccACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGa                                                                                  >  1:55544/1‑142 (MQ=255)
                                                                             cACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCctgct                                                                                                  <  1:436503/125‑1 (MQ=255)
                                                                             cACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCctgct                                                                                                  >  2:436503/1‑125 (MQ=255)
                                                                              gtggccggcAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGc                                                                        <  2:289063/148‑1 (MQ=255)
                                                                              aCGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCCACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGc                                                                        <  2:1418560/150‑1 (MQ=255)
                                                                                ggccggcAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCAc                                                                      <  1:1173399/150‑1 (MQ=255)
                                                                                    ggcAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGc                                                                  >  2:1012930/1‑150 (MQ=255)
                                                                                      cagcaTCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGATGCCCGCCCCAGGCTGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTg                                                                <  2:951396/150‑1 (MQ=255)
                                                                                         caTCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgg                                                                                                                                                <  2:649907/67‑1 (MQ=255)
                                                                                         caTCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgg                                                                                                                                                >  1:649907/1‑67 (MQ=255)
                                                                                           tCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGGCCCCTGGTAATCGCCGATAAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTa                                                           >  1:108782/1‑150 (MQ=255)
                                                                                                  ggCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGc                                                    <  1:710248/150‑1 (MQ=255)
                                                                                                   gCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCg                                                                                                                                                 <  1:367763/56‑1 (MQ=18)
                                                                                                   gCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCg                                                                                                                                                 >  2:367763/1‑56 (MQ=18)
                                                                                                            ccGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAg                                          >  1:1521883/1‑150 (MQ=255)
                                                                                                                   tCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCa                                   <  1:1028434/150‑1 (MQ=255)
                                                                                                                     tgctgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAgg                                 >  2:1435840/1‑150 (MQ=255)
                                                                                                                          cAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGcc                            >  1:309443/1‑150 (MQ=255)
                                                                                                                          cAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGATGcc                            >  2:1588758/1‑150 (MQ=255)
                                                                                                                           aaCAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCt                           >  2:1640617/1‑150 (MQ=255)
                                                                                                                            aCAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCtt                          >  2:1466886/1‑150 (MQ=255)
                                                                                                                              aGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTgg                        <  1:752372/150‑1 (MQ=255)
                                                                                                                               gCCCCGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGc                       >  2:1226127/1‑150 (MQ=255)
                                                                                                                                 cccGGCAAGGGCTGCGGGCGTTGCCGGGGCCAGTCGCACTCCAGGTAATCGCCGATCAGATACTTGGCCTGCTGCAATTCGGTCACTGATTCCCGCCGCGCCATCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCg                     <  1:1435394/150‑1 (MQ=255)
                                                                                                                                 cccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCg                     >  1:1046957/1‑150 (MQ=255)
                                                                                                                                 cccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCg                     >  2:789944/1‑150 (MQ=255)
                                                                                                                                                aGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTCCAAAggg      <  1:294275/150‑1 (MQ=255)
                                                                                                                                                  gCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTCCTCCAAAGGGCCg   >  2:565215/1‑150 (MQ=255)
                                                                                                                                                    gTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTCCAAAGGGCCGa  <  2:1609078/150‑1 (MQ=255)
                                                                                                                                                    gTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGATTGGCAGGCTGCCTTGGCCGTTTTCCTCCAAAGGGCCGa  >  1:1162284/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTCCAAAGGGCCGA  >  NC_002947/5685819‑5686116

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: