Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,698,803 T→C 7.1% intergenic (+37/‑130) hslU → / → PP_5002 protease HslVU ATPase subunit/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,698,8030TC 7.1% 194.9 / 4.3 70intergenic (+37/‑130)hslU/PP_5002protease HslVU ATPase subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base C (3/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

CCTGGCCAGCACCCATGACGAAGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATG  >  NC_002947/5698658‑5698935
                                                                                                                                                 |                                                                                                                                    
ccTGGCCAGCACCCATGACGAAGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCCCCCGCAc                                                                                                                                  >  1:153159/1‑150 (MQ=255)
   ggCCAGCACCCATGACGAAGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCCGCAcgga                                                                                                                               >  2:354157/1‑148 (MQ=255)
     ccAGCACCCATGACGAAGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCCGCAcggatc                                                                                                                             <  1:354157/150‑5 (MQ=255)
       aGCACCCATGACGAAGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGccac                                                                                                                           >  1:537598/1‑150 (MQ=255)
              aTGACGAAGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCTCCTCTGCGCGCGCCACCACAAcc                                                                                                                    >  2:1229852/1‑150 (MQ=255)
                gACGAAGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCt                                                                                                                  >  2:1144163/1‑150 (MQ=255)
                   gAAGCGCCGATCCAGATCGACGCGACCAGAATGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTgg                                                                                                               >  1:358903/1‑150 (MQ=255)
                     aGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGc                                                                                                             <  2:77724/150‑1 (MQ=255)
                     aGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGc                                                                                                             <  1:187118/150‑1 (MQ=255)
                             tCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTgtggtg                                                                                                     <  1:87431/150‑1 (MQ=255)
                              ccAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGa                                                                                                    >  2:1135547/1‑150 (MQ=255)
                                aGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATc                                                                                                  >  1:97751/1‑150 (MQ=255)
                                  aTCGACGCGGCGTATGTGAACAGCCACCTGGGTGCGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTTATCCg                                                                                                >  1:539288/1‑150 (MQ=255)
                                      aCGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCg                                                                                            <  1:1135547/150‑1 (MQ=255)
                                            ggcaggcgggccatgggcagaatctctcgggagctgatgaagcacagcttccagcttatgggggagggattcaagggcctcTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGgggcg                                                                                      <  2:114377‑M1/69‑1 (MQ=255)
                                             gTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGAt                                                                                                   >  2:536275/1‑136 (MQ=255)
                                             gTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGAt                                                                                                   <  1:536275/136‑1 (MQ=255)
                                             gTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGgggcgg                                                                                     <  1:812029/150‑1 (MQ=255)
                                             gTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGgggcgg                                                                                     <  1:812234/150‑1 (MQ=255)
                                             gTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGgggcgg                                                                                     <  1:1144163/150‑1 (MQ=255)
                                             gTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGgggcgg                                                                                     <  2:537598/150‑1 (MQ=255)
                                               aggcgggccatgggcagaatctctcgggagctgatgaagcacagcttccagcttatgggggagggattcaagggcctcTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGgggcgggc                                                                                   >  1:114377‑M1/79‑150 (MQ=255)
                                               aTGTGAACAGCCACCTGGGTGAGCCGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAACCTTCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGgggcgggc                                                                                   <  2:1138546/150‑1 (MQ=255)
                                                  tgAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTg                                                                                >  1:584869/1‑150 (MQ=255)
                                                     aCAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGccc                                                                             >  1:144741/1‑150 (MQ=255)
                                                               gggTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGc                                                                   >  2:1126763/1‑150 (MQ=255)
                                                                ggTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGcc                                                                  >  2:1242132/1‑150 (MQ=255)
                                                                      cTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGaa                                                            >  2:573078/1‑150 (MQ=255)
                                                                       tGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGgggcgggc                                                                                   >  2:581497/1‑126 (MQ=255)
                                                                       tGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGgggcgggc                                                                                   <  1:581497/126‑1 (MQ=255)
                                                                       tGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAt                                                           <  1:1126763/150‑1 (MQ=255)
                                                                       tGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAt                                                           <  2:388639/150‑1 (MQ=255)
                                                                           acaGAATGTAGACCTGTCGCGTGACATTCGGGAAGACCGGTTGCGACAGTTTCGCGGGCAAGCCCGCGCCTGCACGCGCCACCTCGACCCTTGGGCCTGTTGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAtccct                                                       <  2:539288/150‑1 (MQ=255)
                                                                           acaGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAtccct                                                       >  2:269087/1‑150 (MQ=255)
                                                                            gctgatgaagcacagcttccagcttatgggggagggattcaagggcctcTTCGCGGGCAAGCCCGCCCCCGCACGCGCCACCACAACCCTTGGGCCtgtg                                                                                                        <  1:1070547‑M1/51‑1 (MQ=255)
                                                                            gctgatgaagcacagcttccagcttatgggggagggattcaagggcctcTTCGCGGGCAAGCCCGCCCCCGCACGCGCCACCACAACCCTTGGGCCtgtg                                                                                                        >  2:1070547‑M1/50‑100 (MQ=255)
                                                                            caGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAtccctc                                                      >  2:1076664/1‑150 (MQ=255)
                                                                             aGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAtccctcc                                                     >  1:239042/1‑150 (MQ=255)
                                                                              gAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAtccctccc                                                    >  1:1284688/1‑150 (MQ=255)
                                                                               aaCGAAGACCGGTGGCGTTACATTCGGTAAGGCGGGTTTCGACGGTTTCGCGGTCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCGTGCCCGCGAAGAGGCCCTTGAACCCCTcccc                                                   <  1:1229852/150‑1 (MQ=255)
                                                                                     gcacagcttccagcttatgggggagggattcaagggcctcTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTgg                                                                                                               <  2:762741‑M1/44‑1 (MQ=255)
                                                                                     gcacagcttccagcttatgggggagggattcaagggcctcTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTgg                                                                                                               >  1:762741‑M1/41‑84 (MQ=255)
                                                                                        cTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTg                                          <  2:584869/150‑1 (MQ=255)
                                                                                         tGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTgg                                         >  1:225019/1‑150 (MQ=255)
                                                                                           tCGCGTTACATTCTGCAAGACCGGTTTCGACAGTTTGGCGGGCACGCTTCCGATTTCGCTGGCCACCCCAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGaa                                       <  2:358903/150‑1 (MQ=255)
                                                                                                   cATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCt                               <  1:573078/150‑1 (MQ=255)
                                                                                                    aTTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCggg                                                                                          <  1:1180619/90‑1 (MQ=255)
                                                                                                    aTTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCggg                                                                                          >  2:1180619/1‑90 (MQ=255)
                                                                                                       cTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGATGGAAGCTGTGCTtcat                           <  1:43893/150‑1 (MQ=255)
                                                                                                        tGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatc                          >  2:569919/1‑150 (MQ=255)
                                                                                                        actAAGACCGGTTTCGACAGTTTCGCGGTGAAGCCTGCCCCTGCCCGCGCCCCCCCACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTAAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatc                          <  2:144741/148‑1 (MQ=255)
                                                                                                         gTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatca                         <  1:309316/150‑1 (MQ=255)
                                                                                                         gTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatca                         <  2:225019/150‑1 (MQ=255)
                                                                                                         gTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatca                         <  2:1284688/150‑1 (MQ=255)
                                                                                                          tAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCCAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAg                        >  1:749663/1‑150 (MQ=255)
                                                                                                          tAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCAAGCCcgccccc                                                                       >  1:527966/1‑99 (MQ=39)
                                                                                                          tAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCAAGCCcgccccc                                                                       <  2:527966/103‑5 (MQ=39)
                                                                                                           aaGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGc                       <  2:97751/150‑1 (MQ=255)
                                                                                                              aCCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGc                                                                   >  2:692151/1‑103 (MQ=255)
                                                                                                              aCCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGc                                                                   <  1:692151/103‑1 (MQ=255)
                                                                                                                cGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGACTATTTTGTCCATAAGCTGGAAACTGTGCTTCATCAGCTCCCg                  >  1:708055/1‑150 (MQ=255)
                                                                                                                        aCAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCGCCACACCCGCTGGTCCTGTGTTGATCCGTGCGGGGGCGGGCTTTCCCGCGAAGAGACCCTTGAATCCCTCCCCCATAAGCTGGACGCTGTGCTTCATCAGCTCCCGAGAGATTc          >  2:447315/1‑150 (MQ=255)
                                                                                                                        aCAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCtgt                                  >  1:1134903/1‑126 (MQ=255)
                                                                                                                        aCAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCtgt                                  <  2:1134903/126‑1 (MQ=255)
                                                                                                                        aCAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTc          >  2:86729/1‑150 (MQ=255)
                                                                                                                                GCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATg  >  1:269087‑M1/1‑150 (MQ=255)
                                                                                                                                    GCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGc                       >  2:935947‑M1/1‑125 (MQ=255)
                                                                                                                                    GCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGc                       <  1:935947‑M1/125‑1 (MQ=255)
                                                                                                                                         CCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGc                                >  2:349768‑M1/1‑111 (MQ=255)
                                                                                                                                         CCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGc                                <  1:349768‑M1/111‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                    
CCTGGCCAGCACCCATGACGAAGCGCCGATCCAGATCGACGCGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATG  >  NC_002947/5698658‑5698935

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