Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,995,958 A→C 8.3% intergenic (‑204/+81) PP_mr26 ← / ← rfbC C4 ncRNA/dTDP‑4‑dehydrorhamnose 3,5‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,995,9580AC 8.3% 165.6 / 3.1 57intergenic (‑204/+81)PP_mr26/rfbCC4 ncRNA/dTDP‑4‑dehydrorhamnose 3,5‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base C (3/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTAC  >  NC_002947/1995820‑1996104
                                                                                                                                          |                                                                                                                                                  
tCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTa                                                                                                                                           >  6:58584/1‑148 (MQ=255)
tCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTAt                                                                                                                                          >  8:71371/1‑149 (MQ=255)
  gagaCTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTTTTTCACACGCCACTaag                                                                                                                                         >  6:121691/1‑146 (MQ=38)
   agaCTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTTTTTCACACGCCACTaagaa                                                                                                                                       <  5:121691/149‑5 (MQ=38)
    gaCTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGc                                                                                                                                      <  8:254366/149‑1 (MQ=255)
       tAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGc                                                                                                                                      >  6:88255/1‑146 (MQ=255)
          gTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGTGAAATGGCCTACAAGGAAGGGGCAAATTATGATTATTTGCCAGACAATGCAGCGGTTGCTTCTTAGTGGGGTGTGAAAAACGCCACTATCTGCAAGAcc                                                                                                                                >  6:66348/1‑149 (MQ=38)
          gTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCAACTATCTGCAAGAcc                                                                                                                                >  6:66345/1‑149 (MQ=255)
           tctcACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGAccc                                                                                                                               >  4:113447/1‑149 (MQ=255)
             tcACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTTTTTCACACGCCACTaa                                                                                                                                          >  4:321767/1‑135 (MQ=38)
             tcACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTTTTTCACACGCCACTaa                                                                                                                                          <  3:321767/136‑2 (MQ=38)
             tcACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGc                                                                                                                                                >  4:355662/1‑130 (MQ=255)
             tcACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGc                                                                                                                                                <  3:355662/130‑1 (MQ=255)
               aCGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCg                                                                                                                              <  6:57161/146‑1 (MQ=255)
                 gTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCa                                                                                                                                              >  1:251275/1‑128 (MQ=255)
                 gTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCa                                                                                                                                              <  2:251275/128‑1 (MQ=255)
                        gtttgCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAACTGGCCTACAAGGAATGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAtgtg                                                                                                                  <  5:348228/146‑1 (MQ=255)
                         tgtgCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAt                                                                                                                     <  1:71751/145‑1 (MQ=255)
                         tgtgCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAt                                                                                                                     >  2:71751/1‑145 (MQ=255)
                              aaTGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCAc                                                                                                                                             >  7:143053/1‑116 (MQ=255)
                              aaTGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCAc                                                                                                                                             <  8:143053/116‑1 (MQ=255)
                                  gAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGAc                                                                                                         >  1:290791/1‑148 (MQ=255)
                                        aGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGAGTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAACACGCCACTATCGGCAAGACCCCGTATGGCATGTGTTCAGAGACAAGGCtt                                                                                                  >  8:377672/1‑148 (MQ=255)
                                          gcgcAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTgcg                                                                                                >  4:106826/1‑149 (MQ=255)
                                            gcAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCTGTATGGCATGTGTTCAGAGACAAGGCTgc                                                                                                 <  8:222804/146‑1 (MQ=255)
                                            gcAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTgc                                                                                                 >  7:222804/1‑146 (MQ=255)
                                            gcAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCt                                                                                              <  1:255841/149‑1 (MQ=255)
                                                gggTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCg                                                                                          >  7:74694/1‑149 (MQ=255)
                                                gggTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCg                                                                                          >  3:8570/1‑149 (MQ=255)
                                                       tttGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCt                                                                                   >  7:257421/1‑149 (MQ=255)
                                                         tGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGGTAGCGCAGccttg                                                                                  <  7:71371/148‑4 (MQ=255)
                                                           ggAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTtct                                                                                <  5:190652/148‑1 (MQ=255)
                                                               aTGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCt                                                                           <  3:113447/149‑1 (MQ=255)
                                                                       cAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGt                                                                   >  4:310393/1‑149 (MQ=255)
                                                                        aaggaaggGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGtt                                                                  <  5:88255/149‑1 (MQ=255)
                                                                         aggaaggGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTa                                                                 <  2:257584/149‑1 (MQ=255)
                                                                                 gCAAATTCTGAGTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCCGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCgca                                                           >  1:281976/1‑145 (MQ=255)
                                                                                   aaaTTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGt                                                       >  8:26087/1‑149 (MQ=255)
                                                                                         tGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCg                                                 >  3:354613/1‑149 (MQ=255)
                                                                                            ttattTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCa                                              <  2:290791/149‑1 (MQ=255)
                                                                                             tattTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGCCCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAg                                             >  5:292387/1‑149 (MQ=255)
                                                                                             tattTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTccg                                             >  6:127020/1‑147 (MQ=255)
                                                                                               tttGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGc                                            >  6:344007/1‑148 (MQ=255)
                                                                                               tttGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGc                                            <  3:310393/148‑1 (MQ=255)
                                                                                                 tGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCaa                                         >  7:101758/1‑149 (MQ=255)
                                                                                                   ccAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAAt                                        <  4:8570/148‑1 (MQ=255)
                                                                                                   ccAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATg                                       <  7:26087/149‑1 (MQ=255)
                                                                                                         cATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGAccc                                 <  3:106826/149‑1 (MQ=255)
                                                                                                              aGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGc                            >  3:101386/1‑149 (MQ=255)
                                                                                                                 ggTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTgg                         <  2:210525/149‑1 (MQ=255)
                                                                                                                   ttGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTc                       <  4:101386/149‑1 (MQ=255)
                                                                                                                      cttcttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCtt                    <  8:257421/149‑1 (MQ=255)
                                                                                                                         cttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGc                  <  5:344007/148‑1 (MQ=255)
                                                                                                                          ttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTc                       >  5:180402/1‑142 (MQ=255)
                                                                                                                          ttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTc                       <  6:180402/142‑1 (MQ=255)
                                                                                                                                gCGTGTGAAAAACTCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACaa           <  8:74694/148‑1 (MQ=255)
                                                                                                                                  gtgtGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCt        >  5:41676/1‑149 (MQ=255)
                                                                                                                                        aaaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTa                                                                 >  4:313774/1‑86 (MQ=255)
                                                                                                                                        aaaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTa                                                                 <  3:313774/86‑1 (MQ=255)
                                                                                                                                        aaaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGATGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCATTCCCAGACAAGCTTGGTAc  >  1:162502/1‑149 (MQ=255)
                                                                                                                                         aaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGc                  <  4:13123/132‑1 (MQ=255)
                                                                                                                                         aaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGc                  >  3:13123/1‑132 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
TCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTAC  >  NC_002947/1995820‑1996104

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: