Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,266,266 C→A 8.0% intergenic (‑50/+34) zinU ← / ← PP_4640 P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,266,2660CA 8.0% 582.5 / 30.6 174intergenic (‑50/+34)zinU/PP_4640P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (80/80);  new base A (8/6);  total (88/86)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 176 total aligned reads displayed.

GGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACA  >  NC_002947/5266118‑5266414
                                                                                                                                                     |                                                                                                                                                    
ggCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGACGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTAcc                                                                                                                                                     >  2:1017221/1‑150 (MQ=255)
  cGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCa                                                                                                                                                   <  1:3605600/150‑1 (MQ=255)
        tCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTcc                                                                                                                                             <  1:945818/150‑1 (MQ=255)
          aGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACCGGCTCCca                                                                                                                                           <  1:218601/150‑1 (MQ=255)
          aGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCca                                                                                                                                           >  2:18624/1‑150 (MQ=255)
           ggCCATGCTCGGCCTTGAGCATGTACTTGAGCGAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACCGGCTCCcac                                                                                                                                          <  1:352654/150‑1 (MQ=255)
            gCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCcaca                                                                                                                                         >  2:78931/1‑150 (MQ=255)
            gCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCcaca                                                                                                                                         >  1:2384616/1‑150 (MQ=255)
             ccATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAg                                                                                                                                        <  1:3194890/150‑1 (MQ=255)
             ccATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAg                                                                                                                                        >  1:37569/1‑150 (MQ=255)
              cATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCca                                                                                                                                           >  1:3205052/1‑146 (MQ=255)
              cATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCca                                                                                                                                           <  2:3205052/146‑1 (MQ=255)
                 gCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGttt                                                                                                                                    >  1:3595997/1‑150 (MQ=255)
                  cTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACCGGCTCCCACAgtcatt                                                                                                                                   <  1:2736379/150‑6 (MQ=39)
                  cTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTa                                                                                                                                   <  1:132146/150‑1 (MQ=255)
                       ccTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGa                                                                                                                              <  1:2799430/150‑1 (MQ=255)
                               aTGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAgg                                                                                                                      <  2:670295/150‑1 (MQ=255)
                               aTGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAgg                                                                                                                      <  2:1290403/150‑1 (MQ=255)
                                   aCTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCa                                                                                                                  >  2:1306113/1‑150 (MQ=255)
                                   aCTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCa                                                                                                                  >  2:4863/1‑150 (MQ=255)
                                   aCTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCa                                                                                                                  >  1:3958056/1‑150 (MQ=255)
                                    cTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCaa                                                                                                                 >  1:1250963/1‑150 (MQ=255)
                                     ttGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAAc                                                                                                                >  1:1792538/1‑150 (MQ=255)
                                     ttGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAAc                                                                                                                >  2:2116443/1‑150 (MQ=255)
                                       gAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAg                                                                                                              <  1:333797/150‑1 (MQ=255)
                                        aGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGc                                                                                                             >  1:4282012/1‑150 (MQ=255)
                                          cAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCAc                                                                                                           >  1:84713/1‑150 (MQ=255)
                                            aGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACtt                                                                                                         <  1:56706/150‑1 (MQ=255)
                                             ggtggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTg                                                                                                        <  1:1615973/150‑1 (MQ=255)
                                             ggtggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTg                                                                                                        >  1:3866634/1‑150 (MQ=255)
                                             ggtggtCTTGCCGGCACCGAGGAAGCCTGAGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCACAACAGCACTTg                                                                                                        <  2:4501793/150‑1 (MQ=255)
                                              gtggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                       <  1:1094681/150‑1 (MQ=255)
                                              gtggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                       <  1:719555/150‑1 (MQ=255)
                                              gtggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                       <  1:1689302/150‑1 (MQ=255)
                                                ggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGc                                                                                                     <  1:2520691/150‑1 (MQ=255)
                                                 gtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGcc                                                                                                    <  2:1444387/150‑1 (MQ=255)
                                                  tCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGccc                                                                                                   >  2:1779838/1‑150 (MQ=255)
                                                    ttGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCgg                                                                                                 >  1:2516784/1‑150 (MQ=255)
                                                    ttGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCgg                                                                                                 >  1:2851898/1‑150 (MQ=255)
                                                       ccGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACt                                                                                              >  1:4530/1‑150 (MQ=255)
                                                        cGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACtt                                                                                             >  1:1543588/1‑150 (MQ=255)
                                                        cGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACtt                                                                                             <  2:3624403/150‑1 (MQ=255)
                                                         ggCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTg                                                                                            >  1:2919597/1‑150 (MQ=255)
                                                            aCCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCa                                                                                         >  2:2087849/1‑150 (MQ=255)
                                                              cGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCAcc                                                                                       <  1:854208/150‑1 (MQ=255)
                                                                       cTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGc                                                                              <  1:2612236/150‑1 (MQ=255)
                                                                       cTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGc                                                                              >  1:3514068/1‑150 (MQ=255)
                                                                        tgtgAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCt                                                                             >  2:2211598/1‑150 (MQ=255)
                                                                            aGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCt                                                                         >  1:3761044/1‑150 (MQ=255)
                                                                             ggACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTg                                                                        >  2:1532642/1‑150 (MQ=255)
                                                                             ggACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTg                                                                        >  2:2402749/1‑150 (MQ=255)
                                                                                   ttACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTc                                                                  <  1:3488921/150‑1 (MQ=255)
                                                                                    tACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTcg                                                                 <  1:2293689/150‑1 (MQ=255)
                                                                                     aCGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTcgc                                                                >  2:2445115/1‑150 (MQ=255)
                                                                                      cGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTcgcg                                                               <  1:833438/150‑1 (MQ=255)
                                                                                         gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                            <  2:3957694/150‑1 (MQ=255)
                                                                                         gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                            >  2:1719268/1‑150 (MQ=255)
                                                                                         gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                            <  1:314373/150‑1 (MQ=255)
                                                                                         gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                            <  2:4221001/150‑1 (MQ=255)
                                                                                         gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                            <  2:4331468/150‑1 (MQ=255)
                                                                                         gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                            <  2:4343544/150‑1 (MQ=255)
                                                                                         gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                            <  2:490205/150‑1 (MQ=255)
                                                                                          aaTGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaag                                                           <  1:3637697/150‑1 (MQ=255)
                                                                                            tGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaagaa                                                         <  1:191988/150‑1 (MQ=255)
                                                                                             gggCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAAt                                                        <  1:4209404/150‑1 (MQ=255)
                                                                                             gggCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAAt                                                        <  2:2262473/150‑1 (MQ=255)
                                                                                               gCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGc                                                      >  2:342513/1‑150 (MQ=255)
                                                                                                cATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGcc                                                     >  1:1697468/1‑150 (MQ=255)
                                                                                                 aTCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCt                                                    <  1:788528/150‑1 (MQ=255)
                                                                                                  tCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTc                                                   <  1:2110325/150‑1 (MQ=255)
                                                                                                  tCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTc                                                   <  2:165726/150‑1 (MQ=255)
                                                                                                      cACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAg                                               <  2:3373906/150‑1 (MQ=255)
                                                                                                      cACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAg                                               >  2:2414370/1‑150 (MQ=255)
                                                                                                         gCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCtc                                            <  1:1035017/150‑1 (MQ=255)
                                                                                                         gCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCtc                                            >  2:2135789/1‑150 (MQ=255)
                                                                                                                 cgccTGAATCTGACAATGACTGGGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGCACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCGTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCgg                                    <  2:3897414/148‑1 (MQ=255)
                                                                                                                  tCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGc                                   <  1:4863/150‑1 (MQ=255)
                                                                                                                    cTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCtt                                 >  2:2676395/1‑150 (MQ=255)
                                                                                                                     tGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTg                                >  1:3074810/1‑150 (MQ=255)
                                                                                                                     tGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTg                                >  1:3360702/1‑150 (MQ=255)
                                                                                                                     tGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTg                                >  2:1171579/1‑150 (MQ=255)
                                                                                                                     tGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTg                                >  2:3107170/1‑150 (MQ=255)
                                                                                                                        aTCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTcc                             <  2:1742528/150‑1 (MQ=255)
                                                                                                                          cTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCgg                           <  2:2172936/150‑1 (MQ=255)
                                                                                                                           tGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGa                          <  2:2489690/150‑1 (MQ=255)
                                                                                                                              aaaTGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGc                       <  1:650059/150‑1 (MQ=255)
                                                                                                                                aaaccTGTGGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGt                     >  2:49892/5‑150 (MQ=255)
                                                                                                                                  accTGTGGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGttt                   >  2:2262640/3‑150 (MQ=255)
                                                                                                                                     tgtgGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGGCGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCt                >  2:2851131/1‑150 (MQ=255)
                                                                                                                                     tgtgGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCt                >  2:2636500/1‑150 (MQ=255)
                                                                                                                                     tgtgGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCt                >  1:1449085/1‑150 (MQ=255)
                                                                                                                                     tgtgGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCt                <  1:4126077/150‑1 (MQ=255)
                                                                                                                                       tgGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGc              >  1:631315/1‑150 (MQ=255)
                                                                                                                                         ggAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt            >  2:4297933/1‑150 (MQ=255)
                                                                                                                                         ggAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt            <  1:2984405/150‑1 (MQ=255)
                                                                                                                                         ggAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt            <  1:2670020/150‑1 (MQ=255)
                                                                                                                                           aGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTGTCTGCAtgt          >  2:3351797/1‑150 (MQ=255)
                                                                                                                                              cGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGtt       <  2:2243166/150‑1 (MQ=255)
                                                                                                                                                ggTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTc      <  1:1700735/150‑1 (MQ=255)
                                                                                                                                                    aCCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACa  <  2:608316/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACA  >  NC_002947/5266118‑5266414

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