Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,318,674 T→G 15.8% A451A (GCT→GCG U712_06715 → Putative mannosyltransferase ykcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,318,6740TG15.8% 118.8 / 4.2 69A451A (GCT→GCGU712_06715Putative mannosyltransferase ykcB
Reads supporting (aligned to +/- strand):  ref base T (23/34);  new base G (10/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 9.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAAGGAAACGGTGTTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTTGGATGGAGTATCGGGGTA  >  CP006881/1318537‑1318812
                                                                                                                                         |                                                                                                                                          
cAGAAGGAAACGGTGTTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAACTTATCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGGCGctc                                                                                                                                           <  3:397061/139‑1 (MQ=255)
   aaGGAAACGGTGTTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTGTTg                                                                                                                                        >  3:389625/1‑139 (MQ=255)
     ggAAACGGTGTTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGtt                                                                                                                                      <  8:367340/139‑1 (MQ=255)
        aaCGGTGTTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCAt                                                                                                                                   >  5:104741/1‑139 (MQ=255)
           ggTGTTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTg                                                                                                                                <  5:217227/139‑1 (MQ=255)
             tgtTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTGTTGTTCATCTGtg                                                                                                                              >  3:327968/1‑138 (MQ=255)
              gtTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTAc                                                                                                                             <  5:121992/139‑1 (MQ=255)
              gtTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTAc                                                                                                                             <  3:203942/139‑1 (MQ=255)
              gtTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGGCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTAc                                                                                                                             <  1:337764/139‑1 (MQ=255)
               ttCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTAc                                                                                                                             <  7:176291/138‑1 (MQ=255)
                 cTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCg                                                                                                                           <  8:6513/138‑1 (MQ=255)
                   gggTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCaa                                                                                                                        <  2:385221/139‑1 (MQ=255)
                       tGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAAt                                                                                                                       >  4:178427/1‑136 (MQ=255)
                       tGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCGTGGTGATCTGTTCCGCAATCat                                                                                                                    >  3:275380/1‑138 (MQ=255)
                       tGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCGCTTGTTCAACTGTACCGCCAACAg                                                                                                                    >  1:134481/1‑139 (MQ=255)
                       tGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGGCGCGGTTGGTCATCAGTACCGCAAt                                                                                                                       >  8:177221/1‑136 (MQ=38)
                         cGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTTGCACCGCCGATTCCGGCTCTTGTCGGTGGGGGGTGGGTCGCTCTTTGTGAGCTGTACCGCCATCAgc                                                                                                                   >  7:1291/1‑137 (MQ=255)
                         cGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTGATCTGTAGCGGAATCata                                                                                                                   >  5:283274/1‑136 (MQ=255)
                         cGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCAt                                                                                                                                   <  5:82607/122‑1 (MQ=255)
                         cGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGa                                                                                                                   >  5:102515/1‑138 (MQ=255)
                         cGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGa                                                                                                                   >  3:309832/1‑138 (MQ=255)
                         cGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGAc                                                                                                                  <  7:220601/139‑1 (MQ=255)
                         cGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGGGGGTCGCGGTTGGTGATAGGTACCGCAAACAGAc                                                                                                                  >  8:375375/1‑139 (MQ=38)
                         cGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACAGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCAt                                                                                                                                   >  6:82607/1‑122 (MQ=255)
                              cTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCGGTTGGGGATGTGGTCCGCACTCAGACAGGcg                                                                                                             >  4:129783/1‑138 (MQ=38)
                              cTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCGCGTGTTGATCTGTACCGCAATAAGACAGGCt                                                                                                             >  8:220601/1‑139 (MQ=255)
                                 gTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTg                                                                                                            <  7:397196/137‑1 (MQ=255)
                                       aTTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGc                                                                                                      >  2:230360/1‑137 (MQ=255)
                                       aTTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCGTTTTGATCTGTTCCGGCATCAGTCAGGCTGGAAAGCAt                                                                                                    >  4:242965/1‑139 (MQ=38)
                                        ttGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGGGGGTCGCTGGTGGTGATCTGTTCCGGAATCAGGCAGGGGGGGAAGCATg                                                                                                   >  6:228205/1‑139 (MQ=37)
                                         tGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATAGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTGTGTGTCATCAGTAACGCAAGCAGACACGCAGGAAATCAttg                                                                                                  >  1:104046/1‑137 (MQ=255)
                                                     cAGTGTCGCATAATTTTTCCATCATTATTATCTCATCATGCTCGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCTCTCTTGTTCATCTGTACCGCAATCAGACAGTCTGGAAAGCATGGCTTTTGCCggg                                                                                       <  1:396001/138‑1 (MQ=255)
                                                     cAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCggg                                                                                       <  2:34383/138‑1 (MQ=255)
                                                        tgtCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGGTGATCTGTACCGCAAACAGCCAGCCAGGTAAGCCTGGCTTTTGGCGGGTGc                                                                                    >  7:211902/1‑138 (MQ=255)
                                                             cAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCGTGTCGGTGCGGGGTGGGTCGCGGTTGTTAAACTGTTCCCCCAACCGACCGGCGGGAAAGCCAGGGGTTTGGCGGGAGCCattat                                                                              >  4:35544/1‑139 (MQ=255)
                                                                aaTTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCGTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATATGTACCGCAATAAGACAGGCTGGAAAGCATTGCTTTTGCCGGGAGCAattttttc                                                                           >  8:95474/1‑137 (MQ=255)
                                                                       ccATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAAc                                                                     <  6:102515/138‑1 (MQ=255)
                                                                                ttatCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGa                                                            <  3:35544/138‑1 (MQ=255)
                                                                                ttatCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGGCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGa                                                            <  7:177221/138‑1 (MQ=255)
                                                                                     tcatcaTGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTGATCTTTGCCGCAATCAGTCACGCTGGAAAGCAGGGCGTTTTGCGGGTGGAATTATTGCGACAACTTGGTTTGTATTGtt                                                      >  5:268192/1‑139 (MQ=255)
                                                                                      catcatGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCGGTTGTTGAGCTGGTCCGGCATCCGACCGGGTGGAAAGCATGGCTTTTGCCGGGGGGAAATAATGCGTCAACTGGTTTTGTAATGttt                                                     >  4:215024/1‑139 (MQ=255)
                                                                                       atcatGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGt                                                       <  1:294465/136‑1 (MQ=255)
                                                                                       atcatGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTa                                                    <  3:129783/139‑1 (MQ=255)
                                                                                       atcatGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCGGTTGTTAATCTGTTACGCAATAAGAACGGCAGGCAAACATGGCTTTTGCCCGGCGCAAATTTTTCGACAACAGGATTTGTAATGtttt                                                    >  6:336407/1‑138 (MQ=255)
                                                                                        tcatGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTGATCTGTAAAGCAAAAACCCAGGCTGGAAAGCAAGGCTTTTGCCCGGAGCAATTATTTGGACAACAGGATTTTAATATTTTAt                                                   >  2:294465/1‑139 (MQ=255)
                                                                                          atGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTGTTGTTGATCTGTTCCGCAATCAGACCGGCTGGGAAGCATGGCTTTTGGCGGGTGCAATTATTGCGACAAATGGATTTGAATTGTTTAtta                                                 >  2:34257/1‑138 (MQ=255)
                                                                                           tGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGTTCGCTGGTGTTGATATTTTACGCAAACAGAAAGGCGTGTAAGGATGGCTTTTGCCGGGAGGAAATATTTAGACAACGTGATTTGATTTTTTTaata                                                 >  6:359503/1‑135 (MQ=255)
                                                                                             cTGTCACCGCCGTTTGCGGCTCTTGTCGGTGCGGGGGGGGTCGCTCTTGTTCATCTGTCCCGCAATCAGACAGGCTTGAAAGCATGGCTTTTGCCGGGAGTAATTATTGCGACAACTGGATTTGAATTGTTTATtctc                                               <  8:1291/138‑1 (MQ=255)
                                                                                                 cACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGGTGATATGTACCGCAATCAGACAGGCGGGAAAGCATTGCTTTTGCCGGGAGCAATTATGGCGCCAACAGGATAAGAATTGTTTTTTCTCAGcaa                                          >  6:268523/1‑139 (MQ=255)
                                                                                                      ccgATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTAc                                      >  7:211573/1‑138 (MQ=255)
                                                                                                      ccgATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCGGGTGTTGATCTGTTCCCCAATAAGACCGGCTGGAAAGCAAGGCGTTTTCCGGGTGCAAATATTGGGGCCACAGGATTTGAATTGTTTTTTATCCTAAACTAc                                      >  2:138041/1‑138 (MQ=255)
                                                                                                      ccgATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCGGGTGGTCATCTGTTACGCAATCCGGCAGGGGGGCAACCCTGGGTTTTGGCGGGGGCGAATATTTCGACAACTGGAGTTGTTTTGTTTTTTTACAGCAACtcc                                      >  6:81889/1‑136 (MQ=255)
                                                                                                                ctctTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAAt                           <  4:275380/139‑1 (MQ=255)
                                                                                                                ctctTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAAt                           <  8:211573/139‑1 (MQ=255)
                                                                                                                ctctTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAAt                           <  6:104741/139‑1 (MQ=255)
                                                                                                                ctctTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAAt                           <  4:309832/139‑1 (MQ=255)
                                                                                                                ctctTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAAt                           <  5:336407/139‑1 (MQ=255)
                                                                                                                      tCGGTGCGGGGTGGGGCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACTACTTGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGgtg                      <  2:104046/138‑1 (MQ=255)
                                                                                                                          tGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTtgga                 <  6:283274/139‑1 (MQ=255)
                                                                                                                          tGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTtgga                 <  8:211902/139‑1 (MQ=255)
                                                                                                                          tGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTtgga                 <  1:230360/139‑1 (MQ=255)
                                                                                                                          tGCGGGGTGGGGCGCTCTTGTGCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTtgga                 <  6:82238/139‑1 (MQ=255)
                                                                                                                            cGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTTTTCTCAGAAACTACAACGATCAAATCGGTGTtggatg               <  4:389625/139‑1 (MQ=255)
                                                                                                                            cGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTtggatg               <  3:150118/139‑1 (MQ=255)
                                                                                                                            cGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTtggatg               <  5:228205/139‑1 (MQ=255)
                                                                                                                            cGGGGTGGGTCGCGCTTGTTCATCTGGACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTtggatg               <  2:134481/139‑1 (MQ=255)
                                                                                                                                  gggtCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACACCGATCAAATCGGTGTTGGATGGAGTAt         <  5:268523/139‑1 (MQ=255)
                                                                                                                                  ggggCGCGCTTGTTAATCCTTACCGCAATCAGCCAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGCCAACTTGCTTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTTGGATGGAGTAt         <  6:124806/139‑1 (MQ=255)
                                                                                                                                    gtCGCTCTTGTTGATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAAATATTGCGACAACTGGATTTGAATTTGTTATTCTCAGAACATAAAACGATAAAAACGGTGTTGGTTGGAGTATCg       >  5:124806/1‑139 (MQ=255)
                                                                                                                                       gCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTTGGATGGAGTATCgggg    <  5:224971/139‑1 (MQ=255)
                                                                                                                                        ctctTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTTGGATGGAGTATCGGGGt   <  2:91095/139‑1 (MQ=255)
                                                                                                                                         tctTGTTGATCTGTACCGCCATCAGACACGCTGGAAAGCATGGCGTTTTGCGGGAGCAATTATTTCGACAACTTGATTTTAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTTGGATGGAGTATCGGGGTa  >  4:76755/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CAGAAGGAAACGGTGTTCTGGGTTGCGTGGCTTGTACCGATTGCCGGATTCTTCAGTGTCGCAGAATTTTTCCATCATTATTATCTCATCATGCTGGCACCGCCGATTGCGGCTCTTGTCGGTGCGGGGTGGGTCGCTCTTGTTCATCTGTACCGCAATCAGACAGGCTGGAAAGCATGGCTTTTGCCGGGAGCAATTATTGCGACAACTGGATTTGAATTGTTTATTCTCAGAAACTACAACGATCAAATCGGTGTTGGATGGAGTATCGGGGTA  >  CP006881/1318537‑1318812

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: