Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,454 A→C 7.5% intergenic (‑97/+19) cutA ← / ← dcuA divalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,4540AC 7.5% 209.8 / 9.9 85intergenic (‑97/+19)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (42/35);  new base C (4/3);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

CCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTG  >  NC_000913/4365306‑4365601
                                                                                                                                                    |                                                                                                                                                   
ccGAACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGGCTGCGGGAAAGAAAGCGGGGCGTGa                                                                                                                                                     <  1:350349/149‑1 (MQ=255)
 cGTACATTGCCCCACGACAGACGCGGTACTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGTAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                    <  1:173725/149‑1 (MQ=255)
 cGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                    <  1:746743/149‑1 (MQ=255)
 cGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                    <  1:720111/149‑1 (MQ=255)
 cGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAAGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                    <  1:774609/149‑1 (MQ=255)
    aCATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCCTAAACACCTCGCAAGAACAGATGGAACCGCAAGGGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGc                                                                                                                                                 <  2:321740/149‑1 (MQ=255)
    aCATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGc                                                                                                                                                 <  1:224270/149‑1 (MQ=255)
    aCATACCACCACGACAGACTCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGc                                                                                                                                                 <  2:346268/149‑1 (MQ=255)
        ggccccgcGACCGGCGCGCTAATCGATCTTTTTCCATCACTCATATACACCTCCCAGGAACCGTTGCAACTGAAAGAGAGAGGTCACTATCCAACTGACGGGTGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcg                                                                                                                                             <  1:737824/142‑1 (MQ=255)
        aGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcg                                                                                                                                             >  2:524958/1‑149 (MQ=255)
         gCACCACGTCAGACGCGGTATTCGAACGTTTTTCATCAAGCATATACACCTCGCAATAACAGAGGGAACCGCAAGAGAGAGGTCACTACACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgc                                                                                                                                            <  2:28503/149‑1 (MQ=255)
          caccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGCCGGCTGCTAGTCTGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgcg                                                                                                                                           >  1:603826/1‑149 (MQ=255)
          caccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAGGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgcg                                                                                                                                           >  1:668562/1‑149 (MQ=255)
           accacGACAGACGCGGTATTCGAACTTTTTTCTTCTTTTTTAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGa                                                                                                                                          >  2:680421/1‑149 (MQ=255)
           accacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctt                                                                                                                                          <  2:42938/149‑3 (MQ=255)
            ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGGGGAGGCCCCGCGag                                                                                                                                         >  2:182380/1‑148 (MQ=255)
            ccacGACAGACGCGGTATTCGAACCTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGaa                                                                                                                                         >  1:583286/1‑148 (MQ=255)
               cGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttctt                                                                                                                                      <  1:137971/149‑7 (MQ=255)
               cGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCACCGCGATTGa                                                                                                                                      >  1:725503/1‑149 (MQ=255)
                gACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGGGCACGCCCCGCGATTGAt                                                                                                                                     >  2:54675/1‑149 (MQ=255)
                gACAGACGCGGTAGTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGGCACTATACACCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGAt                                                                                                                                     <  1:742934/149‑1 (MQ=255)
                       gcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGGAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACa                                                                                                                              <  1:362103/149‑1 (MQ=255)
                        cgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACCCCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                             <  1:734504/149‑1 (MQ=255)
                         ggTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGc                                                                                                                            <  1:265564/149‑1 (MQ=255)
                          gTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCa                                                                                                                           >  1:721618/1‑149 (MQ=255)
                           tATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCAt                                                                                                                          >  2:762333/1‑149 (MQ=255)
                             ttCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCACGAGAGAGGGCACTATACAACGGACGGGGGAAGGATGTCCGTCGGCAGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGa                                                                                                                        >  1:568960/1‑149 (MQ=255)
                             ttCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGAGGGGCGTGAACGCCCCGCGATTGATCATTACAGCAGGa                                                                                                                        >  1:290297/1‑149 (MQ=255)
                              tCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGaa                                                                                                                       <  2:472366/149‑1 (MQ=255)
                               cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                                                                      <  2:176286/149‑1 (MQ=255)
                                gAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                     >  1:784610/1‑149 (MQ=255)
                                gAACTTTTTTCATCAAGCATAAACACATCGCAAGAACAGATGGAACCGCAAGAGAGAGGGCACTATACAAAGGACGGGGGAAGTATGTTGATCGGAGGGAAAGAAAGCGCGGCGCGAACGCCCCGCGATGGGGCAGTACCGCAGGAAgg                                                                                                                     >  2:336722/1‑148 (MQ=255)
                                 aaCTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCGGCGATTGATCATTACAGCATGAAgat                                                                                                                    >  2:298982/1‑147 (MQ=255)
                                  aCTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGAGCATTACAGCATGAAGCTa                                                                                                                   >  2:521766/1‑149 (MQ=255)
                                           atcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCCTGAAGCTACCCAGCAcg                                                                                                          >  2:288664/1‑149 (MQ=255)
                                           atcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTGCCCAGCAcg                                                                                                          >  1:395235/1‑149 (MQ=255)
                                            tcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcga                                                                                                         >  1:151221/1‑149 (MQ=255)
                                               aGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACGACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAGTGATCATTACAGCATGAGCCTACCCAGCCccaagc                                                                                                      >  1:338331/1‑149 (MQ=255)
                                                 cATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGCAAGTAGACGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccg                                                                                                    >  1:570683/1‑149 (MQ=255)
                                                        acCTCGCAAGAACAGAGGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGAGGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATCGATCATTACAGCATGAAGCTACCCAGCACGAACCCGAAGCaaa                                                                                             >  1:653020/1‑149 (MQ=255)
                                                          ctctCAAGAACAAATGGAACCGCAAGAGACAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACg                                                                                           <  1:1482/149‑1 (MQ=255)
                                                          cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACg                                                                                           <  1:639246/149‑1 (MQ=255)
                                                          cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACg                                                                                           <  1:311604/149‑1 (MQ=255)
                                                            cGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGGAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATGACAGCATGAAGCTACCCAGCAcgaagccgaagccacaggc                                                                                         >  2:336901/1‑142 (MQ=255)
                                                               aaGAACAGAGGGAACCGCAAGAGACAGGTCACTATACAACGGACGGGTGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTAACGCCCCGCGATTGTTCATTACATCATGAAGGTACCCAGCACGAAGCCGAAGCAAACGGCCac                                                                                      <  2:572426/149‑2 (MQ=255)
                                                                     aGATGGAACCGCAAGAGAGAGGTCACTAAACAGCGGACGGGGGAAGGATGATAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAac                                                                                <  1:388352/149‑1 (MQ=255)
                                                                       aTGGAACCGCAAGAGAGAGGTCACTATACATGGGGCGGGGGAAGGATGTTAGTCGGCGGGGAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAacac                                                                              >  1:776912/1‑149 (MQ=255)
                                                                          gAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCCTTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAAGACCCa                                                                           >  1:76487/1‑149 (MQ=255)
                                                                          gAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                                                                           >  1:750882/1‑149 (MQ=255)
                                                                           aaCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGACAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCGAGGACGAAGCCGAAGCAAACGGCCAGGGCGccagccag                                                                          >  1:652371/1‑149 (MQ=255)
                                                                           aaCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGGTGGATCATTACAGCCGGAAGCTACCCAGCACGAAGCCGAAGGAAACGGGCAGGGCAACACGCag                                                                          >  2:594596/1‑149 (MQ=255)
                                                                           aaCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                                                                          >  2:288826/1‑149 (MQ=255)
                                                                            aCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCaga                                                                         >  1:724836/1‑149 (MQ=255)
                                                                                    cgagagGTCACTATACCACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCgg                                                                 <  1:143796/148‑1 (MQ=255)
                                                                                        agGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAAGACCCAGAGTACCCGGGATg                                                             >  1:50849/1‑149 (MQ=255)
                                                                                         ggTGACTATACAACGGACGGGGGAAGGATGTTGGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCCATTGGTCTTTACAGCATGAAGCGACCCAGCAGGAAGGCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATga                                                            <  2:264570/149‑1 (MQ=255)
                                                                                             aCTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaagaa                                                        >  1:646495/1‑149 (MQ=255)
                                                                                             aCTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaagaa                                                        <  2:724836/149‑1 (MQ=255)
                                                                                              ctatcaaaCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAAc                                                       <  1:202069/143‑1 (MQ=255)
                                                                                              cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGGTGAAGAAc                                                       >  1:243473/1‑149 (MQ=255)
                                                                                                 taCAACGGACGGGGGAAGGATGTTAGTCGGTGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGa                                                    <  1:54675/149‑1 (MQ=255)
                                                                                                      cggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGtt                                               <  1:70128/149‑1 (MQ=255)
                                                                                                      cggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGtt                                               <  2:620519/149‑1 (MQ=255)
                                                                                                       ggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGCAACTACCCAGCACGAAGCCGAAGCAAACGGCCAGCGCCACACCCAGAGTACCCGGGATGAAGACCGGATCGTTg                                              >  1:374973/1‑149 (MQ=255)
                                                                                                        gacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGa                                             <  2:784610/149‑1 (MQ=255)
                                                                                                         gcggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGaa                                            <  1:648973/148‑1 (MQ=255)
                                                                                                         ccggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGaa                                            <  2:241483/148‑1 (MQ=255)
                                                                                                         acggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCCACACCCAGAGTCCCCGGGATGAAGAACGGAAGGTTGaa                                            >  1:214256/1‑149 (MQ=255)
                                                                                                           gggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGa                                          <  2:716521/149‑1 (MQ=255)
                                                                                                           gggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGa                                          <  2:531214/149‑1 (MQ=255)
                                                                                                             gggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACg                                        >  1:329721/1‑149 (MQ=255)
                                                                                                               gAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGaa                                      >  1:284869/1‑149 (MQ=255)
                                                                                                                  ggATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCACCACCCAGAGTACCCGGGATGAAGAACGGAAGGTTGAAGAAGAAttt                                   >  2:291313/1‑149 (MQ=255)
                                                                                                                   ggtGTTGGTCGGCGGGAAAGAATGCGGGGCGTGAGCGCCCCGCAATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCTGGGCAACACCCAGAGTACCCGGGTTGAAGAACGGATGGTTGAAGACGAATTTa                                  <  2:429980/147‑1 (MQ=255)
                                                                                                                     cgtgaGTCGGCGGGAAAGAAAGCGGGGGGGGACCGCCCCGCGATTGACCATGACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAcc                                <  1:280130/145‑1 (MQ=255)
                                                                                                                                gggAAAGCAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACACCGAACGAAGACGAAGCAAACGGCCAGTGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGACGAAGAATTTACCGCTCCGGgtag                     >  1:542757/1‑149 (MQ=255)
                                                                                                                                gggAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTCCAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAACCGCGAAGAAGAACGGAAGGTTGAAGACGAATTTACCGATACGGgtaa                     >  1:358493/1‑148 (MQ=255)
                                                                                                                                gggAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCATCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAAGGTTGAAGACGAATTTACCCATCCGGgtag                     >  1:324878/1‑149 (MQ=255)
                                                                                                                                gggAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAGGACGAATTTACCGATACGGgtag                     <  1:413237/149‑1 (MQ=255)
                                                                                                                                  gaaagaaaGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgtagta                   <  2:725503/149‑1 (MQ=255)
                                                                                                                                   aaagaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTAc                  >  1:297914/1‑149 (MQ=255)
                                                                                                                                     agaaagCGGGGCGTAACAGCCAAACCATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTAccc                <  1:336901/149‑1 (MQ=255)
                                                                                                                                      aatcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg               >  1:526579/5‑149 (MQ=255)
                                                                                                                                         agCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtgt            <  2:603826/149‑1 (MQ=255)
                                                                                                                                         agCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGCAGTACCCgtgg            >  2:647602/1‑148 (MQ=255)
                                                                                                                                          gCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGGATTTACCGATACGGGTAGTACCCgttgc           >  2:12435/1‑146 (MQ=255)
                                                                                                                                          gCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATCCGGGTAGTGGCCGGGTc           >  2:546788/1‑149 (MQ=255)
                                                                                                                                              ggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCTGTGTcatcc       <  2:614150/149‑1 (MQ=255)
                                                                                                                                                cGTGACCGCCCCGCGATTGATCATTACAGCATGACGCTACCCAGCATGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatccat     <  1:216806/149‑1 (MQ=255)
                                                                                                                                                   gAACGCCCCGCGATTGATCATTACACCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACACGTGACATCCATCTg  >  1:393016/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTG  >  NC_000913/4365306‑4365601

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