Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 986,209 C→A 20.4% intergenic (+75/‑175) U712_04810 → / → U712_04815 putative oxidoreductase yhdF/putative amino acid permease YhdG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881986,2090CA20.4% 151.5 / 2.9 67intergenic (+75/‑175)U712_04810/U712_04815putative oxidoreductase yhdF/putative amino acid permease YhdG
Reads supporting (aligned to +/- strand):  ref base C (8/40);  new base A (7/3);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 1.64e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTGTTAC  >  CP006881/986079‑986322
                                                                                                                                    |                                                                                                                 
agtataTGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAAtt                                                                                                           <  7:138012/137‑1 (MQ=255)
    tatGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTCAAAGGAATCAAAGAAAAAACAACACAAGCGGGGGGTTTTTGTGTGGTGTAAAATTTAAAAAAAAAAAAATTTTGGATTTTTTTAACAAattttt                                                                                                       >  3:7114/1‑137 (MQ=255)
        aCAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTATTTAAACCTAATAAACCAAAAAACAACACAACACCTGGTTTTTTTTTTTGTGTAAAATTTAAAAAAAAAAAATATTTGTAGTTTTATAACAAAtttttttt                                                                                                    >  2:315694/1‑138 (MQ=255)
          aGGGCAGACCATTCATGTGAATGGCGGCCGTTTTAATTAAACGTAATCAACCAAAAACCAGCTCAAGCCCAGGTTTTTTGTGTGGTGAAAAATTTAAAAAAAAAAAATTTGTGGATATTTAAAACAAGTTTGGTTTAAt                                                                                                 >  8:279509/1‑139 (MQ=255)
          aGGGCAGACCATTCATGTGAATGGCGGCCGTGTTATTTAAACCGAATCAACGAAAAAACACCTCAAGAGCTGGTTTTTTTTGTGGTGAAAAAATTAAAAAAAAAAAAATTTTATTTTTTTTTAAAAATTTTTTTTTAAt                                                                                                 >  2:346191/1‑139 (MQ=255)
            ggCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAACCTAATAAACCAAAAAACACCAAAAAAGGTGGTTTTTTTTGTTGTGAAAAAATTAAAAAAAAAAAATTTGTAGATTTTTGTAAAAAATTTTTTTTTAttt                                                                                               >  2:382724/1‑139 (MQ=255)
            ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTTGGGTGAAAAATTTCAAAATAAAAAATTTGTGGATTTTTGTAAAAATTTTGTTTTCAttt                                                                                               >  6:90915/1‑139 (MQ=255)
            ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAACCAAAAACCACCAGAAGAGCGGGGTTTTTTTGTGGTGAAAAATTTAAAAATAAAAAAATTGTAGAGTTTTTTTAAAATTTTTTTTTCAttt                                                                                               >  6:351365/1‑139 (MQ=255)
            ggCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTAAACTTAATAAAAAAAAAAAAACCACAAAACATGGTGTTTTTTTTGGTTTAAAATTTAAAAAAAAAAAAATTTTATATTTTTTTTAAAAATTTTTTTTTAAtt                                                                                                >  2:380721/1‑139 (MQ=255)
            ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTGAACGTAATCAACGAAAAAACACCTCAAGAGGTGGGTTTTTTTGTGGTGATAAATTTGAAAATAAAAAATTTTTAGAGTTTTTTAAAAATTTTGTTTTGAttt                                                                                               >  3:128614/1‑139 (MQ=255)
            ggCAGACCATTCATGTGAATGGAGGCCGTTTTATTTAAAAGTAATCAACGAAAAACCAACTCAAGAGCGGTTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTTTCGATTTTTATAAAAAAATTTTTTTAAttt                                                                                               >  8:23629/1‑139 (MQ=255)
             gCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATCAACCAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTTTAGATTTTTCTAAAAATTTTTTTTTAAttat                                                                                              >  8:138012/1‑137 (MQ=255)
               cgACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                            <  3:59343/138‑1 (MQ=255)
                 aCCATTCATGTGAATGGCGGCCGTTTTATTTGAAAGGAATAAACCAAAAAACACCGCAAGAGCATGGTTTTTTTGTGGTGAATAATTTAAAAATAAAAAATTTTTGGATTTTTTTTAAAAATTTATTTTAAtgttttg                                                                                           >  2:355788/1‑136 (MQ=255)
                 aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAACCCAGCTCAAGAGCTGGTTTTTTGTGTTGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCa                                                                                          <  3:370217/139‑1 (MQ=255)
                 aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCa                                                                                          <  5:73870/139‑1 (MQ=255)
                 aCCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATCAACCAAAAAACAGCACAATAGCTGGTTTTTTTTGTTGTGTAAAATTTAAAAAAAAAAAATTTTTCGATTTTTATAACAATTTTTTTTTAAttttttg                                                                                           >  6:356816/1‑136 (MQ=255)
                 aCCATTCATGTGAATGGCGGCCGTGTTTTTTTTACGTAAAAAAACAAAAAACACCCCAAAAGCTGGGTTTTTTTTTGGTTTTAAATATAAAAAAAAAAAAATTTTTGATTTTTTTAAAAAAATTTTTTTAtttttttga                                                                                          >  6:26683/1‑136 (MQ=255)
                 aCCATTCATGTGAATGGCGGCCGTGTTATTTAAAAGTAATAAAACAAAAAACAGCCCAAGCGCGGGGTTTTTTTGTGGTGGAAAAATTAAAAAATAAAAATTTTTGTGTTTTTATAAAAAtttttttttttttttttg                                                                                           >  5:358698/1‑136 (MQ=255)
                 aCCATTCATGTGAATGGCGGCCGTGTTAATTAAAACGAATAAACCAAAAACCATCACAAGAACTGGGTTTTTTTTTTGTTAAAAAATTAAAAAAAAAAAATTTTTAGGTTTTTATAAAAAATTTTTTTTGAattttta                                                                                           >  8:297711/1‑136 (MQ=255)
                    aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTTACAATTTTGTTTTCATTTTTGCAAtt                                                                                       <  7:14057/139‑1 (MQ=255)
                       cATGTGAATTGCGGCAGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTTGTTTTTTGTGTTGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                     <  3:364535/138‑1 (MQ=255)
                       cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                     <  6:152251/138‑1 (MQ=255)
                       cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                     <  7:396690/138‑1 (MQ=255)
                       cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                     <  4:303743/138‑1 (MQ=255)
                       cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                     <  1:266372/138‑1 (MQ=255)
                       cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                     <  1:348472/138‑1 (MQ=255)
                       cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                     <  1:349057/138‑1 (MQ=255)
                       cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                     <  3:133104/138‑1 (MQ=255)
                       cATGTGAATGGCGGCAGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                     <  7:240483/138‑1 (MQ=255)
                        aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                                   <  4:128614/139‑1 (MQ=255)
                           tgtattgcgggcggTTTAATTCCAAGTAATCCACGAAAAAACAGCTCAAGAGCTGGTTTTTTTTGTTGTGTAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                                <  1:346191/133‑1 (MQ=255)
                           tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                                <  3:255277/139‑1 (MQ=255)
                           tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                                <  3:260774/139‑1 (MQ=255)
                           tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCACTTTTGCAATTTtatata                                                                                <  8:140216/139‑1 (MQ=255)
                           tgAATGGAGTCCGTTTTATATCAACGTAATAAAAGAAAAACCAGCTCAAGAGCTTGTTTTTTGTGTTGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTTCAATTTtatata                                                                                <  1:134709/139‑1 (MQ=255)
                                 gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                           <  3:345128/138‑1 (MQ=255)
                                 gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                           <  5:356816/138‑1 (MQ=255)
                                 gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                           <  5:231625/138‑1 (MQ=255)
                                    gcCGTTTTATTTCAACGTAATCAACGAAAAAACACCTCAAGAGCAGGTTTTTTGTTTGGTGAAAAATTTTAAAATAAAAAAATTTTGGATTTTTATAACAATTTTGTTTTTATTTTTGAAATTTTTTTTATAAAAtttt                                                                       >  5:114878/1‑137 (MQ=255)
                                     ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                                      <  1:315694/139‑1 (MQ=255)
                                     ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                                      <  3:402435/139‑1 (MQ=255)
                                     ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                                      <  6:146034/139‑1 (MQ=255)
                                             tttCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGgaga                                                               <  1:37855/138‑1 (MQ=255)
                                             tttCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGaa                                                              <  1:380721/139‑1 (MQ=255)
                                               tCAACGTAATCAACGATAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGc                                                                                           <  3:326015/108‑1 (MQ=255)
                                               tCAACGTAATAAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGc                                                                                           >  4:326015/1‑108 (MQ=255)
                                                   cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                           <  2:80192/120‑1 (MQ=255)
                                                   cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                           >  1:80192/1‑120 (MQ=255)
                                                   cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATATTTCTAACAATTTTGTATTCATTTTTGCAATTTTATATATAAAATTGAGTAGGAGAAAGGAAt                                                        >  8:247893/1‑139 (MQ=255)
                                                       aTCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGAtt                                                    >  3:347257/1‑139 (MQ=255)
                                                        tCAACGAAAAACCATCTCAAGTGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTAATTTTTTCAATTTTATATATAATATTGAGTAGGAGAAAGTAATGATTc                                                   <  6:154526/139‑1 (MQ=255)
                                                             gAAAACCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTCGGAGAAAGGAATGATTCGCGAc                                              <  6:131334/139‑1 (MQ=255)
                                                             gAAAAACCAGCTCAAGAGCTGTTTTTTTGTTTTGTTAAAAATTTAAAAATAAAAAATTTGTCTATTTTTCTAAAAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTCAGTTGTATAAAGTAATGATTCGCGAc                                              <  6:114878/139‑1 (MQ=255)
                                                             gAAAAACCAGATCAATAGCTGTTTTTTTTTTTGGTTAAAAATTTCAAAATAAAAAATTTGTCTATTTTTATAAAATTTTTGTTTTCTTTTTTGAAATTTTATATATAATTTTTATTAGTACAATGTAATTTTTCGCGAc                                              <  8:293440/139‑1 (MQ=255)
                                                              aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                             <  3:375426/139‑1 (MQ=255)
                                                              aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTGTCGATTTTTATAACAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                             <  7:16443/139‑1 (MQ=255)
                                                               aaaaCCAGGTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                             <  7:279509/138‑1 (MQ=255)
                                                                aaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTAATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAg                                           <  5:162692/139‑1 (MQ=255)
                                                                   ccAGCTCATTATCTGGTTTTTTGTGTGTTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACATTTTTGTTTTCTTTTTTGCAATTTTATATTTAATATTGTTTAGTAGAAAGTATTTATTCGCGACCAGTAc                                        <  4:55250/139‑1 (MQ=255)
                                                                      gCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACaaa                                     >  3:55250/1‑139 (MQ=255)
                                                                                  ggTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGaa                           >  3:279963/1‑138 (MQ=255)
                                                                                      ttttGTTTGTTTAAAAATTTCAAAATAAAAATTTTGTCTATTTTTAAAACATTTTTGTTTACATTTTTGAAATTTTTTATAATATTTTGCGTACGAGAATGGAATTATCCGCGACCAGTACAAATATAGTAGaaaaaa                       <  3:234235/138‑1 (MQ=255)
                                                                                           tgtgGTTAAAAAATTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTATGTTTTCATTTTTGCAATTTTATATATAATATTGAGTCGGATAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaat                 <  2:145379/139‑1 (MQ=255)
                                                                                               gttaaaaaTTTCAAAATAAAAATTTTGTCGATTTTTCTAACATTTTTGTTTTCATTTTTGCAATTTTAATTATAATATTGAGTAGGAGAAAGGAATTATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataat              <  7:297711/135‑1 (MQ=255)
                                                                                                  aaaaaTTTAAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTAATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAc           <  4:347257/138‑1 (MQ=255)
                                                                                                          cAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTGTTAc  >  8:237616/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                 
CCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTGTTAC  >  CP006881/986079‑986322

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