Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 867,829 G→A 5.4% D312N (GAC→AAC)  U712_04205 → Sensor histidine kinase yfiJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881867,8290GA 5.4% 253.0 / 2.2 74D312N (GAC→AAC) U712_04205Sensor histidine kinase yfiJ
Reads supporting (aligned to +/- strand):  ref base G (42/28);  new base A (2/2);  total (44/30)
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.62e-01

TTTGCGGTCAGCGGTGATCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGATCAAG  >  CP006881/867691‑867964
                                                                                                                                          |                                                                                                                                       
tttGCGGTCAGCGGTGATCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAAC‑GTGCAGGAAGCACTGACCACTGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGa                                                                                                                                        <  3:241428/139‑1 (MQ=255)
  tGCGGGCAGCGGTGATCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGATGCACTGACCAATCCAAAACGGCACGGCGGCGCGGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGa                                                                                                                                        <  7:301411/138‑1 (MQ=255)
        cAGCGGTGATCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCAt                                                                                                                                  <  8:201653/138‑1 (MQ=255)
        cAGCGGTGATCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCAt                                                                                                                                  >  7:201653/1‑138 (MQ=255)
            ggTGATCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGcc                                                                                                                             <  1:155671/139‑1 (MQ=255)
              tGATCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCtt                                                                                                                           <  7:79526/139‑1 (MQ=255)
              tGATCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCtt                                                                                                                           <  3:325420/139‑1 (MQ=255)
                aTCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGCCCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGt                                                                                                                         >  7:248375/1‑139 (MQ=255)
                                gCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGc                                                                                                         <  8:248375/139‑1 (MQ=255)
                                gCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGc                                                                                                         <  6:143885/139‑1 (MQ=255)
                                   gTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGAACAGCATCAGCCTTGTCATTAAGGATGACGGCa                                                                                                        >  3:116201/1‑137 (MQ=255)
                                   gTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGAACAGCATCAGCCTTGTCATTAAGGATGACGGCa                                                                                                        >  3:121221/1‑137 (MQ=255)
                                            tCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCt                                                                                             >  2:70129/1‑139 (MQ=255)
                                              cGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGa                                                                                                                                        <  6:299564/94‑1 (MQ=255)
                                              cGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGa                                                                                                                                        >  5:299564/1‑94 (MQ=255)
                                                aCCCTGATCAGAACGGTGCAGGAAGCTCTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACACCGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAgg                                                                                         >  3:222645/1‑139 (MQ=255)
                                                aCCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAgg                                                                                         >  8:301521/1‑139 (MQ=255)
                                                    tGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGgcagc                                                                                     >  3:168778/1‑139 (MQ=255)
                                                      aTCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGa                                                                                                                                        >  3:104270/1‑86 (MQ=255)
                                                      aTCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGa                                                                                                                                        <  4:104270/86‑1 (MQ=255)
                                                      aTCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCAtt                                                                                                                     >  4:339900/1‑105 (MQ=255)
                                                      aTCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCAtt                                                                                                                     <  3:339900/105‑1 (MQ=255)
                                                       tCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACt                                                                                  <  7:49388/139‑1 (MQ=255)
                                                                tGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGc                                                                                                   >  6:970/1‑113 (MQ=255)
                                                                tGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGc                                                                                                   <  5:970/113‑1 (MQ=255)
                                                                tGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTgg                                                                         >  4:134170/1‑139 (MQ=255)
                                                                tGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTgg                                                                         >  5:206482/1‑139 (MQ=255)
                                                                tGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTgg                                                                         >  5:30865/1‑139 (MQ=255)
                                                                tGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTgg                                                                         >  2:308302/1‑139 (MQ=255)
                                                                     gAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTctct                                                                    >  2:53408/1‑139 (MQ=255)
                                                                       aGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCt                                                                                                                            <  3:210721/81‑1 (MQ=255)
                                                                       aGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCt                                                                                                                            >  4:210721/1‑81 (MQ=255)
                                                                         cACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGaa                                                                >  3:281049/1‑139 (MQ=255)
                                                                          aCTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGCCGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGaa                                                                >  4:231277/1‑138 (MQ=255)
                                                                          aCTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGaa                                                                >  1:225897/1‑138 (MQ=255)
                                                                          aCTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGaa                                                                >  4:107961/1‑138 (MQ=255)
                                                                          aCTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGaa                                                                >  6:119141/1‑138 (MQ=255)
                                                                          aCTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGaa                                                                >  5:132605/1‑138 (MQ=255)
                                                                                   tGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAg                                                      >  7:157560/1‑139 (MQ=255)
                                                                                   tGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAg                                                      >  4:98034/1‑139 (MQ=255)
                                                                                   tGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAg                                                      >  3:322857/1‑139 (MQ=255)
                                                                                       aaaCGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCggg                                                  <  1:142891/139‑1 (MQ=255)
                                                                                        aaCGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCggg                                                  >  5:301816/1‑138 (MQ=255)
                                                                                            gcacggcGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCGGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCg                                             <  6:206482/139‑1 (MQ=255)
                                                                                             cacggcGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCgg                                            >  3:99644/1‑139 (MQ=255)
                                                                                             cacggcGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCgg                                            >  7:327776/1‑139 (MQ=255)
                                                                                               cggcggCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCCGCGGAg                                          >  5:8303/1‑139 (MQ=255)
                                                                                                     cgccgcAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGc                                    >  5:72977/1‑139 (MQ=255)
                                                                                                              cTGTTCGATAAAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATc                           <  7:231863/139‑1 (MQ=255)
                                                                                                                        cAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGaaa                 <  4:99644/139‑1 (MQ=255)
                                                                                                                           cTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCg                                             >  4:160739/1‑108 (MQ=255)
                                                                                                                           cTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCg                                             <  3:160739/108‑1 (MQ=255)
                                                                                                                           cTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGc               >  3:62869/1‑138 (MQ=255)
                                                                                                                           cTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGc               >  1:722/1‑138 (MQ=255)
                                                                                                                           cTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGc               >  2:251038/1‑138 (MQ=255)
                                                                                                                              gCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATc                           <  5:160155/123‑1 (MQ=255)
                                                                                                                              gCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATc                           >  6:160155/1‑123 (MQ=255)
                                                                                                                              gCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGaa            <  3:98034/138‑1 (MQ=255)
                                                                                                                              gCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGaa            <  1:308302/138‑1 (MQ=255)
                                                                                                                              gCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGaa            <  5:119141/138‑1 (MQ=255)
                                                                                                                                cTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACgg         >  6:63378/1‑139 (MQ=255)
                                                                                                                                   cacaACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGAt      >  7:323591/1‑139 (MQ=255)
                                                                                                                                   cacaACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGAt      >  7:73209/1‑139 (MQ=255)
                                                                                                                                   cacaACGAACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGAt      <  4:121221/139‑1 (MQ=255)
                                                                                                                                   cacaACGAACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGAt      <  4:116201/139‑1 (MQ=255)
                                                                                                                                    acaACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGATc     <  1:141076/139‑1 (MQ=255)
                                                                                                                                     caACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATcc                                                                                              >  2:245456/1‑49 (MQ=255)
                                                                                                                                     caACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATcc                                                                                              <  1:245456/49‑1 (MQ=255)
                                                                                                                                     caACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGATCa    <  1:70129/139‑1 (MQ=255)
                                                                                                                                      aaCGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGATCaa   <  2:722/139‑1 (MQ=255)
                                                                                                                                        cGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACggg        <  7:149191/132‑1 (MQ=255)
                                                                                                                                        cGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACggg        >  8:149191/1‑132 (MQ=255)
                                                                                                                                        cGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGATCAAg  >  2:267545/1‑138 (MQ=255)
                                                                                                                                        cGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGATCAAg  >  6:259627/1‑138 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TTTGCGGTCAGCGGTGATCCTGCAATCATCCCGCTGTCTCTCCATCCGACCCTGATCAGAACGGTGCAGGAAGCACTGACCAATGCCAAACGGCACGGCGGCGCCGCAGCCTGTTCGATACAGCTTGCCTGCACAACGGACAGCATCAGCCTTGTCATTAAGGATGACGGCAAAGGCAATCCTGAGGCAGCACTTGGTTTTGGTCTCTTGAATATGAAAAAGCGGGCAGCGGAGCATGGCGGAATGATCCGCTTTGAAAGCGAACGGGATCAAG  >  CP006881/867691‑867964

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