Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 182,661 G→C 5.9% V67L (GTA→CTA)  cdaR → DNA‑binding transcriptional activator CdaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913182,6610GC 5.9% 576.4 / 19.3 169V67L (GTA→CTA) cdaRDNA‑binding transcriptional activator CdaR
Reads supporting (aligned to +/- strand):  ref base G (89/70);  new base C (5/5);  total (94/75)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

Only 100 of 173 total aligned reads displayed.

GTACCATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAATGATGCTG  >  NC_000913/182512‑182810
                                                                                                                                                     |                                                                                                                                                     
gTACCATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCgg                                                                                                                                                       >  1:1046661/1‑150 (MQ=255)
gTACCATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCgg                                                                                                                                                       >  2:755954/1‑150 (MQ=255)
 tACCATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGt                                                                                                                                                      <  2:1385198/150‑1 (MQ=255)
    cATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGc                                                                                                                                                   >  1:2244250/1‑150 (MQ=255)
     aTGCGCACCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCa                                                                                                                                                  <  1:670421/150‑1 (MQ=255)
      tGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCAc                                                                                                                                                 >  2:790732/1‑150 (MQ=255)
      tGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCAc                                                                                                                                                 >  1:158073/1‑150 (MQ=255)
       gcgcATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACg                                                                                                                                                <  2:1130437/150‑1 (MQ=255)
       gcgcATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACg                                                                                                                                                <  1:340868/150‑1 (MQ=255)
         gcATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTc                                                                                                                                              >  1:161156/1‑150 (MQ=255)
          catcatCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCa                                                                                                                                             >  1:8447/1‑150 (MQ=255)
          catcatCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCa                                                                                                                                             >  1:1784347/1‑150 (MQ=255)
           atcatcGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGTGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCAt                                                                                                                                            >  1:1304006/1‑150 (MQ=255)
           atcatcGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCAt                                                                                                                                            <  1:1814166/150‑1 (MQ=255)
               tcGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGTGCTACCGCGTCATcgat                                                                                                                                        <  1:2062625/150‑4 (MQ=255)
                 gATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCAc                                                                                                                                      <  2:1981431/150‑1 (MQ=255)
                   tACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACgg                                                                                                                                    <  1:2228210/150‑1 (MQ=255)
                   tACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACgg                                                                                                                                    <  1:2444641/150‑1 (MQ=255)
                   tACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACgg                                                                                                                                    <  2:2221608/150‑1 (MQ=255)
                     ccAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGgtg                                                                                                                                  >  1:14085/1‑150 (MQ=255)
                     ccAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGgtg                                                                                                                                  <  1:722941/150‑1 (MQ=255)
                     ccAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCAGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGgtg                                                                                                                                  >  2:1261504/1‑150 (MQ=255)
                        atatCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGc                                                                                                                               >  1:2374209/1‑150 (MQ=255)
                          atCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCgg                                                                                                                             <  1:2389570/150‑1 (MQ=255)
                           tCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGAGATCGTGAGCGTATTGGTGAATTTCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGc                                                                                                                            >  2:441585/1‑150 (MQ=255)
                             aaCGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAg                                                                                                                          >  1:1080169/1‑150 (MQ=255)
                             aaCGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCAATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAg                                                                                                                          <  2:1886478/150‑1 (MQ=255)
                              aCGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAgg                                                                                                                         >  2:1343226/1‑150 (MQ=255)
                                gTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAgggg                                                                                                                       >  2:705683/1‑150 (MQ=255)
                                gTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAgggg                                                                                                                       <  2:2288054/150‑1 (MQ=255)
                                    tGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTa                                                                                                                   >  2:245838/1‑150 (MQ=255)
                                       aTGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATc                                                                                                                >  2:598267/1‑150 (MQ=255)
                                          cccGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTAc                                                                                                             <  2:1718080/150‑1 (MQ=255)
                                              ttggCGAATTATCGGCAGCGGCGATCGTGAGCGTATTTGTGAATTTCACGAAGGTTCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGtt                                                                                                         <  1:2126509/148‑1 (MQ=255)
                                                  cGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACgg                                                                                                     >  2:859701/1‑150 (MQ=255)
                                                  cGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACgg                                                                                                     <  1:790732/150‑1 (MQ=255)
                                                  cGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACgg                                                                                                     >  2:2240517/1‑150 (MQ=255)
                                                     aTTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTg                                                                                                  <  1:2074687/150‑1 (MQ=255)
                                                        aTCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGaa                                                                                               >  1:164656/1‑150 (MQ=255)
                                                          cGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAgg                                                                                             <  2:1320911/150‑1 (MQ=255)
                                                           ggCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGTGCTACCGCGTCATc                                                                                                                                           >  1:1970043/1‑103 (MQ=255)
                                                           ggCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGTGCTACCGCGTCATc                                                                                                                                           <  2:1970043/103‑1 (MQ=255)
                                                           ggCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGt                                                                                            >  1:2381947/1‑150 (MQ=255)
                                                                cggcgATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCt                                                                                                               <  2:2213327/126‑1 (MQ=255)
                                                                cggcgATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTAAGGCTGGAAGGTGAAAt                                                                                       >  1:2382961/1‑150 (MQ=255)
                                                                     aTCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCg                                                                                  <  2:1949205/150‑1 (MQ=255)
                                                                        gTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTAGTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCg                                                                               >  1:412528/1‑150 (MQ=255)
                                                                        gTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCg                                                                               >  1:261456/1‑150 (MQ=255)
                                                                           aGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTaa                                                                            >  2:498045/1‑150 (MQ=255)
                                                                           aGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTaa                                                                            >  2:506993/1‑150 (MQ=255)
                                                                             cGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAAtt                                                                          >  2:367434/1‑150 (MQ=255)
                                                                                aTTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGc                                                                       <  2:923059/150‑1 (MQ=255)
                                                                                aTTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGc                                                                       <  2:1080169/150‑1 (MQ=255)
                                                                                 ttGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGcc                                                                      <  1:2189736/150‑1 (MQ=255)
                                                                                 ttGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCCGGAAGGTGAAATTGTCGGCGTAATTGGcc                                                                      <  2:1852171/150‑1 (MQ=255)
                                                                                  tGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCt                                                                     <  2:1003955/150‑1 (MQ=255)
                                                                                   ggTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTg                                                                    >  2:2058021/1‑150 (MQ=255)
                                                                                     tGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGAc                                                                  <  1:598267/150‑1 (MQ=255)
                                                                                     tGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGAc                                                                  >  1:1738830/1‑150 (MQ=255)
                                                                                     tGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGAc                                                                  <  2:1806050/150‑1 (MQ=255)
                                                                                           gCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGa                                                            >  1:1220581/1‑150 (MQ=255)
                                                                                           gCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGa                                                            >  2:1167387/1‑150 (MQ=255)
                                                                                            cACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGaa                                                           >  1:1939737/1‑150 (MQ=255)
                                                                                              cGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAAcc                                                         >  2:1272014/1‑150 (MQ=255)
                                                                                              cGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAAcc                                                         <  1:441585/150‑1 (MQ=255)
                                                                                                  ggTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCagag                                                     >  1:693758/1‑150 (MQ=255)
                                                                                                  ggTGCATTGCTGGTACTTTCACAGGGACGAGTCGTAGATATCGATGACGCGGTAGCACGTAATCTGCACGGTGTGCGGCAGTGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTTATTGGCCTGACAGGTGAACCagag                                                     <  2:1784347/150‑1 (MQ=255)
                                                                                                    tGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGaa                                                   >  1:2054808/1‑150 (MQ=255)
                                                                                                      cATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATc                                                 >  1:2177223/1‑150 (MQ=255)
                                                                                                      cATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATc                                                 >  1:1833122/1‑150 (MQ=255)
                                                                                                        ttGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTg                                               <  1:1379656/150‑1 (MQ=255)
                                                                                                        ttGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTg                                               <  1:1864540/150‑1 (MQ=255)
                                                                                                           cTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGt                                            >  1:1469208/1‑150 (MQ=255)
                                                                                                           cTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGt                                            >  1:355741/1‑150 (MQ=255)
                                                                                                             ggTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTaa                                          >  1:89435/1‑150 (MQ=255)
                                                                                                                aCTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAata                                       <  2:529325/150‑1 (MQ=255)
                                                                                                                 cTTTCACAGGGACGAGTAGTCCAAATCGATCACGCGGTAGCACGTCATCTGCACCGTGTGCGGCAGGGGATTAATCTAGCGTTACGGCTGGAAGGTGAAATTGTCGGAGTTAATGTCCTGACAGATGAACCAGAGAATCTCCGTAAatat                                      >  2:234172/1‑150 (MQ=255)
                                                                                                                     cacaGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCg                                  <  1:403453/150‑1 (MQ=255)
                                                                                                                       caGGGAAGAGACGTCGATATTGATGACGCGGTAGAAAGTAAACTGCACGGTGTCAGGCAGGGCAAAAAAATACCATTACGACAAGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGaa                                >  2:2025063/1‑150 (MQ=255)
                                                                                                                        aGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAAc                               <  2:161156/150‑1 (MQ=255)
                                                                                                                         gggACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACt                              <  1:2281447/150‑1 (MQ=255)
                                                                                                                             cGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTc                          <  1:2240517/150‑1 (MQ=255)
                                                                                                                                gtcgtcGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGc                       >  1:867718/1‑150 (MQ=255)
                                                                                                                                gtcgtcGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGc                       >  2:2202311/1‑150 (MQ=255)
                                                                                                                                gtcgtcGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGc                       >  1:478180/1‑150 (MQ=255)
                                                                                                                                 tcgtcgATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCa                      <  1:1261504/150‑1 (MQ=255)
                                                                                                                                  cgtcgATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCAt                     >  1:53709/1‑150 (MQ=255)
                                                                                                                                  cgtcgATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCAt                     >  1:2049498/1‑150 (MQ=255)
                                                                                                                                     cgATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGAc                  >  2:2165526/1‑150 (MQ=255)
                                                                                                                                      gATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACg                 >  1:1705612/1‑150 (MQ=255)
                                                                                                                                          tCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTg             <  2:1677615/150‑1 (MQ=255)
                                                                                                                                          cagATGACGTGCTACCGCGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTg             >  2:2413761/3‑150 (MQ=255)
                                                                                                                                           cGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGa            >  1:559405/1‑150 (MQ=255)
                                                                                                                                           agATGACGTGCTACCGCGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAata                                       >  1:1782634/2‑123 (MQ=255)
                                                                                                                                              tGACGTGCTACCGCGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAat         >  2:309598/1‑150 (MQ=255)
                                                                                                                                              tGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAat         <  2:1750122/150‑1 (MQ=255)
                                                                                                                                                 cgcgGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAatgat      >  1:68490/1‑150 (MQ=255)
                                                                                                                                                   cgGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAATGATGc    <  2:1151698/150‑1 (MQ=255)
                                                                                                                                                    ggTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATATACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAATGATGCt   >  2:1662863/1‑150 (MQ=255)
                                                                                                                                                     gTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAATGATGCTg  <  1:1683624/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTACCATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAATGATGCTG  >  NC_000913/182512‑182810

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