Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 928,308 T→G 5.1% P296P (CCA→CCC)  cydC ← glutathione/L‑cysteine ABC exporter subunit CydC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913928,3080TG 5.1% 239.6 / 4.5 78P296P (CCA→CCC) cydCglutathione/L‑cysteine ABC exporter subunit CydC
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base G (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCG  >  NC_000913/928161‑928456
                                                                                                                                                   |                                                                                                                                                    
ctgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTgg                                                                                                                                                    >  1:847203/1‑150 (MQ=255)
ctgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTgg                                                                                                                                                    <  2:651979/150‑1 (MQ=255)
   aacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGc                                                                                                                                                 <  2:751333/150‑1 (MQ=255)
    acATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGcc                                                                                                                                                >  2:1458522/1‑150 (MQ=255)
        cccGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTg                                                                                                                                            >  2:1316166/1‑150 (MQ=255)
         ccGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGCTACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                                                                                                                           >  1:775794/1‑150 (MQ=255)
         ccGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                                                                                                                           >  1:416351/1‑150 (MQ=255)
          cgaaaCGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCt                                                                                                                                          >  1:875286/4‑150 (MQ=255)
          cGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCt                                                                                                                                          <  1:1418918/150‑1 (MQ=255)
          cGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCt                                                                                                                                          <  2:1567117/150‑1 (MQ=255)
          cGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCt                                                                                                                                          >  1:875287/1‑150 (MQ=255)
                gTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAc                                                                                                                                    <  2:847203/150‑1 (MQ=255)
                gTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAc                                                                                                                                    >  1:335177/1‑150 (MQ=255)
                gTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAc                                                                                                                                    >  1:341098/1‑150 (MQ=255)
                gTCAGCGAAACGCGATCGGCAACACGAGTTTCGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAc                                                                                                                                    <  2:1460786/150‑1 (MQ=255)
                   aGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAcgcg                                                                                                                                 <  2:1402069/150‑1 (MQ=255)
                     cGAAACGCGATCGGCAACACGAGTTTGGGTATCAGCAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGCTGCCAGTGCTTCAAACGCGGc                                                                                                                               >  1:1390420/1‑150 (MQ=255)
                      gAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCt                                                                                                                              <  2:1341580/150‑1 (MQ=255)
                        aaCGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTaa                                                                                                                            >  1:1367162/1‑150 (MQ=255)
                          cgcgATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAAcg                                                                                                                          <  2:1213607/150‑1 (MQ=255)
                             gATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACgcgc                                                                                                                       <  1:1455350/150‑1 (MQ=255)
                               tCGGCAACACGAGTCTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGGCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATCCACCCGCTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCag                                                                                                                     <  2:1239662/150‑1 (MQ=255)
                                cGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCaga                                                                                                                    >  2:978548/1‑150 (MQ=255)
                                 ggCAACACGAGTTGGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCagaa                                                                                                                   <  1:541792/150‑1 (MQ=255)
                                  gCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCagaag                                                                                                                  <  1:284578/150‑1 (MQ=255)
                                      cacGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACaa                                                                                                              >  1:622520/1‑150 (MQ=255)
                                          aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCGGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAg                                                                                                          >  2:1553601/1‑150 (MQ=255)
                                               gggTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCaa                                                                                                     <  2:923979/150‑1 (MQ=255)
                                                  tATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTa                                                                                                  >  2:373751/1‑150 (MQ=255)
                                                  tATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTa                                                                                                  >  2:1020492/1‑150 (MQ=255)
                                                   aTCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTaa                                                                                                 <  1:250357/150‑1 (MQ=255)
                                                      aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAAcgc                                                                                              >  2:304971/1‑150 (MQ=255)
                                                       ggAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAAcgcg                                                                                             <  1:616097/150‑1 (MQ=255)
                                                        gAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACgcgc                                                                                            >  1:1560179/1‑150 (MQ=255)
                                                          aaGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCg                                                                                          <  2:696526/150‑1 (MQ=255)
                                                          aaGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACACGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCg                                                                                          >  2:231979/1‑150 (MQ=255)
                                                            ggTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCggg                                                                                        >  1:192266/1‑150 (MQ=255)
                                                               gACCTCGGGTTTTTGATCCGTTAAGTCAGCGATACGACCGGCAGAGGCAATGACTTGCCCCAGATTCTGAAATGCACCCGTTTCTTTTGCCAGTGCTTCAAACGCGGCTAACGCGCTGAAGACAAACAGGGCTATTAACGCGCCGGGTTg                                                                                     <  2:1449763/150‑1 (MQ=255)
                                                               gACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTg                                                                                     >  2:252960/1‑150 (MQ=255)
                                                                aCCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGa                                                                                    <  2:579843/150‑1 (MQ=255)
                                                                 ccTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAg                                                                                   <  1:1316166/150‑1 (MQ=255)
                                                                 ccTCCGG‑TTTTGATCCGTTAAGTCATAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGc                                                                                  <  2:1060109/150‑1 (MQ=255)
                                                                  cTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGc                                                                                  >  2:379228/1‑150 (MQ=255)
                                                                         ttttGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTgccg                                                                           >  2:865405/1‑150 (MQ=255)
                                                                           ttGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTgccgcc                                                                         <  1:1510422/150‑1 (MQ=255)
                                                                            tGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCa                                                                        >  1:672215/1‑150 (MQ=255)
                                                                                  gTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAAcgccg                                                                  <  2:775794/150‑1 (MQ=255)
                                                                                   ttAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAAcgccgc                                                                 <  1:304971/150‑1 (MQ=255)
                                                                                    tAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACgccgcc                                                                <  2:54873/150‑1 (MQ=255)
                                                                                    tAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACgccgcc                                                                >  2:651530/1‑150 (MQ=255)
                                                                                    tAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACgccgcc                                                                >  1:878374/1‑150 (MQ=255)
                                                                                     aaGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCa                                                               >  1:806881/1‑150 (MQ=255)
                                                                                           gagaTACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGcc                                                         >  1:86315/1‑150 (MQ=255)
                                                                                               tacgtacgGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATcc                                                     <  1:745941/150‑1 (MQ=255)
                                                                                               tacgtacgGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATcc                                                     <  2:806881/150‑1 (MQ=255)
                                                                                                 cgtacgGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCcac                                                   >  2:580100/1‑150 (MQ=255)
                                                                                                          gagGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAgg                                          <  2:1404819/150‑1 (MQ=255)
                                                                                                            ggCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGAt                                        >  1:1125683/1‑150 (MQ=255)
                                                                                                              cAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGTATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCa                                      >  1:1039767/1‑150 (MQ=255)
                                                                                                                aTGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCAcc                                    >  1:188735/1‑150 (MQ=255)
                                                                                                                    cTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTACCGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTa                                <  2:878374/150‑1 (MQ=255)
                                                                                                                       gCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAAcg                             >  2:874229/1‑150 (MQ=255)
                                                                                                                           cAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGcc                         <  2:416351/150‑1 (MQ=255)
                                                                                                                                 cTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag                   >  2:96312/1‑150 (MQ=255)
                                                                                                                                 cTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag                   <  2:763345/150‑1 (MQ=255)
                                                                                                                                 agcactggCACCAGTAACGGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag                   >  2:1275841/8‑150 (MQ=255)
                                                                                                                                   gAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGagca                 <  1:379228/150‑1 (MQ=255)
                                                                                                                                    actggCACCAGTAACGGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATg                     <  2:1174956/141‑1 (MQ=255)
                                                                                                                                    actggCACCAGTAACGGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATg                     >  1:1174956/5‑145 (MQ=255)
                                                                                                                                    actggCACCAGTAACGGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGagcag                <  2:954385/146‑1 (MQ=255)
                                                                                                                                         cACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTa           <  1:978548/150‑1 (MQ=255)
                                                                                                                                           cccGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATc         >  1:1327696/1‑150 (MQ=255)
                                                                                                                                               ttACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCtt     <  2:335177/150‑1 (MQ=255)
                                                                                                                                               ttACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCtt     <  2:1367162/150‑1 (MQ=255)
                                                                                                                                               ttACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCtt     <  2:341098/150‑1 (MQ=255)
                                                                                                                                                tACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTg    >  1:1628535/1‑150 (MQ=255)
                                                                                                                                                 aCTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGc   <  2:460891/150‑1 (MQ=255)
                                                                                                                                                  gtGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACTCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGAGCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCg  <  2:875286/149‑1 (MQ=255)
                                                                                                                                                  cTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGGTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCg  <  2:875287/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCG  >  NC_000913/928161‑928456

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: ‑