Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 593,472 G→T 5.7% R433R (CGC→CGA cusS ← copper‑sensing histidine kinase in two‑component regulatory system with CusR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913593,4720GT 5.7% 522.2 / 14.6 140R433R (CGC→CGAcusScopper‑sensing histidine kinase in two‑component regulatory system with CusR
Reads supporting (aligned to +/- strand):  ref base G (67/65);  new base T (4/4);  total (71/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 143 total aligned reads displayed.

AAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTCGAAAG  >  NC_000913/593323‑593621
                                                                                                                                                     |                                                                                                                                                     
aaGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcg                                                                                                                                                       >  1:2297352/1‑150 (MQ=255)
   gATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGa                                                                                                                                                    >  1:2257575/1‑150 (MQ=255)
    aTTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGAt                                                                                                                                                   >  1:1779069/1‑150 (MQ=255)
     ttAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATa                                                                                                                                                  <  1:215413/150‑1 (MQ=255)
         gCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaaa                                                                                                                                              >  2:694775/1‑150 (MQ=255)
          cGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAAc                                                                                                                                             <  1:2288934/150‑1 (MQ=255)
           gggTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACg                                                                                                                                            <  2:1277584/150‑1 (MQ=255)
           gggTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACg                                                                                                                                            <  2:1277545/150‑1 (MQ=255)
                 tgtgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCaaa                                                                                                                                      <  2:388210/150‑1 (MQ=255)
                   tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACTCGATAGAAACGGTCAGACa                                                                                                                                    >  2:1749073/1‑150 (MQ=255)
                   tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACTCGATAGAAACGGTCAAACa                                                                                                                                    >  2:158046/1‑150 (MQ=255)
                   tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACTCGATAGAAACGGTCAAACa                                                                                                                                    >  2:1850620/1‑150 (MQ=255)
                   tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACa                                                                                                                                    >  2:273367/1‑150 (MQ=255)
                   tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACa                                                                                                                                    >  2:1479930/1‑150 (MQ=255)
                      tAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATc                                                                                                                                 <  2:1766062/150‑1 (MQ=255)
                      tAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATc                                                                                                                                 >  1:1302959/1‑150 (MQ=255)
                        aCAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATcgc                                                                                                                               >  1:1738701/1‑150 (MQ=255)
                        aCAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAAACGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATcgc                                                                                                                               >  1:361615/1‑150 (MQ=255)
                          aaaCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCgg                                                                                                                             >  1:1409851/1‑150 (MQ=255)
                             ccTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTaa                                                                                                                          <  1:964454/150‑1 (MQ=255)
                              cTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAg                                                                                                                         >  1:1277197/1‑150 (MQ=255)
                                     cccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCgg                                                                                                                  >  2:2294243/1‑150 (MQ=255)
                                     cccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCgg                                                                                                                  >  2:175647/1‑150 (MQ=255)
                                     cccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCgg                                                                                                                  >  2:1727135/1‑150 (MQ=255)
                                      ccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCggg                                                                                                                 >  2:2255785/1‑150 (MQ=255)
                                      ccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACAGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGACAAACAATCGCGGTAAGTGCTCggg                                                                                                                 >  1:1990521/1‑150 (MQ=255)
                                       cgcgCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgc                                                                                                                >  1:1570075/1‑150 (MQ=255)
                                        gcgcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgcg                                                                                                               >  2:1547278/1‑150 (MQ=255)
                                        gcgcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgcg                                                                                                               >  2:1168092/1‑150 (MQ=255)
                                            aTCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAAt                                                                                                           <  2:2208462/150‑1 (MQ=255)
                                            aTCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAAt                                                                                                           <  2:1589728/150‑1 (MQ=255)
                                              ccGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATgg                                                                                                         >  1:1745796/1‑150 (MQ=255)
                                              ccGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATgg                                                                                                         >  1:1218060/1‑150 (MQ=255)
                                              ccGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATgg                                                                                                         >  1:1664589/1‑150 (MQ=255)
                                               cGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATggg                                                                                                        <  2:1414412/150‑1 (MQ=255)
                                               cGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATggg                                                                                                        <  2:1976048/150‑1 (MQ=255)
                                                 aCGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCg                                                                                                      <  1:553665/150‑1 (MQ=255)
                                                   gTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTa                                                                                                    >  2:1255017/1‑150 (MQ=255)
                                                      aCCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCg                                                                                                 <  2:2001172/150‑1 (MQ=255)
                                                        cGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGACGCAATGGGCGTACCggg                                                                                               >  1:2341104/1‑150 (MQ=255)
                                                          cAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGt                                                                                             >  1:1948429/1‑150 (MQ=255)
                                                          cAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGt                                                                                             >  1:1330642/1‑150 (MQ=255)
                                                             ccGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGtttt                                                                                          <  2:2100379/150‑1 (MQ=255)
                                                                   cTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTcgacga                                                                                    >  1:1472272/1‑150 (MQ=255)
                                                                      tATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATaa                                                                                 <  1:273367/150‑1 (MQ=255)
                                                                      tATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATaa                                                                                 <  1:1479930/150‑1 (MQ=255)
                                                                       aTGGGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCACTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAAc                                                                                <  2:265061/150‑1 (MQ=255)
                                                                       aTGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAAc                                                                                <  2:265062/150‑1 (MQ=255)
                                                                                aacaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAg                                                                       >  1:2329947/1‑150 (MQ=255)
                                                                                aacaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAg                                                                       >  1:1143055/1‑150 (MQ=255)
                                                                                 acaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAgg                                                                      <  1:740509/150‑1 (MQ=255)
                                                                                 acaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAgg                                                                      <  1:991201/150‑1 (MQ=255)
                                                                                     tCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGAt                                                                  >  2:1398254/1‑150 (MQ=255)
                                                                                     tCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGAt                                                                  >  2:1398178/1‑150 (MQ=255)
                                                                                      cGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACTCGATAGAAACGGTCAGACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATc                                                                 <  1:1749073/150‑1 (MQ=255)
                                                                                      cGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACTCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATc                                                                 <  1:1850620/150‑1 (MQ=255)
                                                                                      cGATGTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACTCGATAGAAACGGTCTAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATc                                                                 <  1:360946/150‑1 (MQ=255)
                                                                                        aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGa                                                               <  1:1381640/150‑1 (MQ=255)
                                                                                        aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGa                                                               <  1:1711443/150‑1 (MQ=255)
                                                                                                aTCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGa                                                       <  1:694775/150‑1 (MQ=255)
                                                                                                      aGGCCAATACCGCTACCTCCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAgcgc                                                 >  2:774677/1‑150 (MQ=255)
                                                                                                        gCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCAc                                               >  1:757224/1‑150 (MQ=255)
                                                                                                         ccAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACt                                              >  1:1057701/1‑150 (MQ=255)
                                                                                                          cAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTa                                             <  1:2255785/150‑1 (MQ=255)
                                                                                                           aaTAGCGCTACCTTCTGGTTTTCGCTGGGGGGAGGGGTGAACGCGATAGAAACGGTCAAACTATCGCGGTAAGTGGTCGGGCGCAATGGGGGTTCCGGGGTTTTCGTCGTTATGTTGCTGCAGGTGATCGTCCGTGTGTCTgcgctgtag                                            >  2:1917044/1‑145 (MQ=255)
                                                                                                           aaTACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTAc                                            >  1:588973/1‑150 (MQ=255)
                                                                                                              aCCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAAt                                         <  2:283108/150‑1 (MQ=255)
                                                                                                                cGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTg                                       >  1:623519/1‑150 (MQ=255)
                                                                                                                 gCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTGGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGt                                      <  1:2152104/150‑1 (MQ=255)
                                                                                                                 gCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGt                                      <  1:378858/150‑1 (MQ=255)
                                                                                                                 gCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGt                                      <  1:378836/150‑1 (MQ=255)
                                                                                                                 gCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGt                                      <  1:378825/150‑1 (MQ=255)
                                                                                                                  cTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtc                                     >  1:1834225/1‑150 (MQ=255)
                                                                                                                   tACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtct                                    >  1:1583010/1‑150 (MQ=255)
                                                                                                                      cTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtctctc                                 <  2:2297352/150‑1 (MQ=255)
                                                                                                                           ccTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGtt                            >  2:459950/1‑150 (MQ=255)
                                                                                                                           ccTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGtt                            <  2:866172/150‑1 (MQ=255)
                                                                                                                           ccTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGtt                            <  2:868180/150‑1 (MQ=255)
                                                                                                                           ccTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGtt                            >  2:2097652/1‑150 (MQ=255)
                                                                                                                           ccTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGtt                            <  1:1092399/150‑1 (MQ=255)
                                                                                                                           ccTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGtt                            <  1:2108013/150‑1 (MQ=255)
                                                                                                                             ttttCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTgg                          >  1:1655856/1‑150 (MQ=255)
                                                                                                                              tttCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGc                         <  2:692523/150‑1 (MQ=255)
                                                                                                                                 cGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGta                      <  1:1540870/150‑1 (MQ=255)
                                                                                                                                   cTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCTATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGtata                    <  1:1970098/150‑1 (MQ=255)
                                                                                                                                    tGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATaa                   <  2:1218060/150‑1 (MQ=255)
                                                                                                                                        gggAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCa               <  1:2200721/150‑1 (MQ=255)
                                                                                                                                        gggAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCa               <  1:2200610/150‑1 (MQ=255)
                                                                                                                                          gAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAgg             >  1:781229/1‑150 (MQ=255)
                                                                                                                                               gTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGtt        >  1:106615/1‑150 (MQ=255)
                                                                                                                                               gTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGtt        <  2:2031581/150‑1 (MQ=255)
                                                                                                                                               gTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGtt        >  2:400343/1‑150 (MQ=255)
                                                                                                                                               gTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGtt        <  1:522097/150‑1 (MQ=255)
                                                                                                                                               gTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGtt        <  2:2031579/150‑1 (MQ=255)
                                                                                                                                                tCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCCGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTc       <  2:2330952/150‑1 (MQ=255)
                                                                                                                                                tCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTc       >  2:1913118/1‑150 (MQ=255)
                                                                                                                                                   aCGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTCGaa    >  1:1967879/1‑150 (MQ=255)
                                                                                                                                                   aCGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTCGaa    >  1:595366/1‑150 (MQ=255)
                                                                                                                                                     gcgATAGAACCGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTCGAAAg  <  2:75107/150‑1 (MQ=255)
                                                                                                                                                     gcgATAGAAACGGTCACACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTCGAAAg  <  2:1708998/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTCGAAAG  >  NC_000913/593323‑593621

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: