Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,281,252:1 +C 6.0% coding (1657/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,2521.C 6.0% 264.5 / 11.5 67coding (1657/1692 nt)PP_1118recombinase‑like protein
Reads supporting (aligned to +/- strand):  ref base . (36/27);  new base C (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 = 6.81e-01

TGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCAGGTAGAA  >  NC_002947/1281104‑1281399
                                                                                                                                                     |                                                                                                                                                     
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)
                                          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                                                                                     >  5:7844/1‑149 (MQ=255)
                                                                aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAgg                                                                                     >  3:383821/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                                                                          >  2:207159/1‑149 (MQ=255)
                                                                           agagagCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGt                                                                          >  1:223119/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                                 <  4:29224/149‑1 (MQ=255)
                                                                                                                    tttGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAAtt                                 <  7:201283/149‑1 (MQ=255)
                                                                                                                       gATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGaa                              >  6:298563/1‑149 (MQ=255)
                                                                                                                       gATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGaa                              >  3:381188/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:360531/1‑149 (MQ=255)
                                                                                                                          aataTGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAAtct                           >  5:268288/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)
                                                                                                                                                  tAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCAGGTag    >  4:127751/1‑148 (MQ=255)
                                                                                                                                                   aG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCAGGTagaa  <  4:86892/149‑1 (MQ=255)
                                                                                                                                                    g‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCAGGTagaa  <  6:174408/148‑1 (MQ=255)
                                                                                                                                                    g‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCAGGTagaa  <  6:10523/148‑1 (MQ=255)
                                                                                                                                                    g‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCAGGTagaa  <  1:66861/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAG‑CATGT‑G‑AGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCTACGTACTCGTCATACATCGCTTCAAATTGAATCTCCTTACCTCGGTGCATCAGGTAGAA  >  NC_002947/1281104‑1281399

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