Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,304,242 C→A 11.6% T240K (ACA→AAA)  U712_16995 → Arabinose metabolism transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,304,2420CA11.6% 174.1 / 2.0 69T240K (ACA→AAA) U712_16995Arabinose metabolism transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (35/25);  new base A (2/6);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCTCTCGGCCACACGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAAGGTGATTGATATGCTGAG  >  CP006881/3304112‑3304377
                                                                                                                                  |                                                                                                                                       
cTTTCTCTCGGCCACACGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTTAAACGGATTAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACaa                                                                                                                                 <  1:26261/139‑1 (MQ=255)
   tctctcGGCCACACGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAAc                                                                                                                                <  8:54212/137‑1 (MQ=255)
   tctctcGGCCACACGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAAc                                                                                                                                <  8:332120/137‑1 (MQ=255)
    ctctcGGCCACACGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACgg                                                                                                                              >  6:359043/1‑138 (MQ=255)
          gCCACACGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTGCTTCTCCGGATATGATCGTGACATTTACAACGGAAGaaa                                                                                                                        >  1:313550/1‑138 (MQ=255)
          gCCACACGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGCTATGATGGTGACTTTTACAACGGAAGcaa                                                                                                                        >  1:313551/1‑138 (MQ=255)
              cacGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGaa                                                                                                                   >  6:168342/1‑139 (MQ=255)
               acGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAAt                                                                                                                  >  5:316400/1‑139 (MQ=255)
                  cATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCaa                                                                                                               <  5:359043/139‑1 (MQ=255)
                  attatGATGGTTATTTTTAAAGCTGATGACACACATGGCGTGAAACGTATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTACTTCTCCGGATATGATCGTGACATTTACACGGGAAGAAAAAGAATCaa                                                                                                               <  5:302591/137‑1 (MQ=255)
                     atgatgGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATATACAACGGAAGAAAAAGAATCAAAAc                                                                                                            >  3:264808/1‑139 (MQ=255)
                         tgGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTGTCCTACTCCCGATATGAACTTTAAAGTTACAAAGGAACAAAAAGAAAAAAAACTTcg                                                                                                        >  3:214215/1‑138 (MQ=255)
                           ggTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTGCGTCTCCGGATATTATAGTGTCATTTTAAACGGAAATAAAAGAATCAAAAAGTCTgg                                                                                                      >  2:185947/1‑139 (MQ=255)
                           ggTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGCCATTTACAACGGAAGAAAAAGAATCAAAACTTCTg                                                                                                       >  3:87082/1‑138 (MQ=255)
                           ggTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTgg                                                                                                      >  2:368449/1‑139 (MQ=255)
                           ggTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTg                                                                                                       >  8:1040/1‑138 (MQ=255)
                           ggTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCAGGATATGCTCGTGACATTTACAACGGAAGAAAAAGCATAAAAACGTCTgg                                                                                                      >  6:44931/1‑139 (MQ=255)
                             tATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACtt                                                                                                          >  8:39337/1‑133 (MQ=255)
                             tATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACtt                                                                                                          <  7:39337/133‑1 (MQ=255)
                             tATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTAAAACGGAAGAAAAAGAATCAAAACTTCTg                                                                                                       >  6:344723/1‑136 (MQ=255)
                               tttttAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTGCTTGTCCGGATATGATCGTGACATTTACAACGGAACAAAAAGAATCAAAACTTATGgtga                                                                                                   >  2:140836/1‑138 (MQ=255)
                                   tAACGCTGATGACACACACGGCGGTACCCGGATGAACGGTTTTATACAGGCGCACCGGGCGCGTGAGTTGTTTCCTTCTCGGGATATGCTCGGGACATTTACAACGGAAGAAAAACAATCAAAACTTCTGGAGAAAGt                                                                                               <  4:87082/138‑1 (MQ=255)
                                     atgCTGATTACCCACTTGGCGTGAACCTTTTTAACGTTTTTATACAGGCGCACCGGGAGGGTGTGTTGTTTCTTTCTCCGTATATGAGCGTCACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTaaa                                                                                            <  2:313550/137‑1 (MQ=255)
                                       gCTGATGACACACACGGCGTGAACCGGATGAACGGATTTATACCGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAg                                                                                          <  7:95745/139‑1 (MQ=255)
                                           atgaCACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCcac                                                                                      >  7:227112/1‑139 (MQ=255)
                                           atgaCACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCca                                                                                       >  1:58520/1‑138 (MQ=255)
                                           atgaCACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCca                                                                                       >  2:128690/1‑138 (MQ=255)
                                           atgaCACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCca                                                                                       >  1:20895/1‑138 (MQ=255)
                                           atgaCACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAAACAAAACTTCTGGAGAAAGTAAAAGCca                                                                                       >  3:46339/1‑138 (MQ=255)
                                           atgaCACACAAGGCGTGAAACGGATGAACGGATTAATACACGCGCACCGGGAGCGTGAGTTGTTTCCTTATCCGGATATGATCCTGAAATATACAACGGAAGAAAAAGAATCAAAAATTCTGGAGAAAGTAAAAGCca                                                                                       >  6:186129/1‑138 (MQ=255)
                                                     aGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTTGAGAAATTAAAAGCCACACTGGAGaa                                                                             >  4:207732/1‑138 (MQ=255)
                                                     aGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGaa                                                                             >  3:159046/1‑138 (MQ=255)
                                                      ggCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGaaa                                                                            >  7:305278/1‑138 (MQ=255)
                                                            aaaCGTATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGTATATGATCGTTAAATTTACAACGCAACAAAAAGAATAAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACagc                                                                      <  1:185947/138‑1 (MQ=255)
                                                            aaaCGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACagc                                                                      >  4:294675/1‑138 (MQ=255)
                                                            aaaCGGATCAACGGTTTTAAACAGGCCCCCCGGGTGCGGGAGTTGTTTCCTTCTCCGTATATGTGCGTGAAATTTACAACGGAAGAAAAACAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACagc                                                                      <  1:368449/138‑1 (MQ=255)
                                                              aCGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACagcaag                                                                   <  4:159046/139‑1 (MQ=255)
                                                                 gATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGgaga                                                                                                   <  7:22872/104‑1 (MQ=255)
                                                                 gATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGgaga                                                                                                   >  8:22872/1‑104 (MQ=255)
                                                                                    gcgcACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGga                                                                                                     >  4:261699/1‑83 (MQ=255)
                                                                                    gcgcACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGga                                                                                                     <  3:261699/83‑1 (MQ=255)
                                                                                    gcgcACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCttt                                             >  8:72854/1‑139 (MQ=255)
                                                                                    gcgcACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTAAAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAATCCACACCGGAGAAAAACAGCAAGCACATGCCGACCGCCATTCttt                                             >  5:338956/1‑139 (MQ=255)
                                                                                    gcgcACCGGGAGCGTGAGTTGTTTCCTTCGCCGGATATGATCGTGACAGTTTCAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTATAAGCGACACTGGAGAAAAACAGCAAGCACATGCCGGCAGCCATTCttt                                             >  7:365772/1‑139 (MQ=255)
                                                                                        aCCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTta                                         >  7:220496/1‑139 (MQ=255)
                                                                                              aGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGa                                   <  2:20895/139‑1 (MQ=255)
                                                                                              aGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGa                                   <  8:220496/139‑1 (MQ=255)
                                                                                              aGCGTGAGTTGTTTCATTCTCCGGTTATGATCGTTACATTTACAACGGAAAAAAAAGAATAAAAACTTCTGCAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGa                                   <  4:46339/139‑1 (MQ=255)
                                                                                              aGCGTGAGGTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACCGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGa                                   >  3:23101/1‑139 (MQ=255)
                                                                                                    aGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAAt                             >  3:76690/1‑139 (MQ=255)
                                                                                                     gtttttTCATTCTCCGGATCTGATTGTTAAAATAACACCGTAAGAAAAATAATCAAACCTTCTGGAGAAAGTAAACGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAAt                             <  3:312938/138‑1 (MQ=255)
                                                                                                     gttgttTCCTTCTCGGGATATCATCGTTAAATTAACAACGTAATAAAAATAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAAt                             <  1:140836/138‑1 (MQ=255)
                                                                                                        gttTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGt                                                                                               <  6:168443/69‑1 (MQ=255)
                                                                                                        gttTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGt                                                                                               >  5:168443/1‑69 (MQ=255)
                                                                                                         tttCCTTCTCCGGATATGAGCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTgcg                         <  5:59923/138‑1 (MQ=255)
                                                                                                            ccTTCTCCTGATATTATCGTTACATTTACAACGTAACAAAAATAATCAAAAATTCTGTAGAAATTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGa                     <  5:344723/139‑1 (MQ=255)
                                                                                                            ccTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTg                      <  3:294675/138‑1 (MQ=255)
                                                                                                             cTTCTCCGGATATGTTCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGa                     <  8:305278/138‑1 (MQ=255)
                                                                                                               tctcCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCAATCTTTGTTATAACGATGAAATTGCGCTGAAgg                  >  1:312572/1‑139 (MQ=255)
                                                                                                                  cggtatatGATCGTAATTTTTACAACGTAACAAAATTATTCAAAAATTCTGGTGAAATTAAAAGCCACACTGGAGAAAAACATCAAGCCCATGCCGACAGCCATTCTTTGTTATAACGATGAATTTGCGCTGAAGGtga               <  5:186129/135‑1 (MQ=255)
                                                                                                                   cGGTTATGATCTTTAAATTTACAACGGAAGAAAATGATTCAAAACTTCTTGAGAATGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAATGt                 <  6:130829/136‑1 (MQ=255)
                                                                                                                   cGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACCCTGGAGCAAAACAGCAAGCACATGCCGCCAGCCATTCTTTGTTATACCGATGAAATTGCGCTGAAGGtcat              >  7:70143/1‑139 (MQ=255)
                                                                                                                   cGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAAGGtga               <  7:72854/138‑1 (MQ=255)
                                                                                                                     tatttGATCGTTAAATTTAAAACGTAATAAAAATAATAAAAACTTCTGTATAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAAGGtgattg            <  4:264808/135‑1 (MQ=255)
                                                                                                                         tGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAAGGTGATTGata         <  3:78575/138‑1 (MQ=255)
                                                                                                                            tCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAAGGTGATTGATATGCt     <  7:277567/139‑1 (MQ=255)
                                                                                                                              gTGACATTTACAACGGAATAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAAGGTGATTGATATGCTGa   <  8:304264/139‑1 (MQ=255)
                                                                                                                              gTGACATTTACAACGGAACAAAAACAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAAGGTGATTGATATGCTGa   <  4:51881/139‑1 (MQ=255)
                                                                                                                               tGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAAGGTGATTGATATGCTGAg  >  7:266753/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                       
CTTTCTCTCGGCCACACGCATATGATGGGTATTTTTAAAGCTGATGACACACAAGGCGTGAAACGGATGAACGGATTTATACAGGCGCACCGGGAGCGTGAGTTGTTTCCTTCTCCGGATATGATCGTGACATTTACAACGGAAGAAAAAGAATCAAAACTTCTGGAGAAAGTAAAAGCCACACTGGAGAAAAACAGCAAGCACATGCCGACAGCCATTCTTTGTTATAACGATGAAATTGCGCTGAAGGTGATTGATATGCTGAG  >  CP006881/3304112‑3304377

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