Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,216,667 T→C 5.9% intergenic (‑176/‑112) 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,6670TC 5.9% 456.2 / 13.5 151intergenic (‑176/‑112)yqjH/yqjIputative siderophore interacting protein/PadR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (55/87);  new base C (3/6);  total (58/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

Only 100 of 154 total aligned reads displayed.

GATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAACATGAAGG  >  NC_000913/3216518‑3216816
                                                                                                                                                     |                                                                                                                                                     
gATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                       >  1:1507375/1‑150 (MQ=255)
    tatCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                   <  2:144073/150‑1 (MQ=255)
       cTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                <  1:1393462/150‑1 (MQ=255)
       cTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGGCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                <  1:1519985/150‑1 (MQ=255)
        taaataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                               >  1:561349/1‑150 (MQ=255)
        taaataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGAcc                                                                                                                                               <  1:1336248/150‑1 (MQ=255)
          aataaaACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                             >  1:1438507/1‑150 (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)
                    gCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGAt                                                                                                                                   >  1:1598085/1‑150 (MQ=255)
                         tgataatgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTc                                                                                                                              >  2:1349233/1‑150 (MQ=255)
                          gataatgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCATCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTcc                                                                                                                             <  2:841904/150‑1 (MQ=255)
                           ataatgatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCa                                                                                                                            <  2:1135991/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)
                                gatTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGttt                                                                                                                       <  2:639517/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)
                                           aTGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAttt                                                                                                            <  2:256287/150‑1 (MQ=255)
                                             gATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAttttt                                                                                                          >  2:38336/1‑150 (MQ=255)
                                               tAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTac                                                                                                        <  2:714246/150‑1 (MQ=255)
                                                  aTGCAAGGGATTTGTCGATCTGATATGTCTATAGTCTGCTAAGACGAACCGCCTCTTCTCAGGCATGATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACt                                                                                                     <  2:865431/150‑1 (MQ=39)
                                                    gCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTg                                                                                                   <  1:219835/150‑1 (MQ=255)
                                                     cAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGc                                                                                                  >  2:379760/1‑150 (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)
                                                              tGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCt                                                                                         >  1:1221773/1‑150 (MQ=255)
                                                               gTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTc                                                                                        <  2:1199625/150‑1 (MQ=255)
                                                                tAGAACTTATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCa                                                                                       <  1:572753/150‑1 (MQ=255)
                                                                tAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCa                                                                                       <  2:1669348/150‑1 (MQ=255)
                                                                tAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCa                                                                                       <  2:879507/150‑1 (MQ=255)
                                                                 aGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCat                                                                                      <  1:889797/150‑1 (MQ=255)
                                                                   aaCTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCatat                                                                                    >  1:855343/1‑150 (MQ=255)
                                                                    aCTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCAGGCCTACGTGTGAGATGAGTCAACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATc                                                                                   >  2:1039080/1‑150 (MQ=255)
                                                                      tGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACATGCAAGAACGCTCATATCAg                                                                                 <  1:1007723/150‑1 (MQ=255)
                                                                       gATATGTCTATAGTCTGATAAGTCGAACCGCCTCTGCTCAGGCATCATTACTCAACGCCGGATGCTGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTATTACACTTGCAAGAACGCTCATATCAGa                                                                                <  1:267006/150‑1 (MQ=255)
                                                                       gATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGa                                                                                >  1:855201/1‑150 (MQ=255)
                                                                        atatGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGat                                                                               <  2:254123/150‑1 (MQ=255)
                                                                         tatGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGata                                                                              <  1:1123808/150‑1 (MQ=255)
                                                                         tatGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGata                                                                              >  1:446908/1‑150 (MQ=255)
                                                                          atGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGTCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGatat                                                                             <  1:1494399/150‑1 (MQ=255)
                                                                           tGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGAtata                                                                            <  2:1120766/150‑1 (MQ=255)
                                                                            gTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGatatta                                                                           >  1:1594007/1‑148 (MQ=255)
                                                                                ataGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACTTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTa                                                                       <  1:1541616/150‑1 (MQ=255)
                                                                                 taGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAg                                                                      <  1:751749/150‑1 (MQ=255)
                                                                                 taGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAg                                                                      <  2:1644706/150‑1 (MQ=255)
                                                                                  aGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGa                                                                     <  2:1220934/150‑1 (MQ=255)
                                                                                   gTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGat                                                                    >  1:1093427/1‑150 (MQ=255)
                                                                                   gTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGat                                                                    >  1:1675849/1‑150 (MQ=255)
                                                                                    tCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGata                                                                   <  2:429927/150‑1 (MQ=255)
                                                                                    tCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGata                                                                   <  2:175551/150‑1 (MQ=255)
                                                                                      tGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGatata                                                                 >  2:980438/1‑150 (MQ=255)
                                                                                        aTAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATc                                                               <  2:1173743/150‑1 (MQ=255)
                                                                                        aTAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATc                                                               <  1:1535051/150‑1 (MQ=255)
                                                                                        aTAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATc                                                               <  2:721193/150‑1 (MQ=255)
                                                                                            gACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAAt                                                           >  1:533193/1‑150 (MQ=255)
                                                                                              cGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTa                                                         <  2:1586520/150‑1 (MQ=255)
                                                                                               gAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTaa                                                        >  1:718666/1‑150 (MQ=255)
                                                                                                aaCCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGACTTATACGGCCTACTTGTGAGATGAGTCCACGGTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAg                                                       <  2:738278/150‑1 (MQ=255)
                                                                                                    gCCTCTTCTCAGGCATCATTACTCAACGCCGGATTCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCaaa                                                   >  2:903692/1‑150 (MQ=255)
                                                                                                          tctcAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAaggagg                                             <  2:1093427/150‑1 (MQ=255)
                                                                                                               ggCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGat                                        >  2:451679/1‑150 (MQ=255)
                                                                                                                   tcatTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGa                                    >  1:1277052/1‑150 (MQ=255)
                                                                                                                    catTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAg                                   <  1:1330390/150‑1 (MQ=255)
                                                                                                                     atTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGc                                  <  2:1038947/150‑1 (MQ=255)
                                                                                                                     atTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACAATTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGc                                  <  1:389410/150‑1 (MQ=255)
                                                                                                                         cTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatc                              >  1:1011764/1‑150 (MQ=255)
                                                                                                                         cTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatc                              >  2:739514/1‑150 (MQ=255)
                                                                                                                         cTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAATTTTACCCTTGCAAGAACGCTCATAACAGATATAAATTAGATATATCAAATTAAGCAAAAGGAGGCTGATATGAGCcatc                              >  1:908341/1‑150 (MQ=255)
                                                                                                                           cAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCcatcat                            >  2:665929/1‑150 (MQ=255)
                                                                                                                            aaCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAACTCATTTTTACACATGCAAGCACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGGTGAGATGCGCCatcatc                           >  1:1733787/1‑150 (MQ=39)
                                                                                                                             aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCAtcatca                          >  1:652619/3‑150 (MQ=255)
                                                                                                                             aCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCAtcatca                          >  1:516776/1‑150 (MQ=255)
                                                                                                                               gCCGGATTCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACg                        <  2:1596073/150‑1 (MQ=255)
                                                                                                                               gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAg                                   >  1:287609/1‑139 (MQ=255)
                                                                                                                               gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAg                                   <  2:287609/139‑1 (MQ=255)
                                                                                                                                  ggATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAg                     <  2:663693/150‑1 (MQ=255)
                                                                                                                                  ggATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAg                     <  1:451679/150‑1 (MQ=255)
                                                                                                                                   gATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAgg                    <  2:855201/150‑1 (MQ=255)
                                                                                                                                     tGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGgt                  <  1:38336/150‑1 (MQ=255)
                                                                                                                                      gcggcgTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATCTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGgtg                 <  2:1526231/150‑1 (MQ=255)
                                                                                                                                       cggcgTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGgtgt                >  1:813008/1‑150 (MQ=255)
                                                                                                                                       cggcgTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGACATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGgtgt                >  1:1150241/1‑150 (MQ=255)
                                                                                                                                         gcgTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGtgttg              <  2:813008/150‑1 (MQ=255)
                                                                                                                                              aaCGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAAc         <  2:308617/150‑1 (MQ=255)
                                                                                                                                                 gCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAACATg      <  1:1362028/150‑1 (MQ=255)
                                                                                                                                                     tATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAACATGAAgg  <  2:446908/150‑1 (MQ=255)
                                                                                                                                                     tATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAACATGAAgg  >  2:1186131/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAACGCTCATATCAGATATAAATTAGATATATCTAATTAAGCAAAAGGAGGCTGATATGAGCCATCATCACGAAGGGTGTTGTAAACATGAAGG  >  NC_000913/3216518‑3216816

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: