Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,281,194 A→T 9.3% intergenic (+20/+23) PP_5742 → / ← PP_1118 membrane protein/recombinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,281,1940AT 9.3% 211.9 / 9.6 75intergenic (+20/+23)PP_5742/PP_1118membrane protein/recombinase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (30/38);  new base T (3/4);  total (33/42)
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.06e-01

GATGGGGTATGACTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCT  >  NC_002947/1281048‑1281340
                                                                                                                                                  |                                                                                                                                                  
gATGGGGTATGACTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAgg                                                                                                                                                  >  4:358313/1‑149 (MQ=255)
  tattgTAAGACTTGATCTGTGATAATTCTTTTGTTGCAGAAGACTATGTGATTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCg                                                                                                                                                <  8:262529/145‑1 (MQ=255)
    gggTATGACTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCgg                                                                                                                                               <  7:28547/148‑1 (MQ=255)
      gTATGACTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAATCTGAATTAAGAGAGCCAGCCATAAGGCg                                                                                                                                                >  7:262529/1‑145 (MQ=255)
       tATGACTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCTTATGGCgggctc                                                                                                                                           <  1:53286/149‑4 (MQ=255)
            cTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCtatta                                                                                                                                       <  7:137030/148‑1 (MQ=255)
               gATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAg                                                                                                                                   >  7:187443/1‑149 (MQ=255)
                     tgATAATTCTGTTGTTGCAGAATACAATGTGTTTTTGACAGAGATGGGGCGGGAATATGTATTTCATACACAAAAACCGCATATTTGGGAAGCTAGACCAAACTGAATTTAGAGAGCCCGCCAAAAGGCGGGTCTATTTTTAGTAACGt                                                                                                                             >  7:55410/1‑149 (MQ=255)
                     tgATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGt                                                                                                                             >  1:18775/1‑149 (MQ=255)
                          aTTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTAtt                                                                                                                         <  5:296526/148‑1 (MQ=255)
                           ttCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGTAAGCTAGACGAAACTGAATTAAGAGATCCCGACATAAGGCGGGTCTATTATGAGTCAGGTTAtatc                                                                                                                       >  8:327361/1‑146 (MQ=255)
                             cTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGa                                                                                                                      <  8:60267/148‑1 (MQ=255)
                                ttgttgCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTAGGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGata                                                                                                                    <  8:89121/147‑1 (MQ=255)
                                      cAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCt                                                                                                            >  5:358844/1‑149 (MQ=255)
                                        gaagaCTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTgg                                                                                                          >  4:242580/1‑149 (MQ=255)
                                           gaCTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCtt                                                                                                       >  5:304997/1‑149 (MQ=255)
                                             cTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCttt                                                                                                      <  8:187443/148‑1 (MQ=255)
                                              tATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTcc                                                                                                    <  1:95634/149‑1 (MQ=255)
                                                tgtgTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTgt                                                                                                          <  5:129441/141‑2 (MQ=255)
                                                tgtgTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTgg                                                                                                          >  6:129441/1‑141 (MQ=255)
                                                  tgtTTTTGACAGAGATGGTGCGGTAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACAt                                                                                                <  8:341193/149‑1 (MQ=255)
                                                    tttttGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATg                                                                                               <  2:18775/148‑1 (MQ=255)
                                                      tttGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGc                                                                                             >  3:180326/1‑148 (MQ=255)
                                                         gACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCAGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAg                                                                                          >  8:33042/1‑148 (MQ=255)
                                                          aCAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCa                                                                                        >  1:29440/1‑149 (MQ=255)
                                                            agagATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCAtg                                                                                      >  4:153377/1‑149 (MQ=255)
                                                             gagaTGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCAtgt                                                                                     <  5:156550/149‑1 (MQ=255)
                                                                 tGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGATAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAgg                                                                                 <  8:55410/149‑1 (MQ=255)
                                                                    ggCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCgg                                                                              <  2:103634/149‑1 (MQ=255)
                                                                     gctgtaGTATGTCTTGCATCCACAGAAGCCTCTTCTATTGTAAGCTAGACCAAACTGAATTAAGATAGCCCGCGATAAGGCGTGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGt                                                                             <  7:65794/144‑1 (MQ=255)
                                                                     gaggtaGTATCTCTTGAATCCACAGAAGACGCTTGTATGGGAAGCTAGTCAAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCgg                                                                              <  4:52846/143‑1 (MQ=255)
                                                                     gCGGTAGTATGTCTTGCATCCACAGAAGCCGCTTCTATTGGAAGCTAGACCAAACTGAATTAAGATCGCCCGCCATAAGGCGGGTCTATTATGAGTCTCGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGt                                                                             <  7:327361/149‑1 (MQ=255)
                                                                        ggAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGa                                                                          <  4:180326/149‑1 (MQ=255)
                                                                         gAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCAGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGAc                                                                         >  5:10785/1‑149 (MQ=255)
                                                                              tGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTcc                                                                    >  4:113647/1‑149 (MQ=255)
                                                                              tGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTcc                                                                    <  3:233134/149‑1 (MQ=255)
                                                                               gTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTcc                                                                    >  4:169569/1‑148 (MQ=255)
                                                                                tCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTccc                                                                   <  2:29440/148‑1 (MQ=255)
                                                                                 cTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTACCGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTccaca                                                                 >  6:64399/1‑146 (MQ=255)
                                                                                  ttGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAg                                                                >  3:274665/1‑149 (MQ=255)
                                                                                    gCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGa                                                               <  3:358313/148‑1 (MQ=255)
                                                                                    gCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGa                                                               <  5:153848/148‑1 (MQ=255)
                                                                                    gCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGAt                                                              >  6:85317/1‑149 (MQ=255)
                                                                                    gCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGAt                                                              >  6:199702/1‑149 (MQ=255)
                                                                                      aTCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTc                                                            <  3:242580/149‑1 (MQ=255)
                                                                                        ccACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAg                                                          >  3:270513/1‑149 (MQ=255)
                                                                                          acaGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGtt                                                        >  4:286526/1‑149 (MQ=255)
                                                                                             gAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTg                                                     <  7:33042/149‑1 (MQ=255)
                                                                                              aaGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTgg                                                    <  2:47771/149‑1 (MQ=255)
                                                                                                  cGCTTCTATGGGAAGCTAGACCAAACTGAATTAATAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCaaa                                                <  6:10785/149‑1 (MQ=255)
                                                                                                      tCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGt                                                                                                              <  1:59460/83‑1 (MQ=255)
                                                                                                      tCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGt                                                                                                              >  2:59460/1‑83 (MQ=255)
                                                                                                        tATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGAtt                                                             <  6:152421/130‑1 (MQ=255)
                                                                                                        tATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGAtt                                                             >  5:152421/1‑130 (MQ=255)
                                                                                                            ggAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTc                                       <  5:85317/148‑1 (MQ=255)
                                                                                                                 cTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTc                                 <  4:274665/149‑1 (MQ=255)
                                                                                                                   aGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGt                               >  5:365525/1‑149 (MQ=255)
                                                                                                                    gACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGtt                              <  3:169569/149‑1 (MQ=255)
                                                                                                                        aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATgg                                                                                              <  4:327504/81‑1 (MQ=255)
                                                                                                                        aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATgg                                                                                              >  3:327504/1‑81 (MQ=255)
                                                                                                                        aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTg                                                     >  4:179974/1‑122 (MQ=255)
                                                                                                                        aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTg                                                     <  3:179974/122‑1 (MQ=255)
                                                                                                                                  aaGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                 <  6:358844/148‑1 (MQ=255)
                                                                                                                                  aaGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                 <  6:304997/148‑1 (MQ=255)
                                                                                                                                  aaGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                 <  3:153377/148‑1 (MQ=255)
                                                                                                                                   agagagCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTggg                >  7:170986/1‑148 (MQ=255)
                                                                                                                                      agacccGCCTTATGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTcc                                                                    <  3:113232/90‑1 (MQ=255)
                                                                                                                                      agacccGCCTTATGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTcc                                                                    >  4:113232/4‑93 (MQ=255)
                                                                                                                                        acccGCCTTATGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGtt                              >  7:312620/2‑129 (MQ=255)
                                                                                                                                        acccGCCTTATGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGtt                              <  8:312620/128‑1 (MQ=255)
                                                                                                                                            gCCTTATGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                 >  5:141227/1‑138 (MQ=255)
                                                                                                                                            gCCTTATGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                 <  6:141227/138‑1 (MQ=255)
                                                                                                                                            gCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAg      <  3:286526/149‑1 (MQ=255)
                                                                                                                                              cATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGt    >  1:41113/1‑149 (MQ=255)
                                                                                                                                                tAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCt  <  5:199702/149‑1 (MQ=255)
                                                                                                                                                 aaGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCtt                   <  2:297875/131‑1 (MQ=255)
                                                                                                                                                 aaGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCtt                   >  1:297875/1‑131 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GATGGGGTATGACTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTCT  >  NC_002947/1281048‑1281340

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: