Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,216,652 G→A 5.9% intergenic (‑161/‑127) yqjH ← / → yqjI putative siderophore interacting protein/PadR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,216,6520GA 5.9% 455.4 / 11.3 154intergenic (‑161/‑127)yqjH/yqjIputative siderophore interacting protein/PadR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (58/83);  new base A (3/6);  total (62/92)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-02

Only 100 of 166 total aligned reads displayed.

CGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTG  >  NC_000913/3216504‑3216801
                                                                                                                                                    |                                                                                                                                                     
cGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgc                                                                                                                                                      <  1:170348/150‑1 (MQ=255)
 gTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcg                                                                                                                                                     >  2:920485/1‑150 (MQ=255)
 gTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcg                                                                                                                                                     >  1:1586520/1‑150 (MQ=255)
  tGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTTTAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcgg                                                                                                                                                    >  2:694800/1‑150 (MQ=255)
  tGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcgg                                                                                                                                                    <  1:737133/150‑1 (MQ=255)
   gATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTTCTCAACGCCGGATgcggc                                                                                                                                                   <  1:1546326/150‑1 (MQ=255)
   gATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcggc                                                                                                                                                   <  2:1104872/150‑1 (MQ=255)
    aTGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTCTAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcggcg                                                                                                                                                  <  1:310461/150‑1 (MQ=255)
     tGGTAATCATATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGt                                                                                                                                                 <  1:382057/150‑1 (MQ=255)
      ggTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTg                                                                                                                                                >  1:1442515/1‑150 (MQ=255)
      ggTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTg                                                                                                                                                <  2:927888/150‑1 (MQ=255)
       gTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGa                                                                                                                                               >  1:1526231/1‑150 (MQ=255)
        tAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGaa                                                                                                                                              >  1:254123/1‑150 (MQ=255)
         aaTCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAAc                                                                                                                                             >  2:1123808/1‑150 (MQ=255)
         aaTCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGCATCCGGCGttgag                                                                                                                                             <  1:950157/150‑5 (MQ=255)
           tCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACAGATATGTCTATAGTCTGATAATACGCACCGCCTCTTCTCAGTCATCATTACTCAACGCCGGATGCGGCGTGAActc                                                                                                                                           >  1:1158132/1‑148 (MQ=255)
             aGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCt                                                                                                                                         >  2:1330390/1‑150 (MQ=255)
                  tatCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    <  2:144073/150‑1 (MQ=255)
                     cTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGCACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                 <  1:1392800/150‑1 (MQ=255)
                     cTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                 <  1:1393462/150‑1 (MQ=255)
                     cTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGGCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                 <  1:1519985/150‑1 (MQ=255)
                      aaaataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                <  2:814016/149‑1 (MQ=255)
                           aaaaCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGcgttgag                                                                                                                           <  2:889879/150‑7 (MQ=38)
                           aaaaCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATAAGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGcgttgag                                                                                                                           <  2:602349/150‑7 (MQ=38)
                            aaaCTCGCAAATTATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATACGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGcgttgagt                                                                                                                          <  2:1423565/150‑8 (MQ=38)
                            aaaCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGcgttgagt                                                                                                                          <  2:595018/150‑8 (MQ=38)
                              aCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACACCGGATGCGGCGTGAACGCCTTATCCGGCCTACGtgtg                                                                                                                        <  1:1725601/150‑1 (MQ=255)
                                tCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTgag                                                                                                                      <  1:1248644/150‑1 (MQ=255)
                                 cGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTgaga                                                                                                                     >  2:868534/1‑150 (MQ=255)
                                  gCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGAt                                                                                                                    >  1:1598085/1‑150 (MQ=255)
                                      atgataatgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGt                                                                                                                <  2:1600866/150‑1 (MQ=255)
                                          taatgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCAc                                                                                                            <  2:631960/150‑1 (MQ=255)
                                          taatgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCAc                                                                                                            >  1:308617/1‑150 (MQ=255)
                                            atgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGt                                                                                                          <  2:885601/150‑1 (MQ=255)
                                                 tGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTGTaaa                                                                                                     >  2:1626685/1‑150 (MQ=255)
                                                  gTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaaaa                                                                                                    >  2:1445736/1‑150 (MQ=255)
                                                  gTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTGTaaaa                                                                                                    >  1:1729466/1‑150 (MQ=255)
                                                   ttAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCGCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaaaaa                                                                                                   <  2:1279829/150‑1 (MQ=255)
                                                    tAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTaaa                                                                                                     >  2:821365/1‑147 (MQ=255)
                                                      aTCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCa                                                                                                <  1:815085/150‑1 (MQ=255)
                                                       tCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAt                                                                                               >  2:899403/1‑150 (MQ=255)
                                                        cATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAtt                                                                                              >  2:117570/1‑150 (MQ=255)
                                                         aTGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAttt                                                                                             <  2:256287/150‑1 (MQ=255)
                                                           gATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAttttt                                                                                           >  2:38336/1‑150 (MQ=255)
                                                           gATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAttttt                                                                                           <  2:733628/150‑1 (MQ=255)
                                                                aTGCAAGGGATTTGTCGATCTGATATGTCTATAGTCTGCTAAGACGAACCGCCTCTTCTCAGGCATGATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACt                                                                                      <  2:865431/150‑1 (MQ=39)
                                                                  gCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTg                                                                                    <  2:1652174/150‑1 (MQ=255)
                                                                   cAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGc                                                                                   <  2:1507375/150‑1 (MQ=255)
                                                                     aGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATc                                                                                                 <  2:1575868/134‑1 (MQ=255)
                                                                     aGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCaa                                                                                 >  1:1669348/1‑150 (MQ=255)
                                                                      gCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAg                                                                                >  2:1168623/1‑150 (MQ=255)
                                                                            tGTAGAACTTATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAATAATCATTTTTACACTTGCAAGAACGCt                                                                          <  2:890538/150‑1 (MQ=255)
                                                                             gTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTc                                                                         <  2:1199625/150‑1 (MQ=255)
                                                                              tAGAACTTATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCa                                                                        <  1:572753/150‑1 (MQ=255)
                                                                               aGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCat                                                                       <  1:889797/150‑1 (MQ=255)
                                                                                 aaCTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCatat                                                                     >  1:855343/1‑150 (MQ=255)
                                                                                  aCTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGCCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCatata                                                                    >  2:850377/1‑149 (MQ=255)
                                                                                    tGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACATGCAAGAACGCTCATATCAg                                                                  <  1:1007723/150‑1 (MQ=255)
                                                                                     gATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGa                                                                 >  1:855201/1‑150 (MQ=255)
                                                                                      atatGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGat                                                                <  1:1333306/150‑1 (MQ=255)
                                                                                       tatGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGata                                                               <  1:1123808/150‑1 (MQ=255)
                                                                                        atGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGTCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGatat                                                              <  1:1494399/150‑1 (MQ=255)
                                                                                         tGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGAtata                                                             >  1:577358/1‑150 (MQ=255)
                                                                                          gTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGatatta                                                            >  1:1594007/1‑148 (MQ=255)
                                                                                           tCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATaaa                                                           >  1:277103/1‑150 (MQ=255)
                                                                                               taGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAg                                                       <  1:751749/150‑1 (MQ=255)
                                                                                                aGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGa                                                      <  2:1220934/150‑1 (MQ=255)
                                                                                                 gTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGat                                                     >  1:1675849/1‑150 (MQ=255)
                                                                                                 gTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGat                                                     <  1:920485/150‑1 (MQ=255)
                                                                                                 gTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGat                                                     >  1:1093427/1‑150 (MQ=255)
                                                                                                 gTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCAGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTAATATCAGATATAAATTAGat                                                     >  2:1540518/1‑150 (MQ=255)
                                                                                                  tCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGata                                                    <  2:175551/150‑1 (MQ=255)
                                                                                                  tCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGata                                                    <  2:429927/150‑1 (MQ=255)
                                                                                                    tGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGatata                                                  >  2:980438/1‑150 (MQ=255)
                                                                                                      aTAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATc                                                <  1:1535051/150‑1 (MQ=255)
                                                                                                      aTAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATc                                                <  2:721193/150‑1 (MQ=255)
                                                                                                            cGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTa                                          <  2:1586520/150‑1 (MQ=255)
                                                                                                             gAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTaa                                         >  1:718666/1‑150 (MQ=255)
                                                                                                                 cGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCaa                                     >  2:1127237/1‑150 (MQ=255)
                                                                                                                   ccTCTTCTCAGGCATCATTACTCACCGCCGGATGCGGCGTGAACGCCTTATACGGCCTCCGTGTGAGATGAGTCCACTTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCaaaa                                   <  2:319870/150‑1 (MQ=255)
                                                                                                                       ttctCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCATATATAAATTAGATATATCTAATTAAGCAAAaggag                               <  1:1499064/150‑1 (MQ=255)
                                                                                                                        tctcAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAaggagg                              <  2:1093427/150‑1 (MQ=255)
                                                                                                                        tctcAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAaggagg                              <  1:5554/150‑1 (MQ=255)
                                                                                                                             ggCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGat                         >  2:451679/1‑150 (MQ=255)
                                                                                                                                 tcatTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGa                     >  1:1277052/1‑150 (MQ=255)
                                                                                                                                   atTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGc                   <  2:1038947/150‑1 (MQ=255)
                                                                                                                                   atTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACAATTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGc                   <  1:389410/150‑1 (MQ=255)
                                                                                                                                    ttACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGcc                  <  2:423203/150‑1 (MQ=255)
                                                                                                                                       cTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatc               >  1:1011764/1‑150 (MQ=255)
                                                                                                                                       cTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatc               >  2:739514/1‑150 (MQ=255)
                                                                                                                                         cAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatcat             >  2:665929/1‑150 (MQ=255)
                                                                                                                                          aaCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAACTCATTTTTACACATGCAAGCACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGGTGAGATGCGCCatcatc            >  1:1733787/1‑150 (MQ=39)
                                                                                                                                          aaCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCatcatc            <  2:180695/150‑1 (MQ=255)
                                                                                                                                           aCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCAtcatca           >  1:516776/1‑150 (MQ=255)
                                                                                                                                            cGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCAc          >  1:250256/1‑150 (MQ=255)
                                                                                                                                             gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAg                    <  2:287609/139‑1 (MQ=255)
                                                                                                                                             accGGATAAGGCGTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACg         >  1:422964/2‑150 (MQ=255)
                                                                                                                                               cGGATAAGGCGTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGaa       >  1:470894/1‑150 (MQ=255)
                                                                                                                                                ggATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAg      <  1:451679/150‑1 (MQ=255)
                                                                                                                                                ggATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAg      <  2:663693/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTG  >  NC_000913/3216504‑3216801

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