Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,506 C→G 5.1% R194P (CGT→CCT)  ftsX ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,5060CG 5.1% 851.8 / 26.9 234R194P (CGT→CCT) ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (101/121);  new base G (9/3);  total (110/124)
Fisher's exact test for biased strand distribution p-value = 7.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

ACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGA  >  NC_000913/3601358‑3601655
                                                                                                                                                    |                                                                                                                                                     
acTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACg                                                                                                                                                      >  2:1111996/1‑150 (MQ=255)
acTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACg                                                                                                                                                      >  1:3207777/1‑150 (MQ=255)
        cGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCa                                                                                                                                              >  1:602378/1‑150 (MQ=255)
         gATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAg                                                                                                                                             >  1:2749921/1‑150 (MQ=255)
              ngAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTGCGTCAATGCCGTTAATCTGCGTGATACGAGCACGCAGCGTAt                                                                                                                                        >  1:3305155/2‑150 (MQ=255)
              cGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTAt                                                                                                                                        >  1:4245896/1‑150 (MQ=255)
               gAGGAACACGGCCGCCACCATCAACACGCCGCTCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTAtt                                                                                                                                       <  2:2902035/150‑1 (MQ=255)
                aGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATCGTATCACGCAgattaacg                                                                                                                                      <  2:3283942/150‑8 (MQ=255)
                aGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTc                                                                                                                                      <  2:1145517/150‑1 (MQ=255)
                 ggAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCa                                                                                                                                     <  2:3026932/150‑1 (MQ=255)
                    acacGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTg                                                                                                                                  >  2:102610/1‑150 (MQ=255)
                    acacGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTg                                                                                                                                  >  2:3263722/1‑150 (MQ=255)
                    acacGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTg                                                                                                                                  >  2:723460/1‑150 (MQ=255)
                     cacGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGa                                                                                                                                 <  2:297038/150‑1 (MQ=255)
                         gccgccACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGAAGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTcc                                                                                                                             <  2:1722970/150‑1 (MQ=255)
                          ccgccACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCg                                                                                                                            <  1:1844923/150‑1 (MQ=255)
                             ccaccaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTcc                                                                                                                         <  1:1222948/150‑1 (MQ=255)
                              caccaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTccc                                                                                                                        >  2:3173487/1‑150 (MQ=255)
                              caccaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTccc                                                                                                                        <  2:2912254/150‑1 (MQ=255)
                              caccaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTccc                                                                                                                        >  2:351762/1‑150 (MQ=255)
                              caccaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTccc                                                                                                                        >  2:1390760/1‑150 (MQ=255)
                                 caTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTg                                                                                                                     >  1:2158170/1‑150 (MQ=255)
                                  aTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGTGTATTCAGTGATTCCGTCCCCTgg                                                                                                                    >  1:2147628/1‑150 (MQ=255)
                                  aTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTgg                                                                                                                    >  1:3448190/1‑150 (MQ=255)
                                  aTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTgg                                                                                                                    >  1:780032/1‑150 (MQ=255)
                                   tCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGa                                                                                                                   >  1:3003127/1‑150 (MQ=255)
                                   tCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGa                                                                                                                   >  1:2174426/1‑150 (MQ=255)
                                        acGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCg                                                                                                              >  2:2109552/1‑150 (MQ=255)
                                             gATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGttt                                                                                                         >  2:3658387/1‑150 (MQ=255)
                                               tcatcGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCg                                                                                                       <  2:1099877/150‑1 (MQ=255)
                                                catcGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCgg                                                                                                      <  2:369212/150‑1 (MQ=255)
                                                catcGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCgg                                                                                                      <  1:877013/150‑1 (MQ=255)
                                                   cGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGAt                                                                                                   <  2:815174/150‑1 (MQ=255)
                                                    gCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATc                                                                                                  >  2:2317696/1‑150 (MQ=255)
                                                     ccGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATca                                                                                                 <  1:3249255/150‑1 (MQ=255)
                                                     ccGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATca                                                                                                 >  2:2811255/1‑150 (MQ=255)
                                                        aaaCGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATcacca                                                                                              <  2:3388344/150‑1 (MQ=255)
                                                         aaCGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATcaccac                                                                                             <  1:1138296/150‑1 (MQ=255)
                                                          aCGCGCCCGCCCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCaccacc                                                                                            <  2:638602/150‑1 (MQ=255)
                                                               gccGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCAc                                                                                       <  2:554051/149‑1 (MQ=255)
                                                                  gACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCCTCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGAGTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGc                                                                                    <  2:939678/150‑1 (MQ=255)
                                                                  gACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGc                                                                                    >  2:1631061/1‑150 (MQ=255)
                                                                   aCCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGcc                                                                                   <  1:1821977/150‑1 (MQ=255)
                                                                         ccGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGc                                                                             <  1:770687/150‑1 (MQ=255)
                                                                           ggTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCCTCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCgg                                                                           <  2:1280258/150‑1 (MQ=255)
                                                                           ggTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCgg                                                                           <  1:2527453/150‑1 (MQ=255)
                                                                           ggTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCgg                                                                           <  1:2692274/150‑1 (MQ=255)
                                                                            gTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCggg                                                                          <  1:968805/150‑1 (MQ=255)
                                                                                 cgccgcCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTttc                                                                     <  2:3450582/150‑1 (MQ=255)
                                                                                 cgccgcCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGCTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTttc                                                                     <  1:2603693/150‑1 (MQ=255)
                                                                                  gccgccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTttct                                                                    >  1:3973194/1‑150 (MQ=255)
                                                                                  gccgccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTttct                                                                    >  1:3761133/1‑150 (MQ=255)
                                                                                  gccgccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTttct                                                                    <  1:2971259/150‑1 (MQ=255)
                                                                                     gccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTcc                                                                 >  2:4033356/1‑150 (MQ=255)
                                                                                     gccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTcc                                                                 >  1:1626684/1‑150 (MQ=255)
                                                                                     gccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTcc                                                                 >  2:2036665/1‑150 (MQ=255)
                                                                                     gccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTcc                                                                 >  2:2508254/1‑150 (MQ=255)
                                                                                      ccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCa                                                                <  2:1579724/150‑1 (MQ=255)
                                                                                      ccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCa                                                                <  1:144134/150‑1 (MQ=255)
                                                                                      ccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCa                                                                <  2:2256880/150‑1 (MQ=255)
                                                                                            gggCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCatat                                                          <  1:3116391/150‑1 (MQ=255)
                                                                                            gggCAAACCAGCTGTCATCCATCCGCACTTCGACAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCatat                                                          >  2:3815784/1‑150 (MQ=255)
                                                                                             ggCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATc                                                         >  1:2037088/1‑150 (MQ=255)
                                                                                              gCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATcc                                                        >  1:4235130/1‑150 (MQ=255)
                                                                                               cAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCa                                                       <  1:817463/150‑1 (MQ=255)
                                                                                                  aCCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAgcg                                                    <  2:2749921/150‑1 (MQ=255)
                                                                                                     aGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac                                                 <  1:3074324/150‑1 (MQ=255)
                                                                                                       cTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcacca                                               <  1:301983/150‑1 (MQ=255)
                                                                                                        tGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcaccac                                              >  1:3332385/1‑150 (MQ=255)
                                                                                                         gTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCaccacc                                             <  1:351762/150‑1 (MQ=255)
                                                                                                         gTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCaccacc                                             <  1:2415224/150‑1 (MQ=255)
                                                                                                            atccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGACGCGGGTTTTCTTCCAGCATATCCAGCGCACCCCCaaa                                          >  2:2529626/1‑150 (MQ=255)
                                                                                                            atccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCaaa                                          >  1:437339/1‑150 (MQ=255)
                                                                                                            atccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCaaa                                          <  1:2796405/150‑1 (MQ=255)
                                                                                                             tccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCaaaa                                         <  1:78677/150‑1 (MQ=255)
                                                                                                              ccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCACAac                                        >  2:213623/1‑150 (MQ=255)
                                                                                                              ccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAac                                        >  2:244585/1‑150 (MQ=255)
                                                                                                              ccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAac                                        <  2:2166712/150‑1 (MQ=255)
                                                                                                              ccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAac                                        <  2:2196963/150‑1 (MQ=255)
                                                                                                              ccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAac                                        <  2:604903/150‑1 (MQ=255)
                                                                                                               catccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAacc                                       >  2:3765626/1‑150 (MQ=255)
                                                                                                               catccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTACCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAacc                                       >  1:1710354/1‑150 (MQ=255)
                                                                                                                atccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAacca                                      >  1:332188/1‑150 (MQ=255)
                                                                                                                 tccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccag                                     >  1:1788616/1‑150 (MQ=255)
                                                                                                                  ccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccaga                                    <  2:780032/150‑1 (MQ=255)
                                                                                                                   cGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagac                                   >  1:2059933/1‑150 (MQ=255)
                                                                                                                   cGCACTTCGTCAATGCCGTTAATCTGAGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCTGCGCACCACCAAAaccagac                                   <  2:2230291/150‑1 (MQ=255)
                                                                                                                     cACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagacca                                 >  2:3981193/1‑150 (MQ=255)
                                                                                                                            tCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACgg                          <  1:3981543/150‑1 (MQ=255)
                                                                                                                                  ccGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCa                    <  2:3619830/150‑1 (MQ=255)
                                                                                                                                  ccGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCa                    <  2:1510423/150‑1 (MQ=255)
                                                                                                                                        aTCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGt              <  2:1206350/150‑1 (MQ=255)
                                                                                                                                          cTGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCAAGAATATCCAGCGCACCACCAAAACCATACCAGTTACGGAACTCACCCAGTGc            >  1:1013098/1‑150 (MQ=255)
                                                                                                                                           tGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCg           >  2:497419/1‑150 (MQ=255)
                                                                                                                                           tGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCg           >  2:3800329/1‑150 (MQ=255)
                                                                                                                                           tGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCg           >  2:1003404/1‑150 (MQ=255)
                                                                                                                                             cGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTGTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTc         >  1:4197963/1‑150 (MQ=255)
                                                                                                                                                 aTAGGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCTCCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCa     >  1:3411655/1‑150 (MQ=255)
                                                                                                                                                 aTACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCa     >  2:1779197/1‑150 (MQ=255)
                                                                                                                                                  tACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCAc    <  2:3207777/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGA  >  NC_000913/3601358‑3601655

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: