Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,508,400 C→G 5.3% R93R (CGC→CGG)  ftsI → Peptidoglycan D,D‑transpeptidase FtsI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,508,4000CG 5.3% 255.9 / 4.2 75R93R (CGC→CGG) ftsIPeptidoglycan D,D‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base C (28/42);  new base G (2/2);  total (31/44)
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.76e-01

GCAGCCTGCGTCACATCCCCATCCCGGCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGTGCTCGACCT  >  ACB122/1508267‑1508534
                                                                                                                                     |                                                                                                                                      
gCAGCCTGCGTCACATCCCCATCCCGGCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCt                                                                                                                                       <  7:26217/135‑1 (MQ=255)
     cTGCGTCACATCCCCATCCCGGCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGcc                                                                                                                                   >  4:228733/1‑134 (MQ=255)
           cacaTCCCCATCCCGGCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCCACCCCAAAGAGCTGCAGGCCTCCAAGGAGCGCTGGCCGCAATTg                                                                                                                            >  1:119372/1‑135 (MQ=255)
           cacaTCCCCATCCCGGCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTg                                                                                                                            >  5:369152/1‑135 (MQ=255)
            acaTCGCCATCCCGGCGCACCGCGGCTGGATCACCGACCCCAAGGGTGAGCCGTTGGCAGAGAGTACCCCGGGCACCACGTTGTGGGCCAACCCCACAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCATTTgg                                                                                                                           <  7:268347/135‑1 (MQ=255)
             catccccatcccGGCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGc                                                                                                                          <  3:78474/135‑1 (MQ=255)
                 cccatcccGGCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCcgct                                                                                                                      >  8:229449/1‑135 (MQ=255)
                  ccatcccGGCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCcgctg                                                                                                                     <  7:260694/135‑1 (MQ=255)
                         ggCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTggg                                                                                                              <  3:228733/135‑1 (MQ=255)
                           cgcACCGCGGCTTGATCACCGCCCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGc                                                                                                            <  5:324207/135‑1 (MQ=255)
                           cgcACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGc                                                                                                            <  5:326923/135‑1 (MQ=255)
                             cACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGGTGGCCGCAATTGGCCGCTGCGCTGGGGCa                                                                                                           >  5:394896/1‑134 (MQ=255)
                             cACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGGTGGCCGCAATTGGCCGCTGCGCTGGGGCa                                                                                                           >  5:387981/1‑134 (MQ=255)
                             cACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCa                                                                                                           >  8:207361/1‑134 (MQ=255)
                             cACCGCGGCTTGATCACCGACCGCAACGGTGAGCAGTTGGCCGTGACTACCCCCGTCAACCCCTTGTGGGGCCACCGCAAAGAGATGCAAGACTTCAACGACCGAAGGCCGCAATTGCCCGCTGCCCTGGGGCa                                                                                                           >  2:223555/1‑134 (MQ=255)
                                cgcgGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAAc                                                                                                       <  5:327342/135‑1 (MQ=255)
                                cgcgGATTGATCACCGACCGCACCGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTACAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAAc                                                                                                       <  1:223555/135‑1 (MQ=255)
                                   ggCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCg                                                                                                    >  8:149523/1‑135 (MQ=255)
                                    gCTTGATCACCGACCGCACCTGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCGGCGCTGGGGCAGAACCCgc                                                                                                   <  2:63874/135‑1 (MQ=255)
                                                  cgCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAGCCCGCAGCAGTTGATCGAg                                                                                     >  1:208272/1‑135 (MQ=255)
                                                   gCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAgc                                                                                    <  3:146185/135‑1 (MQ=255)
                                                     aaCGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAgcgc                                                                                  >  1:192718/1‑135 (MQ=255)
                                                           gAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGGTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGAcc                                                                            <  6:387981/135‑1 (MQ=255)
                                                           gAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGGTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGAcc                                                                            <  6:394896/135‑1 (MQ=255)
                                                               cGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAAAACCCGCAGCAGTTGATCGAGCGCCTGACCca                                                                          >  2:415240/1‑133 (MQ=255)
                                                                gTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCcagca                                                                       <  5:312534/135‑1 (MQ=255)
                                                                     ccGTGAGTACCCCGGTCACCACGTTGGGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCGGGCCGCAATTGGCCGCGGCGCTGGGGCAGAACCCGCAGCAGTTGATCGACCGCCTGACCCAGCAGGCca                                                                  <  3:87852/135‑1 (MQ=255)
                                                                      cGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCca                                                                  >  7:8604/1‑134 (MQ=255)
                                                                        tGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCcaaca                                                               <  1:272678/135‑1 (MQ=255)
                                                                        tGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCcaaca                                                               <  3:215607/135‑1 (MQ=255)
                                                                          aGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACaaa                                                             <  7:33051/135‑1 (MQ=255)
                                                                           gTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCcaacaa                                                              >  7:418595/1‑133 (MQ=255)
                                                                           gTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACaaa                                                             <  8:164530/134‑1 (MQ=255)
                                                                           gTACCCCGGGCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACaaa                                                             <  3:45414/134‑1 (MQ=255)
                                                                               cccGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCGGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGt                                                         <  7:149523/134‑1 (MQ=255)
                                                                                   gTCACGACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATc                                                    <  6:79630/135‑1 (MQ=255)
                                                                                    tCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCt                                                   <  3:226360/135‑1 (MQ=255)
                                                                                     caccacGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTa                                                  <  8:361462/135‑1 (MQ=255)
                                                                                      accacGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTAc                                                 >  1:180036/1‑135 (MQ=255)
                                                                                      accacGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGGTGATCGAGCGCCTGACCCAGCAGGCCAACACAGAGTTCATCTAc                                                 >  5:394300/1‑135 (MQ=255)
                                                                                         acGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCt                                                                                                                                       <  1:186216/46‑1 (MQ=255)
                                                                                         acGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCt                                                                                                                                       >  2:186216/1‑46 (MQ=255)
                                                                                          cGTTGTGGGCCAACCCCACAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGGCAACAAAGAGTTCATCTACCTgg                                             <  7:236772/135‑1 (MQ=255)
                                                                                             tgtgGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTcc                                          >  5:158096/1‑135 (MQ=255)
                                                                                                gggCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCggg                                       <  4:216081/135‑1 (MQ=255)
                                                                                                gggCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCggg                                       >  8:142961/1‑135 (MQ=255)
                                                                                                 ggCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCgggg                                      <  2:119372/135‑1 (MQ=255)
                                                                                                   ccAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGc                                    <  7:229449/135‑1 (MQ=255)
                                                                                                   ccAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGc                                    >  1:450073/1‑135 (MQ=255)
                                                                                                   ccAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGc                                    <  6:369152/135‑1 (MQ=255)
                                                                                                     acccccAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAgc                                                                                    >  2:99678/3‑85 (MQ=255)
                                                                                                     aaCCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAgc                                                                                    <  1:99678/85‑1 (MQ=255)
                                                                                                     aaCCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTg                                  >  3:173508/1‑135 (MQ=255)
                                                                                                       ccccAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGAc                                <  2:321331/135‑1 (MQ=255)
                                                                                                         ccAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCcagca                                                                       <  4:22778/94‑1 (MQ=255)
                                                                                                         ccAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCcagca                                                                       >  3:22778/1‑94 (MQ=255)
                                                                                                           aacgagaTGCAGGCCTCCAAGGCGCGCGGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGATCATCTACCTGGTCCGGGGGCTGAccccc                            <  1:232040/132‑1 (MQ=255)
                                                                                                           aaaGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAAAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGAccccc                            <  1:415240/135‑1 (MQ=255)
                                                                                                           aaaGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCGGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGAccccc                            <  2:449518/135‑1 (MQ=255)
                                                                                                                    cAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGCACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATATACCTGGTCCGGGGGCTGACCCCCGAGCAGGGc                   >  8:314545/1‑135 (MQ=255)
                                                                                                                    cAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGCGCAgcg                    >  3:203194/1‑132 (MQ=255)
                                                                                                                    cAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACATGGTCCGGGGGCTGCCCCCCGAGCAGGGc                   >  1:181670/1‑135 (MQ=255)
                                                                                                                     aGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAggg                    >  1:186521/1‑133 (MQ=255)
                                                                                                                     aGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGcc                  >  2:298389/1‑135 (MQ=255)
                                                                                                                      ggCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCa                 <  7:340530/135‑1 (MQ=255)
                                                                                                                          tCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCAc             >  7:406438/1‑135 (MQ=255)
                                                                                                                            cAATGAGCGCTGGACGCTATTGGCAGCTGAGCTGGGGCAGAACCCGCATCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGt           <  8:399039/135‑1 (MQ=255)
                                                                                                                            cAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGt           <  7:323880/135‑1 (MQ=255)
                                                                                                                             aaGGAGCGCTGGCCGCAATTGGCCGCTGTGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGt           <  6:410268/134‑1 (MQ=255)
                                                                                                                             aaGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGt           <  5:169707/134‑1 (MQ=255)
                                                                                                                             aaGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGt           <  2:344877/134‑1 (MQ=255)
                                                                                                                              aGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGTGc         >  3:459805/1‑135 (MQ=255)
                                                                                                                               ggAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCg                           <  7:234272/116‑1 (MQ=255)
                                                                                                                               ggAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCg                           >  8:234272/1‑116 (MQ=255)
                                                                                                                                gAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACTCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGTGCTc       <  2:347216/135‑1 (MQ=255)
                                                                                                                                     cTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGTCAAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGTGCTCGACCt  <  5:86856/135‑1 (MQ=255)
                                                                                                                                     cTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGTGCTCGACCt  <  1:10113/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
GCAGCCTGCGTCACATCCCCATCCCGGCGCACCGCGGCTTGATCACCGACCGCAACGGTGAGCCGTTGGCCGTGAGTACCCCGGTCACCACGTTGTGGGCCAACCCCAAAGAGATGCAGGCCTCCAAGGAGCGCTGGCCGCAATTGGCCGCTGCGCTGGGGCAGAACCCGCAGCAGTTGATCGAGCGCCTGACCCAGCAGGCCAACAAAGAGTTCATCTACCTGGTCCGGGGGCTGACCCCCGAGCAGGGCCAGCACGTGCTCGACCT  >  ACB122/1508267‑1508534

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: ‑