Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,635 T→G 5.1% intergenic (+33/‑117) ybgE → / → ybgC putative inner membrane protein in cydABX‑ybgE operon/acyl‑CoA thioester hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,6350TG 5.1% 323.8 / 4.8 110intergenic (+33/‑117)ybgE/ybgCputative inner membrane protein in cydABX‑ybgE operon/acyl‑CoA thioester hydrolase
Reads supporting (aligned to +/- strand):  ref base T (62/41);  new base G (3/3);  total (65/45)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 110 total aligned reads displayed.

GGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCG  >  NC_000913/774489‑774783
                                                                                                                                                  |                                                                                                                                                    
ggTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCa                                                                                                                                                   >  2:1289181/1‑150 (MQ=255)
   gtgATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGc                                                                                                                                                <  1:1130839/150‑1 (MQ=255)
    tgATTCACGGCGTGGGCTTGCGTCCGCAGAAGGTTCTTTGGCAAGGGATACTGAGCCCATAGCTTGCCGATAGTGTTCTCATTGTCGTGCGGATTCTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCGCCAGGcc                                                                                                                                               <  1:284342/150‑1 (MQ=255)
      aTTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCa                                                                                                                                             >  1:1894836/1‑150 (MQ=255)
       ttCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCAt                                                                                                                                            >  2:657418/1‑150 (MQ=255)
       ttCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCAt                                                                                                                                            >  2:499370/1‑150 (MQ=255)
        tCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATa                                                                                                                                           >  2:212574/1‑150 (MQ=255)
         cACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATaa                                                                                                                                          >  2:1290867/1‑150 (MQ=255)
          aCGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTGTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGACCATaaa                                                                                                                                         >  2:276587/1‑150 (MQ=255)
           cGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAt                                                                                                                                        >  1:868293/1‑150 (MQ=255)
              cGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATGTTTTGCCCAGTGCTTGCCGATAGTGTGCGCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATATAtttt                                                                                                                                     >  2:291735/1‑150 (MQ=255)
                tGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTGCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTa                                                                                                                                   >  1:474744/1‑150 (MQ=255)
                tGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGAAATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTa                                                                                                                                   >  2:1364519/1‑150 (MQ=255)
                      ttCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGTCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTcc                                                                                                                             >  2:1045741/1‑150 (MQ=255)
                      ttCGTCCGCAGAAGGATCTTTGGCAAGGGATTTTTTGCCCATTGCGTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTaaagaga                                                                                                                             <  2:467520/150‑7 (MQ=37)
                        cttCCGCAGAAGGTTCTTTTGCAAGTGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCt                                                                                                                           <  2:1330938/148‑1 (MQ=255)
                         gTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCAtt                                                                                                                          >  2:1419653/1‑150 (MQ=255)
                           ccGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTaa                                                                                                                        >  1:1391246/1‑150 (MQ=255)
                            cGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTaaa                                                                                                                       >  1:1428519/1‑149 (MQ=255)
                                 aaGGTTCTTTGTCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCATTTAAATCAGAATTCTCTTTAAATAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGcc                                                                                                                  <  1:904848/150‑1 (MQ=255)
                                  agcttCTTTGGCAAGCGATTTTATGCCCATTGCTAGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAAACAGAATTCTCTTAAAAAAATTATGTGCCGCTCCCGGCCCATAAATTTTTACGCTCCCGTAACTAGccc                                                                                                                 >  2:220600/4‑150 (MQ=255)
                                  aGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGccc                                                                                                                 >  1:1835137/1‑150 (MQ=255)
                                  aGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGccc                                                                                                                 >  1:1711348/1‑150 (MQ=255)
                                   ggTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCt                                                                                                                >  2:1856069/1‑150 (MQ=255)
                                   ggTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCt                                                                                                                <  1:745900/150‑1 (MQ=255)
                                   ggTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCt                                                                                                                <  1:1461396/150‑1 (MQ=255)
                                     ttCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATAAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCa                                                                                                              >  2:1732496/1‑150 (MQ=255)
                                       cTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCAtt                                                                                                            <  2:551582/150‑1 (MQ=255)
                                         ttGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCATTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCATTAACTTGCCCTCAATcc                                                                                                          >  2:336449/1‑150 (MQ=255)
                                          tGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTccc                                                                                                         >  2:1664236/1‑150 (MQ=255)
                                           ggCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCa                                                                                                        >  2:1644310/1‑150 (MQ=255)
                                             cAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCaac                                                                                                      >  2:421753/1‑149 (MQ=255)
                                             cAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCaaa                                                                                                      >  2:870460/1‑150 (MQ=255)
                                               cgggATGTTTTGCCCATTGCTTGCCGATATTGCTCTCAATGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAAcc                                                                                                    <  2:856203/149‑1 (MQ=255)
                                               aGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAAcc                                                                                                    <  1:463492/150‑1 (MQ=255)
                                                 ggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGGCCTCATTCCCAAACCTc                                                                                                  >  2:1202039/1‑150 (MQ=255)
                                                 ggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCATTAACCTGCCCTCATTCCCAAACCTc                                                                                                  >  2:1196030/1‑150 (MQ=255)
                                                  gATTTCTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCACCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCa                                                                                                 <  2:476473/150‑1 (MQ=255)
                                                   aTTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCaa                                                                                                >  1:1874555/1‑150 (MQ=255)
                                                   aTTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCaa                                                                                                >  1:1867468/1‑150 (MQ=255)
                                                    ttttttGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAAt                                                                                               <  2:971710/150‑1 (MQ=255)
                                                       tttGCCCATTGCTTGCCGATATTGTTCTCATTGTAGGGATGATTTTCTTCTTCTTTTAAATCATAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgc                                                                                            <  2:1050705/150‑1 (MQ=255)
                                                        ttGCCCATTGCTTGCCGAAATTGTTCTCCTCGACGCGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAAATATCGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcg                                                                                           <  2:1035653/150‑1 (MQ=255)
                                                          gCCCATTGCTTGCCGATATTGGTCTCCTTGTCGGGCTGAGTTTCTTCTTCTTTTAACTCAGAATTCTCTTTAAAAACTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTACCTTGCCCTCATTCCCAAACCTCACTcgcgcg                                                                                         <  2:1394988/150‑1 (MQ=255)
                                                             cATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGt                                                                                      <  2:868293/150‑1 (MQ=255)
                                                                ttcttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTATTCTTTTAAATCAGAAATCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCATAACTTGCCCTCATTACCAAACCTCAATCGCGCGCGtata                                                                                   >  2:435822/3‑150 (MQ=255)
                                                                tgcttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGtata                                                                                   <  1:1851665/150‑1 (MQ=255)
                                                                 gcttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtag                                                                                  >  1:552583/1‑150 (MQ=255)
                                                                 gcttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCATTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGtttat                                                                                  >  1:135691/1‑149 (MQ=255)
                                                                  cttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagt                                                                                 <  2:1863112/150‑1 (MQ=255)
                                                                   ttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagta                                                                                <  1:1654425/150‑1 (MQ=255)
                                                                   ttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGACCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagta                                                                                >  1:1269015/1‑150 (MQ=255)
                                                                     gcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTagc                                                                              <  1:1755958/150‑1 (MQ=255)
                                                                      ccGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTagca                                                                             >  2:1366467/1‑150 (MQ=255)
                                                                       cGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTagcag                                                                            >  2:424522/1‑150 (MQ=255)
                                                                             tGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCTCGCGCGTATAGTAGCAGCGTTTa                                                                      >  2:131721/1‑150 (MQ=255)
                                                                              gagctcATAGTCTGGCTGACATTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTAACGGCCGCTCCAGGCCCATAAATTTATACTCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                     <  2:1245738/147‑1 (MQ=255)
                                                                              gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                     >  2:617384/1‑150 (MQ=255)
                                                                                tctcATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaaaa                                                                   >  2:417350/1‑150 (MQ=255)
                                                                                tctcATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaaaa                                                                   >  2:1075967/1‑150 (MQ=255)
                                                                                 ctcATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAg                                                                  >  1:1125139/1‑150 (MQ=255)
                                                                                    aTTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGctc                                                               <  1:1528064/150‑1 (MQ=255)
                                                                                       gTCGGGCTGATTTTCTTCTTCTTTTAAATCAGCATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATACTGTTTTACGCTCCCTTAACTGGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCTGCGTTTAAAAGctctta                                                            >  2:1438094/1‑148 (MQ=255)
                                                                                       gTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCACAACCCTCAATCGCGCGCGTATTGTAGCAGCGTTTTAAAGCTCTaa                                                            >  2:1156696/1‑150 (MQ=255)
                                                                                        tCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAAc                                                           <  2:1405329/150‑1 (MQ=255)
                                                                                            gCTGATTTTCTTCTTCATTTAAATCAGAATTCTCTATAAAAACTACACGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACtttt                                                       <  2:676881/150‑1 (MQ=255)
                                                                                              tGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttgt                                                     <  2:1711348/150‑1 (MQ=255)
                                                                                               gATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttgtt                                                    <  1:1268540/150‑1 (MQ=255)
                                                                                               gATTTTCTACTTCTTTTAAATCAGAATTCTCTTTAAAAAATAATGTGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttgtt                                                    <  1:35848/150‑1 (MQ=255)
                                                                                                   ttcttcttcTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCATTAACTTGCACTCATTCCCAAACATCAATCGAGCGCGTATAGTATAAGCGTTTAAAAGCTCTAACTTTTGTTGCAt                                                >  1:1306757/1‑150 (MQ=255)
                                                                                                      ttcttcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTAc                                             >  2:886685/1‑150 (MQ=255)
                                                                                                       tcttcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTAcc                                            >  1:654620/1‑150 (MQ=255)
                                                                                                          tcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCggg                                         >  2:774569/1‑150 (MQ=255)
                                                                                                          tcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCggg                                         <  1:1856069/150‑1 (MQ=255)
                                                                                                           cttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGa                                        <  1:554343/150‑1 (MQ=255)
                                                                                                           cttTTAAATCAGAATTCTCATTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGa                                        <  2:1161949/150‑1 (MQ=255)
                                                                                                             tttAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATg                                      <  1:1180061/150‑1 (MQ=255)
                                                                                                              ttAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGt                                     <  2:1769945/150‑1 (MQ=255)
                                                                                                                  aTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAg                                 >  2:479694/1‑150 (MQ=255)
                                                                                                                     aGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAACAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGa                              >  2:1438825/1‑150 (MQ=255)
                                                                                                                       aaTTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAAt                            <  2:556496/150‑1 (MQ=255)
                                                                                                                           gagcgtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACaa                        >  1:1480163/6‑150 (MQ=37)
                                                                                                                            tctTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAAc                       <  1:1732496/150‑1 (MQ=255)
                                                                                                                              tttAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCTGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACgg                     >  1:223829/1‑149 (MQ=255)
                                                                                                                              cgtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCTATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCt                    >  2:1722510/3‑150 (MQ=38)
                                                                                                                               ttAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTGATTCCCCACCCTCAGTCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTCACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGcg                    >  1:1887465/1‑149 (MQ=255)
                                                                                                                               ttAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCt                    >  2:1192572/1‑150 (MQ=255)
                                                                                                                               ttAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTAAGCAAACTTAACTAGCCCTCATACCCAAACATCAATAGCGCGCGTAAAGTAGCAGCGTTAAAAATCAATAAATATAGTTGCATTACCGGGATGTAAAGTGATTACACCGCt                    >  2:1516487/1‑150 (MQ=255)
                                                                                                                                  aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGtt                 >  2:30388/1‑150 (MQ=255)
                                                                                                                                   aaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGttt                >  2:1798210/1‑150 (MQ=255)
                                                                                                                                    aaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCGCTTAACTTGCCCTCATTCCCTAACCTCAATCGCGCGCGTATAGAAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCTTTACCGGGATGTAAAGTGAATACAACGCTGTTTc               <  1:336449/150‑1 (MQ=255)
                                                                                                                                    aaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAGCCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTc               >  1:1364937/1‑150 (MQ=255)
                                                                                                                                    aaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTc               <  2:1391246/150‑1 (MQ=255)
                                                                                                                                    aaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTc               <  1:562885/150‑1 (MQ=255)
                                                                                                                                       ttATGGGCAGCTACAGTCCCATAAATTTTTACTCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGATCGTATAGTAGCAGCGTTTAAAAGCTCTAACATTTGTTGCATTACCGGGATGTAAAGCGAATACAACGCTGTTTCGAt            <  2:700492/150‑1 (MQ=255)
                                                                                                                                        tATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATg           <  2:1894836/150‑1 (MQ=255)
                                                                                                                                         aTGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATgg          <  1:1644370/150‑1 (MQ=255)
                                                                                                                                            ggCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCg       >  2:673424/1‑150 (MQ=255)
                                                                                                                                              ccGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGt     <  1:842806/150‑1 (MQ=255)
                                                                                                                                                 cTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTcg  >  1:1118794/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCG  >  NC_000913/774489‑774783

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