Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,445 A→T 8.4% intergenic (‑34/+25) dtpC ← / ← cadA dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,4450AT 8.4% 360.3 / 19.2 109intergenic (‑34/+25)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (55/44);  new base T (6/3);  total (61/48)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 112 total aligned reads displayed.

GATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGG  >  NC_000913/4356300‑4356593
                                                                                                                                                 |                                                                                                                                                    
gATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGc                                                                                                                                                   <  3:175096/149‑1 (MQ=255)
      tGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACtt                                                                                                                                              <  1:322777/148‑1 (MQ=255)
      tGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                                                                              <  2:455622/148‑1 (MQ=255)
         gaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                                                                           >  2:264402/1‑148 (MQ=255)
          aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCt                                                                                                                                          <  2:105145/148‑1 (MQ=255)
          aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCt                                                                                                                                          <  1:528825/148‑1 (MQ=255)
           tAGAGAATGGGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTATCATGTGTTCTCCTTATGAGCAAAAAAGGGATGTGGCAATCCCATTCCCtt                                                                                                                                         >  2:415507/1‑148 (MQ=255)
            agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCGACTTCCCTTGt                                                                                                                                       >  4:320472/1‑149 (MQ=255)
            agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTg                                                                                                                                        >  2:609030/1‑148 (MQ=255)
             gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                                                                                                                       >  2:217413/1‑148 (MQ=255)
                 aTGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGa                                                                                                                                   >  1:852990/1‑148 (MQ=255)
                     gTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTaa                                                                                                                              <  4:322918/149‑1 (MQ=255)
                      tAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTaa                                                                                                                              >  1:918551/1‑148 (MQ=255)
                      tAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAat                                                                                                                             >  4:61861/1‑149 (MQ=255)
                         ggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatta                                                                                                                           >  2:447941/1‑148 (MQ=255)
                             cGCTTGCCGTAAAAACTGAAGTACCCGTAGATTTGGAGCGCCGCGATATAGTAGATCGCGCGCTGCTGTGAGGGTGTTTTAATGTGTTCTCCTTGTGAGCAAAAAAGGTAAGTGGCGAGCTTCTTCCGTTGTACGAGCTAATTAtttt                                                                                                                       <  2:745819/148‑1 (MQ=255)
                              gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                                                                     <  3:320472/149‑1 (MQ=255)
                              gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttttt                                                                                                                      <  1:689124/148‑1 (MQ=255)
                                aTGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGc                                                                                                                   >  4:413653/1‑149 (MQ=255)
                                 tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATACAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGc                                                                                                                   >  1:441048/1‑148 (MQ=255)
                                      tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttctt                                                                                                             >  3:201592/1‑149 (MQ=255)
                                          aaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTcttctt                                                                                                          <  1:34823/148‑1 (MQ=255)
                                            cTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                                                                        >  1:29404/1‑148 (MQ=255)
                                             tGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                                                                      <  3:317686/149‑1 (MQ=255)
                                              gAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                                                                      >  2:581119/1‑148 (MQ=255)
                                              gAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                                                                      >  2:620997/1‑148 (MQ=255)
                                               aaGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAAt                                                                                                     >  2:656663/1‑148 (MQ=255)
                                                aGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAc                                                                                                   <  4:876293/149‑1 (MQ=255)
                                                aGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAc                                                                                                   >  4:462818/1‑149 (MQ=255)
                                                 gTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                                                                  >  4:217513/1‑149 (MQ=255)
                                                      cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                                                                              <  2:796395/148‑1 (MQ=255)
                                                       ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatt                                                                                                                            <  3:423568/117‑1 (MQ=255)
                                                       ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatt                                                                                                                            >  4:423568/1‑117 (MQ=255)
                                                       ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                           >  4:691038/1‑134 (MQ=255)
                                                       ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                           <  3:691038/134‑1 (MQ=255)
                                                        cAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                            >  2:697645/1‑148 (MQ=255)
                                                        cAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                            >  1:235496/1‑148 (MQ=255)
                                                               ggATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCg                                                                                    <  3:462818/149‑1 (MQ=255)
                                                                 aTCGCCACGATATAGTATATCGCACGCGGCGGTGGGGGTGTTTTCATGTGTTCTCCTTATGATCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTATTCTTTCAATACCTTAACGGTATAGCGGc                                                                                  >  4:535405/1‑149 (MQ=255)
                                                                  tCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTAAGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                                                                 <  4:810270/149‑1 (MQ=255)
                                                                      cacaATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATc                                                                              <  1:415507/148‑1 (MQ=255)
                                                                      cACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATc                                                                              >  1:873623/1‑148 (MQ=255)
                                                                      cACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATc                                                                              <  1:581119/148‑1 (MQ=255)
                                                                       aCGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                                                                             >  1:796829/1‑148 (MQ=255)
                                                                             taGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGAc                                                                      >  3:32349/1‑149 (MQ=255)
                                                                                tataTCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                                                                   >  3:184998/1‑149 (MQ=255)
                                                                                 atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGtat                                                                                        <  4:242269/127‑1 (MQ=255)
                                                                                 atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGtat                                                                                        >  3:242269/1‑127 (MQ=255)
                                                                                 atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                                                                   <  1:620997/148‑1 (MQ=255)
                                                                                 atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                                                                   >  2:226182/1‑148 (MQ=255)
                                                                                   atCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg                                                                >  3:909667/1‑149 (MQ=255)
                                                                                           cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGa                                                        <  3:32913/149‑1 (MQ=255)
                                                                                                   aGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGt                                                <  4:801804/149‑1 (MQ=255)
                                                                                                   aGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGt                                                <  3:950685/149‑1 (MQ=255)
                                                                                                    gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGtt                                               <  3:357598/149‑1 (MQ=255)
                                                                                                    gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGtt                                               <  4:761755/149‑1 (MQ=255)
                                                                                                       tgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGa                                                        >  1:419106/1‑137 (MQ=255)
                                                                                                       tgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGa                                                        <  2:419106/137‑1 (MQ=255)
                                                                                                             cATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGcc                                       <  2:298634/148‑1 (MQ=255)
                                                                                                              aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCg                                     >  3:203238/1‑149 (MQ=255)
                                                                                                              aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCg                                     >  4:263247/1‑149 (MQ=255)
                                                                                                              aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCg                                     >  4:275220/1‑149 (MQ=255)
                                                                                                              aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCg                                     >  4:399771/1‑149 (MQ=255)
                                                                                                                gtgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCgg                                    <  2:918551/148‑1 (MQ=255)
                                                                                                                gtgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCgg                                    <  2:319098/148‑1 (MQ=255)
                                                                                                                   ttCTCCTTATGAGCAAAAAAGGGAAGTGCCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTACGCACCGTGAATATCGGTTTCAAAGCCCGGATAg                                >  4:391851/1‑149 (MQ=255)
                                                                                                                    tctcCTTATGAGCTAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTGCTTCTTTCAATACCTTAACCGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGt                               <  4:606791/149‑1 (MQ=255)
                                                                                                                    tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                           <  1:958804/73‑1 (MQ=255)
                                                                                                                    tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                           >  2:958804/1‑73 (MQ=255)
                                                                                                                    tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGt                               <  3:755957/149‑1 (MQ=255)
                                                                                                                       ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGa                             <  2:441048/148‑1 (MQ=255)
                                                                                                                         ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAgcg                          >  3:693634/1‑149 (MQ=255)
                                                                                                                          tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATACTGAgcg                          >  1:259266/1‑148 (MQ=255)
                                                                                                                             gAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCg                       >  1:839078/1‑148 (MQ=255)
                                                                                                                              aGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAg                            >  2:506341/1‑142 (MQ=255)
                                                                                                                              aGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGa                      <  2:796829/148‑1 (MQ=255)
                                                                                                                              aGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCAGGATAGTGAg                            <  1:506341/142‑1 (MQ=255)
                                                                                                                               gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGAtt                    <  3:61861/149‑1 (MQ=255)
                                                                                                                                 aaaaaaGGGAAGTGGCACGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTc                  >  3:92502/1‑149 (MQ=255)
                                                                                                                                  aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCg                                                           >  2:838592/1‑107 (MQ=255)
                                                                                                                                  aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCg                                                           <  1:838592/107‑1 (MQ=255)
                                                                                                                                    aaaGGTAAGTGGCAAGCCAGTTCCCTTGTACGAGATAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATAGGTTTCAAAGCCCGGATAGGGAGCGCCGATTTcaca               >  4:538459/1‑149 (MQ=255)
                                                                                                                                     aaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcaca               >  2:29310/1‑148 (MQ=255)
                                                                                                                                      aGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcacaca             >  3:117361/1‑149 (MQ=255)
                                                                                                                                      aGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcacac              >  2:379802/1‑148 (MQ=255)
                                                                                                                                      aGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGaa                                                       >  3:495882/1‑107 (MQ=255)
                                                                                                                                      aGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGaa                                                       <  4:495882/107‑1 (MQ=255)
                                                                                                                                       gggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATTCCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcacaca             >  2:569870/1‑148 (MQ=255)
                                                                                                                                       gggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAg                                >  1:391355/1‑129 (MQ=255)
                                                                                                                                       gggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAg                                <  2:391355/129‑1 (MQ=255)
                                                                                                                                        ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGc           >  4:329469/1‑149 (MQ=255)
                                                                                                                                         gAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCa          >  3:197340/1‑149 (MQ=255)
                                                                                                                                           aGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCAt         <  1:880111/148‑1 (MQ=255)
                                                                                                                                            gTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATTCCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATc        <  1:569870/148‑1 (MQ=255)
                                                                                                                                               gCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGc     >  1:811301/1‑148 (MQ=255)
                                                                                                                                               gCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCa    <  4:184998/149‑1 (MQ=255)
                                                                                                                                                cAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCa    <  2:266013/148‑1 (MQ=255)
                                                                                                                                                 aaGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCa    <  4:79370/147‑1 (MQ=255)
                                                                                                                                                 aaGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCa    >  3:79370/1‑147 (MQ=255)
                                                                                                                                                 aaGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAgg  <  4:257109/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGG  >  NC_000913/4356300‑4356593

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: