Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,382,332 A→T 6.2% intergenic (+12/+15) dinB → / ← PP_1204 DNA polymerase IV/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,382,3320AT 6.2% 251.7 / 5.9 80intergenic (+12/+15)dinB/PP_1204DNA polymerase IV/lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base T (2/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAGAGTTATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTC  >  NC_002947/1382187‑1382478
                                                                                                                                                 |                                                                                                                                                  
tGGAGAGTTATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGGTAGTGTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATCCa                                                                                                                                                >  1:1321254/1‑150 (MQ=255)
  gagaGTTATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGc                                                                                                                                              <  2:1280729/150‑1 (MQ=255)
   agagTTATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGcc                                                                                                                                             <  1:1941453/150‑1 (MQ=255)
    gagTTATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCATTGCAGccc                                                                                                                                            <  1:1742789/150‑1 (MQ=255)
       ttATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCtt                                                                                                                                         <  2:1957056/150‑1 (MQ=255)
       ttATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCtt                                                                                                                                         >  1:91708/1‑150 (MQ=255)
        tATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCttt                                                                                                                                        <  1:238260/150‑2 (MQ=255)
         aTCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCt                                                                                                                                       >  2:647455/1‑150 (MQ=255)
               aaTTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCATTGCAGcccc                                                                                                                                           >  1:306041/1‑140 (MQ=255)
               aaTTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCATTGCAGcccc                                                                                                                                           <  2:306041/140‑1 (MQ=255)
                 ttGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGc                                                                                                                               <  2:1686148/150‑1 (MQ=255)
                 ttGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTCCGGCGGCCGGCCTGGCTTGTCGCGCATGAACAGTTGGAGCTGTTTACGACCAAATGAAAGGTGCTGCAATGCAGCCCCTTCTATTAATGc                                                                                                                               <  1:1695414/150‑1 (MQ=255)
                   gCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCaa                                                                                                                             >  2:1695414/1‑149 (MQ=255)
                   gCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCAc                                                                                                                             <  1:647455/150‑1 (MQ=255)
                     tGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGc                                                                                                                           <  2:770872/150‑1 (MQ=255)
                          cAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGa                                                                                                                      <  2:1610912/150‑1 (MQ=255)
                             gCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATcc                                                                                                                   >  1:1631376/1‑150 (MQ=255)
                                ttCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGcc                                                                                                                >  1:342444/1‑150 (MQ=255)
                                ttCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGcc                                                                                                                >  1:1121731/1‑150 (MQ=255)
                                  cGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGccac                                                                                                              >  2:1594018/1‑150 (MQ=255)
                                     ccGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCaa                                                                                                           >  1:608427/1‑150 (MQ=255)
                                       gTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAAt                                                                                                                                                   <  2:1288875/108‑1 (MQ=255)
                                       gTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAAt                                                                                                                                                   >  1:1288875/1‑108 (MQ=255)
                                           cggcAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATAAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGcc                                                                                                     >  1:2308049/1‑150 (MQ=255)
                                            ggcAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTc                                                                                                                                        <  2:1239721/114‑1 (MQ=255)
                                            ggcAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTc                                                                                                                                        >  1:1239721/1‑114 (MQ=255)
                                                 gCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTc                                                                                                                                        <  1:1535844/109‑1 (MQ=255)
                                                 gCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTc                                                                                                                                        >  2:1535844/1‑109 (MQ=255)
                                                     gTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTg                                                                                           <  2:1507338/150‑1 (MQ=255)
                                                     gTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTg                                                                                           >  2:2143896/1‑150 (MQ=255)
                                                      tGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTgg                                                                                          <  2:820479/150‑1 (MQ=255)
                                                           ttGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGt                                                                                                  <  2:2000481/137‑1 (MQ=255)
                                                           ttGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGt                                                                                                  >  1:2000481/1‑137 (MQ=255)
                                                                 ggggtaggggtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTAt                                                                                                                                     <  2:826555/96‑1 (MQ=255)
                                                                 ggggtaggggtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTAt                                                                                                                                     >  1:826555/1‑96 (MQ=255)
                                                                  gggtaggggtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCTCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCa                                                                              <  1:1615085/150‑1 (MQ=255)
                                                                      aggggtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGaa                                                                          <  2:2144881/150‑1 (MQ=255)
                                                                        gggtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACt                                                                        >  1:333921/1‑150 (MQ=255)
                                                                          gtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTc                                                                                                                                        >  1:1930013/1‑84 (MQ=255)
                                                                          gtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTc                                                                                                                                        <  2:1930013/84‑1 (MQ=255)
                                                                          gtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTcc                                                                      >  2:621342/1‑150 (MQ=255)
                                                                          gtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTcc                                                                      >  1:491562/1‑150 (MQ=255)
                                                                           tGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCct                                                                     >  1:1516563/1‑150 (MQ=255)
                                                                             aGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGAGTGCAGGAACTCctac                                                                   >  2:1623793/1‑148 (MQ=255)
                                                                                 tGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGcca                                                                                                               <  2:1358982/102‑1 (MQ=255)
                                                                                 tGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGcca                                                                                                               >  1:1358982/1‑102 (MQ=255)
                                                                                    gTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGTATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGCAACTCCTGCTGCAGCt                                                            >  2:826141/1‑150 (MQ=255)
                                                                                         ttGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGa                                                       <  2:91708/150‑1 (MQ=255)
                                                                                         ttGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACAGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGa                                                       >  2:2176690/1‑150 (MQ=255)
                                                                                              gggCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCctgc                                                                   >  2:2180725/1‑133 (MQ=255)
                                                                                              gggCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCctgc                                                                   <  1:2180725/133‑1 (MQ=255)
                                                                                                gcgcgcATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGc                                                <  1:1510657/150‑1 (MQ=255)
                                                                                                  gcgcATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGc                                                             >  2:298018/1‑135 (MQ=255)
                                                                                                  gcgcATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGc                                                             <  1:298018/135‑1 (MQ=255)
                                                                                                      aTGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCa                                          <  2:1516563/150‑1 (MQ=255)
                                                                                                      aTGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCa                                          >  1:138281/1‑150 (MQ=255)
                                                                                                              ttGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCa                                                                              <  2:1753275/106‑1 (MQ=255)
                                                                                                              ttGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCa                                                                              >  1:1753275/1‑106 (MQ=255)
                                                                                                              ttGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATc                                  >  1:1995566/1‑150 (MQ=255)
                                                                                                                 gAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCctgctg                                                                 <  2:1838984/116‑1 (MQ=255)
                                                                                                                 gAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCctgctg                                                                 >  1:1838984/1‑116 (MQ=255)
                                                                                                                   gCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGct                             <  2:333921/150‑1 (MQ=255)
                                                                                                                      gTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTg                          >  2:1328191/1‑150 (MQ=255)
                                                                                                                        ttCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGtt                        <  1:2143896/150‑1 (MQ=255)
                                                                                                                           cGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCa                     <  2:342444/150‑1 (MQ=255)
                                                                                                                             cccAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGt                                                                                                  >  2:609356/1‑71 (MQ=255)
                                                                                                                             cccAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGt                                                                                                  <  1:609356/71‑1 (MQ=255)
                                                                                                                               ttaaTAAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTc                 >  1:1736983/3‑150 (MQ=255)
                                                                                                                               cAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTc                 <  2:1631376/150‑1 (MQ=255)
                                                                                                                                  aTGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGt                                            <  2:2247106/120‑1 (MQ=255)
                                                                                                                                  aTGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGt                                            >  1:2247106/1‑120 (MQ=255)
                                                                                                                                   tGAAAGGGGCTGCAATGCATCCCCTTCTATTAATGCACGCCGGGATTCGCCACCAACCGGCCATTGGCGTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGCTCAGGCGCTGTTCACTTGCACTg             >  2:331051/1‑150 (MQ=255)
                                                                                                                                     gaaGGGGCTGCATTGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGc           >  1:130473/2‑150 (MQ=255)
                                                                                                                                      aaGGGGCTGCATTGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCt          <  2:130473/150‑1 (MQ=255)
                                                                                                                                      aaGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCt          >  2:126732/1‑150 (MQ=255)
                                                                                                                                          ggCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCct      <  2:138281/150‑1 (MQ=255)
                                                                                                                                           gCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCctc     >  1:2090252/1‑150 (MQ=255)
                                                                                                                                             tGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCtt   >  1:1928830/1‑150 (MQ=255)
                                                                                                                                              gCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCATCTTc  >  1:267885/1‑150 (MQ=255)
                                                                                                                                              gCAATGAAGCCCCGTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTATTc  >  2:1644799/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGGAGAGTTATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTC  >  NC_002947/1382187‑1382478

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: