Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,281,247 Δ1 bp 6.0% coding (1662/1692 nt) PP_1118 ← recombinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,281,2470G. 6.0% 216.4 / 6.7 67coding (1662/1692 nt)PP_1118recombinase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (36/27);  new base . (2/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

TTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCA  >  NC_002947/1281100‑1281392
                                                                                                                                                   |                                                                                                                                                    
tttttGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGACTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCAAATgg                                                                                                                                                     >  2:228201/1‑149 (MQ=255)
    tGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAg                                                                                                                                                 >  2:201413/1‑149 (MQ=255)
      aCAGCGACGGGGCGGGAGTATGTCTTGCCTCCAGACAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑Ca                                                                                                                                              <  2:199923/149‑1 (MQ=255)
      aCAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑c                                                                                                                                               <  4:1377/148‑1 (MQ=255)
       cAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CAt                                                                                                                                             <  3:189864/149‑1 (MQ=255)
        agagATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACGAAACTGAATTATGAGAGCCCGCCATAAGGCGGGGCTATTATGAGTCACGTTATTTGAGAATATGTCTGGCTTTCCACATGGCTAG‑CAtg                                                                                                                                            >  6:15359/1‑149 (MQ=255)
          agaTGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CAtgt‑g                                                                                                                                         <  4:257625/149‑1 (MQ=255)
            aTGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑Ag                                                                                                                                      <  8:175610/149‑1 (MQ=255)
             tGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑Agg                                                                                                                                     <  3:97577/149‑1 (MQ=255)
               gggCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTTAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACAT‑GCTAGCCAtgt‑ggaaa                                                                                                                                     <  7:331654/148‑5 (MQ=255)
                        tATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGt                                                                                                                          <  1:201413/149‑1 (MQ=255)
                                tcATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGTGTCTATTATGAGTCACTTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCaa                                                                                                                     <  8:355208/145‑1 (MQ=255)
                                gCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGa                                                                                                                   <  4:236942/148‑1 (MQ=255)
                                      acaGAAGCCGCTTCTATGGGAAGCTAGACCTAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGtt                                                                                                            >  1:31363/1‑149 (MQ=255)
                                        aGAAGCCTCTTCTATGTTAAGCTAGACGAAACTGTATTTAGAGAGCGCGCCATAAGGCGGGTCTATTATGAGTCTGGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGtttt                                                                                                          <  7:343292/149‑1 (MQ=255)
                                            gCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATgg                                                                                                                                                     <  6:245589/105‑1 (MQ=255)
                                            gCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATgg                                                                                                                                                     >  5:245589/1‑105 (MQ=255)
                                              cGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCaaa                                                                                                    <  3:2031/149‑1 (MQ=255)
                                               cgatgtATGGGAATCTAGAACAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCaaa                                                                                                    <  1:390837/143‑1 (MQ=255)
                                                      tGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCGAAGATTCAGTTTTGGCaaa                                                                                                    >  2:390837/1‑141 (MQ=255)
                                                      tGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCtt                                                                                            >  7:102027/1‑149 (MQ=255)
                                                         gAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCa                                                                                         >  3:29224/1‑149 (MQ=255)
                                                                gACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGtt                                                                                  >  3:147649/1‑149 (MQ=255)
                                                                 aCCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGtt                                                                                  >  8:27098/1‑148 (MQ=255)
                                                                    aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAgg                                                                              >  3:383821/1‑149 (MQ=255)
                                                                    aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAgg                                                                              >  5:7844/1‑149 (MQ=255)
                                                                         gAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGcc                                                                         >  8:201283/1‑149 (MQ=255)
                                                                          aaTTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGc                                                                          <  8:102027/147‑1 (MQ=255)
                                                                             tAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCCCGTTCAGGTCGCCTTgg                                                                     >  3:214242/1‑149 (MQ=255)
                                                                              aaGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                                                                     <  2:31363/148‑1 (MQ=255)
                                                                               agagagCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGt                                                                   >  1:223119/1‑149 (MQ=255)
                                                                               agagagCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGt                                                                   >  2:207159/1‑149 (MQ=255)
                                                                                gagagCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGt                                                                   >  5:165261/1‑148 (MQ=255)
                                                                                  gagCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCaa                                                                <  6:165261/149‑1 (MQ=255)
                                                                                   gacccGCCTTATGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCAGCTCGTTCAGGTCGCCTTGGGTCAAc                                                               >  7:124088/3‑149 (MQ=255)
                                                                                   agCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAAc                                                               >  4:345552/1‑149 (MQ=255)
                                                                                     cccGCCTTATGGCGGGTCTATTATGAGTCACGTTATTTGAAAATATGTATGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGAATGTCCAAAGATTCATTTTTGGCAAAGGCGGCTTCAACCTAGTTCAGGTCGCCATGAGTCAACt                                                              >  3:74426/1‑148 (MQ=255)
                                                                                       cGCCATCAGGCTGGTCTATTATGAGTAACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTTG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGaa                                                           <  6:37710/149‑1 (MQ=255)
                                                                                       cGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTg                                                             >  5:10523/1‑147 (MQ=255)
                                                                                       cGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGaa                                                           <  7:324363/149‑1 (MQ=255)
                                                                                           aTAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTc                                                       >  2:266052/1‑149 (MQ=255)
                                                                                                   ggTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCGTTGGGTCAACTGAAGGTCTACGTACt                                               >  4:51333/1‑149 (MQ=255)
                                                                                                    gTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTc                                              <  6:7844/149‑1 (MQ=255)
                                                                                                         ttatGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCa                                          >  7:149235/1‑148 (MQ=255)
                                                                                                           atGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTc                                                                                                                         >  3:349613/1‑67 (MQ=255)
                                                                                                           atGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTc                                                                                                                         <  4:349613/67‑1 (MQ=255)
                                                                                                               gtgaCGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATc                                    <  7:355066/145‑1 (MQ=255)
                                                                                                               gTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTACCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCg                                   >  5:174408/1‑149 (MQ=255)
                                                                                                                  aCGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCCCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCt                                 <  4:214242/148‑1 (MQ=255)
                                                                                                                        tttGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAAtt                          <  7:201283/149‑1 (MQ=255)
                                                                                                                        tttGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAAtt                          <  4:29224/149‑1 (MQ=255)
                                                                                                                           gATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGaa                       >  3:381188/1‑149 (MQ=255)
                                                                                                                           gATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGaa                       >  6:298563/1‑149 (MQ=255)
                                                                                                                            ataataTGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGaa                       <  5:298563/148‑1 (MQ=255)
                                                                                                                            ataataTGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAAt                      >  4:143635/1‑149 (MQ=255)
                                                                                                                              aataTGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAAtct                    >  5:268288/1‑149 (MQ=255)
                                                                                                                              aataTGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAAtct                    >  5:360531/1‑149 (MQ=255)
                                                                                                                               atatGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCtt                                <  6:123000/136‑1 (MQ=255)
                                                                                                                               atatGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCtt                                >  5:123000/1‑136 (MQ=255)
                                                                                                                                 atGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTGGTCATACATCGCTTCAAATTGATTCTCCt                 >  5:320076/1‑149 (MQ=255)
                                                                                                                                 atGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCt                 >  7:235056/1‑149 (MQ=255)
                                                                                                                                    tCTGGCTTTCCTCATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTa               <  1:228201/148‑1 (MQ=255)
                                                                                                                                    tCTGGCTTTCCACATGGCTAG‑CATGTGG‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTa               >  1:3843/1‑149 (MQ=255)
                                                                                                                                    tCTGGCTTTCCACATGGCTAG‑CATGT‑GTAGGGGGTAGAGTGTCCAAAGATTCAGTTTTGGCAAAGGCGGGTTCCCCATCGTTCAGGTGGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTa               >  1:217392/1‑149 (MQ=255)
                                                                                                                                      tGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCCGCTTCCACCTCGTTCACGTCGCCTTGGGTCAACTGAAGGTCTAAGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTTCCt            >  6:123791/1‑149 (MQ=255)
                                                                                                                                          tttCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGt        >  2:66861/1‑149 (MQ=255)
                                                                                                                                          gcctcacaT‑GCTAGCCATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCgg         >  6:133146/5‑148 (MQ=255)
                                                                                                                                          gcctcacaT‑GCTAGCCATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCgg         <  5:133146/144‑1 (MQ=255)
                                                                                                                                          gcctcacaT‑GCTAGCCATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGt        >  5:252575/5‑149 (MQ=255)
                                                                                                                                                caTGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCa  <  4:381188/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCA  >  NC_002947/1281100‑1281392

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