Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,676,196 A→C 5.2% intergenic (+88/+94) yhjE → / ← yhjG putative MFS transporter; membrane protein/putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,676,1960AC 5.2% 407.8 / 11.2 116intergenic (+88/+94)yhjE/yhjGputative MFS transporter; membrane protein/putative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base A (64/46);  new base C (3/3);  total (67/49)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 118 total aligned reads displayed.

GGCTGGCTTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAA  >  NC_000913/3676048‑3676344
                                                                                                                                                    |                                                                                                                                                    
ggctggctTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATaa                                                                                                                                                      <  3:822923/149‑1 (MQ=255)
 gctggctAGACGTTAATCGCACTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTCAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAg                                                                                                                                                     <  1:690047/149‑1 (MQ=255)
    ggctTGACGTTAATCGCCCTGATGCTGACCCATGAGCCGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCTTAAGGATTTGTAGGCCGGATAAGTCg                                                                                                                                                  <  2:198015/149‑1 (MQ=39)
     gctTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGt                                                                                                                                                 <  1:698536/149‑1 (MQ=255)
      ctTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAAGCCGGATAAGGCGtt                                                                                                                                                >  1:463519/1‑149 (MQ=255)
      cctGACGTTAATCGCCCTGCTGCTGACCCAAGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGtt                                                                                                                                                <  3:718467/147‑1 (MQ=255)
           cGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                           <  1:71042/149‑1 (MQ=255)
            gTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                          <  4:583729/149‑1 (MQ=255)
            gTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                          >  2:507900/1‑149 (MQ=255)
               aaTCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATGCGGCGTGAAcgcctta                                                                                                                                       >  3:128401/1‑146 (MQ=39)
               aaTCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                                                       >  4:1601/1‑149 (MQ=255)
                aTCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                      <  4:90905/149‑5 (MQ=39)
                aTCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                      <  4:128401/149‑5 (MQ=39)
                 tCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                     >  2:765124/1‑149 (MQ=255)
                 tCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                     >  3:837558/1‑149 (MQ=255)
                       tgctgctgACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                               <  3:81920/149‑1 (MQ=255)
                       tgctgctgACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                               <  2:842516/149‑1 (MQ=255)
                        gctgctgACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                              <  2:298585/149‑1 (MQ=255)
                          tgctgACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                                                                                                                            <  4:477263/149‑1 (MQ=255)
                               aCCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAg                                                                                                                       <  2:744571/149‑1 (MQ=255)
                                 ccATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGca                                                                                                                     >  3:697708/1‑149 (MQ=255)
                                 ccATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGca                                                                                                                     >  1:123638/1‑149 (MQ=255)
                                 ccATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGca                                                                                                                     >  4:68589/1‑149 (MQ=255)
                                  cATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGcac                                                                                                                    >  1:123564/1‑149 (MQ=255)
                                  cATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGcac                                                                                                                    >  1:96207/1‑149 (MQ=255)
                                  cATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGcac                                                                                                                    >  1:484364/1‑149 (MQ=255)
                                    tGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCacac                                                                                                                  >  1:294255/1‑149 (MQ=255)
                                     gagaCGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACt                                                                                                                 <  4:131561/149‑1 (MQ=255)
                                         cgcgACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGc                                                                                                             <  2:432380/149‑1 (MQ=255)
                                             aCATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaac                                                                                                         >  1:463841/1‑149 (MQ=255)
                                             aCATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaac                                                                                                         >  1:105910/1‑149 (MQ=255)
                                                  aGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCt                                                                                                    >  4:610680/1‑149 (MQ=255)
                                                   gTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTa                                                                                                   >  4:525350/1‑149 (MQ=255)
                                                   gTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTa                                                                                                   >  2:49524/1‑149 (MQ=255)
                                                   gTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTa                                                                                                   >  2:674727/1‑149 (MQ=255)
                                                      gTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAAt                                                                                                >  1:285788/1‑149 (MQ=255)
                                                      gTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAAt                                                                                                >  2:74522/1‑149 (MQ=255)
                                                      gTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAAt                                                                                                >  1:145109/1‑149 (MQ=255)
                                                      gTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAAt                                                                                                <  1:674727/149‑1 (MQ=255)
                                                          taataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATGTTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTCTGCCAACAATCTAAATCTGt                                                                                            <  1:713916/149‑1 (MQ=255)
                                                          taataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGt                                                                                            >  2:542144/1‑149 (MQ=255)
                                                          taataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGt                                                                                            >  4:87440/1‑149 (MQ=255)
                                                            ataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGttt                                                                                          >  4:173087/1‑149 (MQ=255)
                                                             taaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTa                                                                                         >  3:54364/1‑149 (MQ=255)
                                                                gagTGAAAAGTCTGATGCGCTACGCTTAGCAGGCCGGAATTTACATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGACAACAATCTAAATCTGTTTAAtt                                                                                      <  2:305146/149‑1 (MQ=25)
                                                                gagTGAAAAGCTTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAAtt                                                                                      >  2:843022/1‑149 (MQ=255)
                                                                gagTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAAtt                                                                                      >  2:773787/1‑149 (MQ=255)
                                                                   tGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGcc                                                                                   <  2:119586/149‑1 (MQ=255)
                                                                          ccTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTc                                                                            <  2:484364/149‑1 (MQ=255)
                                                                          ccTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTc                                                                            <  2:96207/149‑1 (MQ=255)
                                                                           cTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCa                                                                           <  2:46475/149‑1 (MQ=255)
                                                                            tGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAg                                                                          >  3:256809/1‑149 (MQ=255)
                                                                             gATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGa                                                                         >  3:601920/1‑149 (MQ=255)
                                                                                  gcTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATcc                                                                    >  2:395928/1‑149 (MQ=255)
                                                                                  gcTACGCTTATCAGGCCTGAATTTCCATACAATATAGTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATcc                                                                    >  4:835170/1‑149 (MQ=255)
                                                                                      cGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCa                                                                <  3:1601/149‑1 (MQ=255)
                                                                                       gCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCaa                                                               <  4:440498/149‑1 (MQ=255)
                                                                                       gCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCaa                                                               <  4:117088/149‑1 (MQ=255)
                                                                                         ttATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAt                                                              >  1:383106/1‑148 (MQ=255)
                                                                                         ttATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAt                                                              <  2:383106/148‑1 (MQ=255)
                                                                                           aTCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTat                                                           >  3:854350/1‑149 (MQ=255)
                                                                                            tCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatc                                                          <  2:105910/149‑1 (MQ=255)
                                                                                              aGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatcat                                                        <  2:830790/149‑1 (MQ=255)
                                                                                              aGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatcat                                                        <  3:458188/149‑1 (MQ=255)
                                                                                                gCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTAtcatca                                                      >  4:639579/1‑149 (MQ=255)
                                                                                                gCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTGGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCACTTAtaatca                                                      >  1:782315/1‑149 (MQ=255)
                                                                                                    gAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCActtc                                                  >  4:676429/1‑149 (MQ=255)
                                                                                                      aTTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCActtctt                                                <  3:68589/149‑1 (MQ=255)
                                                                                                       tttCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACttcttc                                               <  3:173087/149‑1 (MQ=255)
                                                                                                            aTACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGAt                                                                               <  1:495214/112‑1 (MQ=255)
                                                                                                            aTACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGAt                                                                               >  2:495214/1‑112 (MQ=255)
                                                                                                             tACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGAt                                                                               <  2:403287/111‑1 (MQ=255)
                                                                                                             tACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGAt                                                                               >  1:403287/1‑111 (MQ=255)
                                                                                                             tACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGc                                         >  4:218589/1‑149 (MQ=255)
                                                                                                             tACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGc                                         >  3:660991/1‑149 (MQ=255)
                                                                                                             tACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGc                                         <  3:610680/149‑1 (MQ=255)
                                                                                                             tACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGc                                         >  4:443946/1‑149 (MQ=255)
                                                                                                                 atatatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACa                                     >  2:777030/1‑149 (MQ=255)
                                                                                                                    tatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACttcttc                                               >  1:358142/1‑136 (MQ=255)
                                                                                                                    tatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACaaaa                                  <  2:276160/149‑1 (MQ=255)
                                                                                                                      ttaaaTTCGCAAGGATTTCTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGTTGTCAGACAGCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATc                                >  4:69763/4‑149 (MQ=255)
                                                                                                                      ttGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATc                                >  1:300717/1‑149 (MQ=255)
                                                                                                                        gAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGt                              <  3:639579/149‑1 (MQ=255)
                                                                                                                         aaTTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTc                             <  2:729706/149‑1 (MQ=255)
                                                                                                                          aTTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTcc                            >  4:284124/1‑149 (MQ=255)
                                                                                                                               cAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGacac                       >  4:750085/1‑149 (MQ=255)
                                                                                                                                aaGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACt                      >  2:413537/1‑149 (MQ=255)
                                                                                                                                  ggATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCAGCGACAAAATCGTCCGACACTGa                    >  1:198025/1‑149 (MQ=255)
                                                                                                                                   gATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGAt                   >  1:164024/1‑149 (MQ=255)
                                                                                                                                   gATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGAt                   <  2:123638/149‑1 (MQ=255)
                                                                                                                                     tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTa                 >  2:628754/1‑149 (MQ=255)
                                                                                                                                      ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAg                <  1:88326/149‑1 (MQ=255)
                                                                                                                                          atgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTc            >  1:646413/3‑149 (MQ=39)
                                                                                                                                           tgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTcc           >  1:575224/2‑149 (MQ=39)
                                                                                                                                               ggATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCtt       >  4:56755/1‑149 (MQ=255)
                                                                                                                                               ggATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATAACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCtt       >  4:313125/1‑149 (MQ=255)
                                                                                                                                                gATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTc      >  3:156406/1‑149 (MQ=255)
                                                                                                                                                 aTAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCg     >  4:202559/1‑149 (MQ=255)
                                                                                                                                                  tAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCgg    >  2:225797/1‑149 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCAGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGaa  <  4:601920/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCTGGCTTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAA  >  NC_000913/3676048‑3676344

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