Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,281,191 A→T 8.9% intergenic (+17/+26) 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,1910AT 8.9% 211.2 / 8.3 75intergenic (+17/+26)PP_5742/PP_1118membrane protein/recombinase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base T (3/4);  total (34/41)
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

ATTGATGGGGTATGACTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGT  >  NC_002947/1281045‑1281338
                                                                                                                                                  |                                                                                                                                                   
aTTGATGGGGTATGACTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCAt                                                                                                                                                    >  1:103634/1‑148 (MQ=255)
   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                                                             <  5:153848/148‑1 (MQ=255)
                                                                                       gCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGa                                                             <  3:358313/148‑1 (MQ=255)
                                                                                       gCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGAt                                                            >  6:199702/1‑149 (MQ=255)
                                                                                       gCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGAt                                                            >  6:85317/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                                                                                            >  3:327504/1‑81 (MQ=255)
                                                                                                                           aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATgg                                                                                            <  4:327504/81‑1 (MQ=255)
                                                                                                                           aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTg                                                   <  3:179974/122‑1 (MQ=255)
                                                                                                                           aaaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTg                                                   >  4:179974/1‑122 (MQ=255)
                                                                                                                                     aaGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg               <  3:153377/148‑1 (MQ=255)
                                                                                                                                     aaGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg               <  6:304997/148‑1 (MQ=255)
                                                                                                                                     aaGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg               <  6:358844/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               <  6:141227/138‑1 (MQ=255)
                                                                                                                                               gCCTTATGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg               >  5:141227/1‑138 (MQ=255)
                                                                                                                                               gCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAg    <  3:286526/149‑1 (MQ=255)
                                                                                                                                                 cATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGt  >  1:41113/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATTGATGGGGTATGACTTGATCTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGT  >  NC_002947/1281045‑1281338

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: