Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,049,411 G→T Q97K (CAG→AAG)  glnL ← sensory histidine kinase in two‑component regulatory system with GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,049,4110GT96.8% 169.0 / NA 63Q97K (CAG→AAG) glnLsensory histidine kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (28/33);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTA  >  CP009273/4049265‑4049553
                                                                                                                                                  |                                                                                                                                              
aaGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTtgg                                                                                                                                              <  1:694203/149‑1 (MQ=255)
      aTTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTca                                                                                                                                        >  1:701466/1‑149 (MQ=255)
          ttAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGa                                                                                                                                    >  1:298116/1‑149 (MQ=255)
              tctcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAgg                                                                                                                                <  2:531208/149‑1 (MQ=255)
                 cATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGCata                                                                                                                             >  2:156860/1‑149 (MQ=255)
                      gCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCCTAGCCGTCCCAGAAAGGATATGCGa                                                                                                                        >  1:199264/1‑149 (MQ=255)
                       ccAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATCGGCGAg                                                                                                                       >  2:401777/1‑149 (MQ=255)
                       ccAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCGGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCCTACGCTTGGCCGTCACCACAAGGATATGCgcg                                                                                                                       >  2:687385/1‑147 (MQ=255)
                         aGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGTCCGTCACAGAAAGGATATGCAAgcg                                                                                                                     >  1:421727/1‑149 (MQ=255)
                          ggCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAgcgc                                                                                                                    <  1:204506/149‑1 (MQ=255)
                                           cGGGGAGCAACGTGGTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGt                                                                                                   <  2:646664/149‑1 (MQ=255)
                                           cGGGCAGCAACCTGCTGGGCGTGCTTTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCTAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGt                                                                                                   <  1:81414/149‑1 (MQ=255)
                                              gcagcaACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCAc                                                                                                <  1:382515/149‑1 (MQ=255)
                                               cagcaACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACt                                                                                               <  2:701466/149‑1 (MQ=255)
                                                  caACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCg                                                                                            <  1:395056/149‑1 (MQ=255)
                                                  caACCTGCTGGGCGTGCTGTAGCTGGTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACCCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCg                                                                                            >  1:444118/1‑149 (MQ=255)
                                                   aaCCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGt                                                                                           <  1:81793/149‑1 (MQ=255)
                                                      ctgctgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTAt                                                                                        >  1:672207/1‑149 (MQ=255)
                                                       ggctgGGCGGGCTGTAGCTGTTCCTGACTTAGGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGGAAGGATATGCGGGCGCCCGTCGATGACCAGCGTCACTTCGTTATc                                                                                       <  2:292690/148‑1 (MQ=255)
                                                       ggcctggCGGGCTGTAGCGTTTCCTGACTTAAGCGGCGGTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGGGCGCCCGTCGATGACCAGCGTCACTTCGTTATc                                                                                       <  2:292671/144‑1 (MQ=255)
                                                          tgGGCGTGCTGTCGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGTCATCTCCAGCAGGATCATGCCGTCCGGCATACGCGTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGt                                                                                    <  1:299030/149‑1 (MQ=255)
                                                          ggggCGTGCTGTAGCTGTTCCTGACGTAACCGGCGCTGCTTATCCATCGGAGTCAGCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGt                                                                                    <  1:159744/148‑1 (MQ=255)
                                                              tgTGCTGTACCTGTTCCTGACTTAAGCGGCCCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTaaaa                                                                                <  2:613235/148‑1 (MQ=255)
                                                                 gCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCt                                                                             >  2:45107/1‑149 (MQ=255)
                                                                  cTGTAGCTGTTCCTCACTTAAGCGGCGATGGTTATCCATCAGAGCCATCGCCAGCAGGATCATGCCGTCCGGCATGCGCTTGGCCGTCACAGACAGGATATTCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCtt                                                                            >  1:543402/1‑149 (MQ=255)
                                                                    ggcgcctttCCTGCCTTAAGCGGCGCTGGGTATCCATCGGAGACAGCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACGTCGTTATCGGTAAAACCTTGc                                                                          <  2:693277/142‑1 (MQ=255)
                                                                    gcaGCCGTTCCTGACTTAAGCGCCGCTGGTTATCCATCGGGGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGCAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGc                                                                          <  1:318908/147‑1 (MQ=255)
                                                                       gCTGTTCCTGACTTAAGCGGCGCTGATTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGcccc                                                                       >  1:413747/1‑149 (MQ=255)
                                                                        cTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGcccgg                                                                      >  2:592530/1‑147 (MQ=255)
                                                                        cTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCg                                                                      <  1:30888/149‑1 (MQ=255)
                                                                               tGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAg                                                               >  1:310494/1‑149 (MQ=255)
                                                                                gACTTAAGCGGCGTTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCATAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGa                                                              <  1:549705/149‑1 (MQ=255)
                                                                                 aCTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGAc                                                             >  1:603807/1‑149 (MQ=255)
                                                                                 aCTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCAGTCACAGAAAGGATATGCGAACCACCGACCTTGCCCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAAAc                                                             >  1:21113/1‑149 (MQ=255)
                                                                                  cTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCGGTCCGGCATACGCTTGGCCGTCACAGAAGGGATATGCGGGCGCCCGGCGGTGCCCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCCCCAgcac                                                            >  2:713858/1‑146 (MQ=255)
                                                                                   ttAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCCGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCCTAGATGACCAGCGTCACTTCGGTATCGGGAAAACCTTGCCCCGCCCCCAGAcgt                                                           >  2:405756/1‑147 (MQ=255)
                                                                                      aGCGGCGGTGGTTATCCATCGGTGCCATCTCCAGCAGGATCATGCCGCCCGGCATACGCTTGGCCGTCACAGAAAGGATAGGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCt                                                        <  2:173869/149‑1 (MQ=255)
                                                                                      aGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCt                                                        <  1:650139/149‑1 (MQ=255)
                                                                                      aGCGGCGCTGGTTAGCCATCGAAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCTTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCt                                                        <  2:257701/149‑1 (MQ=255)
                                                                                       gcggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCtt                                                       <  1:659950/149‑1 (MQ=255)
                                                                                       gcggcgCTGGGTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGCCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCtt                                                       <  1:635981/149‑1 (MQ=255)
                                                                                         ggctgtGGTTCTCCTTCGGGGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGc                                                     <  2:424388/144‑1 (MQ=255)
                                                                                           cgcTGGTTCACCATCGGAGCCACCCCCCGCAGGATCATGCCGTCCGGTATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTCTCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCAt                                                   <  1:703739/149‑1 (MQ=255)
                                                                                              tGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGGCCCGCCTCCAGACTTTCTTGCATCAg                                                >  2:285429/1‑149 (MQ=255)
                                                                                                    tCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGat                                          <  2:413747/149‑1 (MQ=255)
                                                                                                        tCGGTGCCTTCTCAAGCAGGACTACGCCGCCCGGCATACGCTTGGCCGTCACAGCAAGGATATGCGAGCGCCCGTCGATCACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAttt                                      <  2:107444/149‑1 (MQ=255)
                                                                                                        tCGGCGCCATCTCCAGCAGGATCATGCCGCCCGGCATACGCTTGGCCGTCACAGAAATGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAttt                                      <  2:470696/149‑1 (MQ=255)
                                                                                                             gCCATCTCCAGCAGGATCATGCCGTCCGGCAGACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATga                                 <  2:466121/149‑1 (MQ=255)
                                                                                                              ccATCTCCAGCAGGTTCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACTATCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATgag                                <  1:285461/149‑1 (MQ=255)
                                                                                                                aTCTCCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGaa                              >  2:513400/1‑149 (MQ=255)
                                                                                                                aTCTCCAGCAGGATCATGCCGCCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGaa                              <  1:710838/149‑1 (MQ=255)
                                                                                                                  ctcCAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGt                            >  2:472909/1‑149 (MQ=255)
                                                                                                                    ccAGCAGGATCATGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAg                          <  2:389059/149‑1 (MQ=255)
                                                                                                                    caagcagGATCATGCCGGCCGGCATACGCTGGGCCGTCACAGAAAGGATAGGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCGCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAg                          <  1:248392/147‑1 (MQ=255)
                                                                                                                        caagcgtaTGCCGTCCGGCATACGCTTGGTCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGGTCa                      >  1:707642/8‑149 (MQ=255)
                                                                                                                         aagcgtatgccggacGGCATACGCTTGGCCGTCACAGAAAGGGTATGCGAGCGCCCGCCGATGAACAGCGTCACCTCGTTATCGGTAAAACCTTGCCCCGCCGCCAGACGTTCTTGCATCAGCTCGGTATTTAATGAAAAGTAGCTCaa                     >  2:301888/15‑149 (MQ=255)
                                                                                                                                tGCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTcc              >  2:573930/1‑149 (MQ=255)
                                                                                                                                 gCCGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCg             <  1:407017/149‑1 (MQ=255)
                                                                                                                                  ccGTCCGGCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCgg            >  1:427500/1‑149 (MQ=255)
                                                                                                                                        ggCATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACgg      >  1:396019/1‑149 (MQ=255)
                                                                                                                                        ggCATACGCTTGGCCGTCACAGAAAGGAGATGCGAGCGCCCGTCGATGGCCAGCGTCACTTGGGTAGCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGGGAAGTGGCTCAACAGTTCCGGTAACgg      >  2:157074/1‑149 (MQ=255)
                                                                                                                                          cATACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCCCACCAGTTCCGGTAACGgtg    >  2:198613/1‑149 (MQ=255)
                                                                                                                                           aTACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGgtgt   <  1:658978/149‑1 (MQ=255)
                                                                                                                                            tACGCTTGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCTTCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTa  <  1:694050/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTA  >  CP009273/4049265‑4049553

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