Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,418,628 G→A 8.5% intergenic (+17/‑414) yhdV → / → yhdW putative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,418,6280GA 8.5% 359.3 / 16.4 117intergenic (+17/‑414)yhdV/yhdWputative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (64/43);  new base A (3/7);  total (67/50)
Fisher's exact test for biased strand distribution p-value = 9.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 118 total aligned reads displayed.

GGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTT  >  NC_000913/3418480‑3418773
                                                                                                                                                    |                                                                                                                                                 
ggTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAg                                                                                                                                                   >  1:1087417/1‑149 (MQ=255)
    gtcTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTcc                                                                                                                                               <  4:505517/149‑1 (MQ=255)
      cTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCcgc                                                                                                                                             <  1:1009364/149‑1 (MQ=255)
       tGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCcgcg                                                                                                                                            >  4:566767/1‑149 (MQ=255)
       tGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGGTGcgcg                                                                                                                                            >  3:817025/1‑149 (MQ=255)
         gTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCgcgcc                                                                                                                                          <  2:990293/149‑2 (MQ=255)
          ttCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGctc                                                                                                                                         >  3:851436/1‑149 (MQ=255)
          ttCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGctc                                                                                                                                         >  4:237648/1‑149 (MQ=255)
           tCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                                                                        <  4:657115/149‑1 (MQ=255)
           tCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGAGCGCGGAACTCCGCGccct                                                                                                                                        <  3:657510/149‑1 (MQ=255)
           tCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGAGCGCGGAACTCCGCGccct                                                                                                                                        <  4:102843/149‑1 (MQ=255)
            cacaCCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCtt                                                                                                                                       <  4:771763/149‑1 (MQ=255)
             acacCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttt                                                                                                                                      >  1:1063739/1‑149 (MQ=255)
               acCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGTGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttttt                                                                                                                                    >  1:760989/1‑149 (MQ=255)
               acCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttttt                                                                                                                                    <  1:53941/149‑1 (MQ=255)
                 cAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATCGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCCCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTaa                                                                                                                                  <  4:791793/149‑1 (MQ=255)
                   aTACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATAATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACg                                                                                                                                >  1:191131/1‑149 (MQ=255)
                    tACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGa                                                                                                                               <  1:556563/149‑1 (MQ=255)
                     aCCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGAt                                                                                                                              <  2:1087417/149‑1 (MQ=255)
                      ccTGTATGAAAGGCACGACACAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGGACAGGGCGCGGAGTTCCGCGCTCTTTTTAGCGATa                                                                                                                             <  2:961962/149‑1 (MQ=255)
                        tGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACa                                                                                                                           >  3:565117/1‑149 (MQ=255)
                         gTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAg                                                                                                                          >  4:445437/1‑149 (MQ=255)
                             gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtc                                                                                                                      >  2:809723/1‑149 (MQ=255)
                             gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtc                                                                                                                      >  3:918807/1‑149 (MQ=255)
                             gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtc                                                                                                                      >  3:764017/1‑149 (MQ=255)
                                aGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctct                                                                                                                   >  3:1032830/1‑149 (MQ=255)
                                aGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctct                                                                                                                   <  3:239229/149‑1 (MQ=255)
                                         tCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGAt                                                                                                          >  2:545760/1‑149 (MQ=255)
                                           aGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaa                                                                                                        >  2:681137/1‑149 (MQ=255)
                                             gATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                      <  3:237648/149‑1 (MQ=255)
                                               tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTa                                                                                                                        >  2:1000482/1‑129 (MQ=255)
                                               tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTa                                                                                                                        <  1:1000482/129‑1 (MQ=255)
                                               tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAata                                                                                                    >  4:706373/1‑149 (MQ=255)
                                               tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                                                                        >  1:678940/1‑113 (MQ=255)
                                               tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                                                                        <  2:678940/113‑1 (MQ=255)
                                                  cTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTg                                                                                                 <  4:52729/149‑1 (MQ=255)
                                                   tGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATAtttc                                                                                                >  4:623542/1‑147 (MQ=255)
                                                     ggAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                      <  4:778407/141‑1 (MQ=255)
                                                     ggAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCta                                                                                              >  3:432149/1‑149 (MQ=255)
                                                     ggAAACGGCTGAGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCta                                                                                              >  3:614190/1‑149 (MQ=255)
                                                      gAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                      >  4:302580/1‑140 (MQ=255)
                                                      gAAACGGATGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                      <  3:302580/140‑1 (MQ=255)
                                                         aCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATAtt                                                                                          >  3:57911/1‑149 (MQ=255)
                                                         aCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATAtt                                                                                          >  4:992185/1‑149 (MQ=255)
                                                             cTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctgc                                                                                      <  3:566767/149‑1 (MQ=255)
                                                              tGGGGCTATTGCCGGCGGGGCAGCTACTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGCCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctgct                                                                                     <  3:191239/149‑1 (MQ=255)
                                                                gggCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTgg                                                                                   <  2:1035412/149‑1 (MQ=255)
                                                                gggCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTgg                                                                                   >  2:193310/1‑149 (MQ=255)
                                                                    tATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTaa                                                                               >  3:1023793/1‑149 (MQ=255)
                                                                          cggcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCct                                                                         <  1:1080959/149‑1 (MQ=255)
                                                                            gcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgc                                                                       <  1:746672/149‑1 (MQ=255)
                                                                              ggggCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAGCTGAACAGGGCTCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgctg                                                                     <  4:1034117/149‑1 (MQ=255)
                                                                              ggggCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgctg                                                                     >  1:648634/1‑149 (MQ=255)
                                                                                 gCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAtt                                                                  >  3:1030135/1‑149 (MQ=255)
                                                                                   aGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTa                                                                <  4:1032830/149‑1 (MQ=255)
                                                                                   aGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTa                                                                <  2:366004/149‑1 (MQ=255)
                                                                                         tgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGt                                                          >  3:71864/1‑149 (MQ=255)
                                                                                               tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCAt                                                            <  2:196200/141‑1 (MQ=255)
                                                                                               tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCAt                                                            >  1:196200/1‑141 (MQ=255)
                                                                                               tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATc                                                    >  1:294249/1‑149 (MQ=255)
                                                                                                   cTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGt                                                >  3:995224/1‑149 (MQ=255)
                                                                                                    tGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGtt                                               <  1:681137/149‑1 (MQ=255)
                                                                                                     gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTa                                              <  1:25900/149‑1 (MQ=255)
                                                                                                     gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTa                                              >  1:947793/1‑149 (MQ=255)
                                                                                                      aCGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTAt                                             >  1:883163/1‑149 (MQ=255)
                                                                                                        gATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTAttt                                           >  4:285763/1‑149 (MQ=255)
                                                                                                            gggATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgct                                       >  2:696647/1‑149 (MQ=255)
                                                                                                              gATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgctgc                                     <  1:1051918/149‑1 (MQ=255)
                                                                                                               atcatcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgctgct                                    >  3:678683/1‑149 (MQ=255)
                                                                                                                tcatcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCtt                                   >  2:48168/1‑149 (MQ=255)
                                                                                                                  atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATaa                                                                                                             <  2:315088/73‑1 (MQ=255)
                                                                                                                  atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTcc                                                                           <  1:690160/107‑1 (MQ=255)
                                                                                                                  atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTcc                                                                           >  2:690160/1‑107 (MQ=255)
                                                                                                                  atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTc                                 >  4:840391/1‑149 (MQ=255)
                                                                                                                  atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTc                                 >  3:12444/1‑149 (MQ=255)
                                                                                                                  atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTc                                 <  4:12444/149‑1 (MQ=255)
                                                                                                                  atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTc                                 >  1:736812/1‑149 (MQ=255)
                                                                                                                      ctctATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATGGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGaaa                             <  2:760989/149‑1 (MQ=255)
                                                                                                                       tctATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAAt                            <  3:202276/149‑1 (MQ=255)
                                                                                                                         tATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATaa                          >  1:353246/1‑149 (MQ=255)
                                                                                                                          aTCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAAc                         >  3:790730/1‑149 (MQ=255)
                                                                                                                           tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                        >  3:122637/1‑149 (MQ=255)
                                                                                                                           tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                        >  4:362121/1‑149 (MQ=255)
                                                                                                                           tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                        >  3:642877/1‑149 (MQ=255)
                                                                                                                           tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                        >  3:726837/1‑149 (MQ=255)
                                                                                                                             gAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCa                      <  4:307128/149‑1 (MQ=255)
                                                                                                                              aaaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCAt                     >  4:494095/1‑149 (MQ=255)
                                                                                                                                aTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCa                   >  1:228918/1‑149 (MQ=255)
                                                                                                                                 tAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaa                  <  1:193310/149‑1 (MQ=255)
                                                                                                                                  aaCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaat                 >  2:897438/1‑149 (MQ=255)
                                                                                                                                   aCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatt                <  4:842111/149‑1 (MQ=255)
                                                                                                                                      gAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaa             <  2:769819/149‑1 (MQ=255)
                                                                                                                                      gAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaa             >  3:89640/1‑149 (MQ=255)
                                                                                                                                        aaagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaatt           >  1:518772/6‑149 (MQ=255)
                                                                                                                                         cAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaattt          >  4:952161/1‑149 (MQ=255)
                                                                                                                                          agagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTg                                                                                                 >  3:808397/4‑61 (MQ=38)
                                                                                                                                          aGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAAtttt         >  4:814033/1‑149 (MQ=255)
                                                                                                                                            ggCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGACTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTca       <  3:623542/149‑1 (MQ=255)
                                                                                                                                               gcgGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCacac    >  1:934578/1‑149 (MQ=255)
                                                                                                                                                 ggAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACtt  >  4:556867/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTT  >  NC_000913/3418480‑3418773

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: