Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,902,788 T→A 6.2% K102M (AAG→ATG) ‡ PP_5175 ← HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,902,7880TA 6.2% 256.7 / 6.9 83K102M (AAG→ATG) ‡PP_5175HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base T (37/40);  new base A (3/2);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

CTTGAGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACG  >  NC_002947/5902640‑5902937
                                                                                                                                                     |                                                                                                                                                        
cTTGAGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCtt                                                                                                                                                        >  2:1684063/1‑150 (MQ=255)
  tGAGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTgg                                                                                                                                                     <  1:504793/150‑1 (MQ=255)
    aGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTAGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGcc                                                                                                                                                   <  1:1209968/150‑1 (MQ=255)
      tGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGTCGa                                                                                                                                                 >  2:1212123/1‑150 (MQ=255)
        cgcTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGTCGAGt                                                                                                                                               <  1:1212123/150‑1 (MQ=255)
         gcTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTc                                                                                                                                              <  2:1722324/150‑1 (MQ=255)
         gcTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTc                                                                                                                                              <  1:200760/150‑1 (MQ=255)
           tctGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGa                                                                                                                                            <  2:682443/150‑1 (MQ=255)
              gATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACAt                                                                                                                                        >  1:181477/1‑150 (MQ=255)
                ttGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCTCTTGTGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTg                                                                                                                                      <  2:1423342/150‑1 (MQ=255)
                ttGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACTCTCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTg                                                                                                                                      <  1:1807179/150‑1 (MQ=255)
                    tgcGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGGCGACATTGTTCt                                                                                                                                  >  1:1246443/1‑150 (MQ=255)
                       gcAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCa                                                                                                                               >  1:1531827/1‑150 (MQ=255)
                       gcAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCa                                                                                                                               <  2:181477/150‑1 (MQ=255)
                       gcAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCGCGGCCTTGGCCGAGTCGACATTGTTCTGCa                                                                                                                               <  2:168048/150‑1 (MQ=255)
                        cAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAg                                                                                                                              <  2:1564819/150‑1 (MQ=255)
                          gcgcGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGt                                                                                                                            <  1:123124/150‑1 (MQ=255)
                                  tCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCggg                                                                                                                    >  1:463521/1‑150 (MQ=255)
                                              tACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCAAGGCCGAggtatt                                                                                                                                           >  2:1097989/1‑111 (MQ=255)
                                              tACTAGCTCTGGCTGTTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCAAGGCCGAggtatt                                                                                                                                           <  1:1097989/116‑6 (MQ=255)
                                                cTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAg                                                                                                      <  1:190920/150‑1 (MQ=255)
                                                 tCGCTCTGGCTGGTGTAGCCCTTCTTCGATAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAATCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGc                                                                                                     >  1:1718402/1‑150 (MQ=255)
                                                   gCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCAc                                                                                                   <  1:1670978/150‑1 (MQ=255)
                                                     tctGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCt                                                                                                 >  1:1785279/1‑150 (MQ=255)
                                                       tggctggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTgg                                                                                               <  2:1246443/150‑1 (MQ=255)
                                                         gctggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAg                                                                                                                              <  2:312172/117‑1 (MQ=255)
                                                         gctggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAg                                                                                                                              >  1:312172/1‑117 (MQ=255)
                                                         gctggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGtt                                                                                             >  1:324315/1‑150 (MQ=255)
                                                         gctggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGtt                                                                                             >  1:324325/1‑150 (MQ=255)
                                                         gctggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGtt                                                                                             >  1:105888/1‑150 (MQ=255)
                                                           tggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGc                                                                                                         <  1:1017007/136‑1 (MQ=255)
                                                           tggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGc                                                                                                         >  2:1017007/1‑136 (MQ=255)
                                                           tggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTgg                                                                                           <  1:57094/150‑1 (MQ=255)
                                                            ggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGc                                                                                          <  2:1531827/150‑1 (MQ=255)
                                                            ggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGc                                                                                          <  2:604218/150‑1 (MQ=255)
                                                              tgtAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCt                                                                                        >  1:1118415/1‑150 (MQ=255)
                                                                   cccTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTc                                                                                                                       <  1:1899336/114‑1 (MQ=255)
                                                                   cccTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTc                                                                                                                       >  2:1899336/1‑114 (MQ=255)
                                                                     cTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTg                                                                                 >  1:662308/1‑150 (MQ=255)
                                                                      ttGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGa                                                                                >  2:787440/1‑150 (MQ=255)
                                                                          gcagcaGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGc                                                                            >  1:1445690/1‑150 (MQ=255)
                                                                           cagcagTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCAAGGCCGAggta                                                                                                                                             >  1:1141796/1‑82 (MQ=255)
                                                                           cagcagTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCGCGGCCAAGGCCGAggta                                                                                                                                             <  2:1141796/85‑4 (MQ=255)
                                                                            agcagTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCCCGCCGAATACCTCGCCCCTGGTCGGGTCGGCATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGccc                                                                          >  1:1205998/1‑150 (MQ=255)
                                                                                gTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGAt                                                                               <  2:939411/141‑1 (MQ=255)
                                                                                gTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGAt                                                                               >  1:939411/1‑141 (MQ=255)
                                                                                gTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCGGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCATGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGt                                                                      <  1:318095/150‑1 (MQ=255)
                                                                                                  aGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGa                                                    >  1:1248017/1‑150 (MQ=255)
                                                                                                    gCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATc                                                  <  1:982841/150‑1 (MQ=255)
                                                                                                    gCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATc                                                  >  1:1079288/1‑150 (MQ=255)
                                                                                                    gCCGCACTGCCCTGGGTGCGCCGCGCCTGCGCGGCGAATACATCGGGCCTATATGCCGA‑‑CGGAGTGTGTTGGATGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATc                                                  <  1:1684063/150‑1 (MQ=255)
                                                                                                             gTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGa                                                                                     <  1:1891100/106‑1 (MQ=255)
                                                                                                             gTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGa                                                                                     >  2:1891100/1‑106 (MQ=255)
                                                                                                              tCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATg                                                     <  2:1111872/137‑1 (MQ=255)
                                                                                                              tCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATg                                                     >  1:1111872/1‑137 (MQ=255)
                                                                                                               cTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCTCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGAc                                       <  2:1118415/150‑1 (MQ=255)
                                                                                                                  ggTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGAc                                                                                    <  2:1525779/102‑1 (MQ=255)
                                                                                                                  ggTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGAc                                                                                    >  1:1525779/1‑102 (MQ=255)
                                                                                                                   gTGACCCGCGCCTGCGCGGCGAATACCTCGGGCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCgg                                   >  2:1194759/1‑150 (MQ=255)
                                                                                                                    tGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGcc                                                                                                                                                   >  1:1110033/1‑38 (MQ=255)
                                                                                                                    tGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGcc                                                                                                                                                   <  2:1110033/38‑1 (MQ=255)
                                                                                                                    tGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGa                                  <  2:662308/150‑1 (MQ=255)
                                                                                                                        ccGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGa                              >  2:800047/1‑150 (MQ=255)
                                                                                                                        ccGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGa                              >  1:262670/1‑150 (MQ=255)
                                                                                                                         cgcgCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCttt                                                                                                                        <  1:1397184/59‑1 (MQ=255)
                                                                                                                         cgcgCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCttt                                                                                                                        >  2:1397184/1‑59 (MQ=255)
                                                                                                                                  gcggcgAATACCTCGGCCTTGGCCGAGCGAGATTTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTc                    >  1:542678/1‑150 (MQ=255)
                                                                                                                                   cggcgAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCt                   <  2:324315/150‑1 (MQ=255)
                                                                                                                                   cggcgAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCt                   <  2:324325/150‑1 (MQ=255)
                                                                                                                                     gcgAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGc                 >  2:1705481/1‑150 (MQ=255)
                                                                                                                                       gAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGAcgc                                     <  1:1495868/128‑1 (MQ=255)
                                                                                                                                       gAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGAcgc                                     >  2:1495868/1‑128 (MQ=255)
                                                                                                                                           gacTCGGCCAAGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCcgc           >  1:1875129/3‑150 (MQ=255)
                                                                                                                                             cTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGcc                                                                           <  2:292824/84‑1 (MQ=255)
                                                                                                                                             cTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGcc                                                                           >  1:292824/1‑84 (MQ=255)
                                                                                                                                             cTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCgcgc         <  1:787440/150‑1 (MQ=255)
                                                                                                                                             cTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCgcgc         >  2:454875/1‑150 (MQ=255)
                                                                                                                                             cTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCgcgc         >  1:1381777/1‑150 (MQ=255)
                                                                                                                                              tCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGcc        <  1:1194759/150‑1 (MQ=255)
                                                                                                                                               cGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCt       <  2:463521/150‑1 (MQ=255)
                                                                                                                                               cGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCt       >  1:1719500/1‑150 (MQ=255)
                                                                                                                                                    cTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGaaa                                                       >  2:659852/1‑98 (MQ=255)
                                                                                                                                                    cTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGaaa                                                       <  1:659852/98‑1 (MQ=255)
                                                                                                                                                    cTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCAc   >  2:744815/1‑150 (MQ=255)
                                                                                                                                                     ttGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACg  >  2:568533/1‑150 (MQ=255)
                                                                                                                                                     ttGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCAGAAGAGGACTGACAACTCGGTCTAGACCCGCGccgacccg  >  2:112476/1‑144 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
CTTGAGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACG  >  NC_002947/5902640‑5902937

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