Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,757 T→C 7.1% intergenic (‑140/+24) pdeB ← / ← ylaC c‑di‑GMP phosphodiesterase PdeB/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,7570TC 7.1% 549.9 / 22.9 183intergenic (‑140/+24)pdeB/ylaCc‑di‑GMP phosphodiesterase PdeB/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (98/72);  new base C (11/2);  total (109/74)
Fisher's exact test for biased strand distribution p-value = 7.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 185 total aligned reads displayed.

GTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAGATGG  >  NC_000913/478610‑478905
                                                                                                                                                   |                                                                                                                                                    
   ctcACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGc                                                                                                                                                 >  1:1003122/1‑150 (MQ=255)
    tcACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                                                                <  1:1518164/150‑1 (MQ=255)
      aCCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                                                                              >  1:1218435/1‑150 (MQ=255)
        cAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                                                                            <  2:1887054/150‑1 (MQ=255)
         aGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                                                                           <  1:265654/150‑1 (MQ=255)
          gCATATCCTTATTAACACAGTGTAGTGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTAt                                                                                                                                          >  2:2461877/1‑150 (MQ=255)
          ccATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                                                                          <  1:1444706/149‑3 (MQ=255)
             tatCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                                                                       >  1:480431/1‑150 (MQ=255)
             tatCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                                                                       >  1:253787/1‑150 (MQ=255)
             tatCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                                                                       >  1:1831410/1‑150 (MQ=255)
               tCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                                                                     >  1:1808155/1‑150 (MQ=255)
                      taACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACt                                                                                                                              >  2:108734/1‑150 (MQ=255)
                         cacaGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                                                                           <  2:538329/150‑1 (MQ=255)
                           caGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                                                                         <  1:1317116/150‑1 (MQ=255)
                           caGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                                                                         >  1:1979805/1‑150 (MQ=255)
                            aGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGAt                                                                                                                        <  1:1893570/150‑1 (MQ=255)
                             gtgtAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATg                                                                                                                       >  1:370605/1‑150 (MQ=255)
                               gtAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                                                                     >  1:1946837/1‑150 (MQ=255)
                                tAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                                                                    >  1:1291704/1‑150 (MQ=255)
                                tAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                                                                    <  1:967286/150‑1 (MQ=255)
                                 aGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                                                                                   >  1:277630/1‑150 (MQ=255)
                                 aGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                                                                                   >  2:1527612/1‑150 (MQ=255)
                                 aGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                                                                                   >  1:1476816/1‑150 (MQ=255)
                                  gggAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                                                                  >  2:1720985/1‑150 (MQ=255)
                                    gAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                                                                                >  2:1669600/1‑150 (MQ=255)
                                            atCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                                                                        >  2:365552/1‑150 (MQ=255)
                                            atCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                                                                        <  1:1157486/150‑1 (MQ=255)
                                               tGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                                                                     >  2:997214/1‑150 (MQ=255)
                                               tGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                                                                     >  1:703114/1‑150 (MQ=255)
                                                gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                                                                    <  2:19917/150‑1 (MQ=255)
                                                  ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                                                                  <  2:119563/150‑1 (MQ=255)
                                                  ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                                                                  >  2:2102511/1‑150 (MQ=255)
                                                   tgtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                                                                 >  2:804745/1‑150 (MQ=255)
                                                   tgtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                                                                 >  2:2426782/1‑150 (MQ=255)
                                                      ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                                                                              >  2:1363815/1‑150 (MQ=255)
                                                        aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaaa                                                                                            <  2:1808155/150‑1 (MQ=255)
                                                        aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaaa                                                                                            <  2:1476816/150‑1 (MQ=255)
                                                             ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                                                                       <  1:1267708/150‑1 (MQ=255)
                                                              gctgggtgGAGATAAGTAGGTAGAGCTTGTATTTCATGAAGTGATGGGGGTAAACGATCATTGTCATTTGACCATGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGTTAAAAATATCGt                                                                                      <  2:1575323/150‑1 (MQ=255)
                                                              gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                                                                                      <  1:1882707/150‑1 (MQ=255)
                                                              gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                                                                                      <  2:277630/150‑1 (MQ=255)
                                                                tggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaa                                                                                    >  1:579792/1‑150 (MQ=255)
                                                                 ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaa                                                                                   >  1:1571078/1‑150 (MQ=255)
                                                                      gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                                                                              >  2:683979/1‑150 (MQ=255)
                                                                        gaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                                                                            >  1:997557/1‑150 (MQ=255)
                                                                          tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                                                                          >  2:11248/1‑150 (MQ=255)
                                                                           aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                                                                         <  1:357185/150‑1 (MQ=255)
                                                                            aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                                                                        >  1:1702050/1‑150 (MQ=255)
                                                                            aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                                                                        >  1:2339931/1‑150 (MQ=255)
                                                                             ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                                                                       <  2:1003122/150‑1 (MQ=255)
                                                                                   agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                                                                 <  2:1226382/150‑1 (MQ=255)
                                                                                    gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                                                                <  2:1571078/150‑1 (MQ=255)
                                                                                      gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa                                                              >  2:817052/1‑150 (MQ=255)
                                                                                      gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa                                                              >  1:1964566/1‑150 (MQ=255)
                                                                                       cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc                                                             <  2:2339931/150‑1 (MQ=255)
                                                                                        tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg                                                            >  1:1379542/1‑150 (MQ=255)
                                                                                        tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg                                                            <  1:2455998/150‑1 (MQ=255)
                                                                                          gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                                                                     >  2:2263355/1‑139 (MQ=255)
                                                                                          gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                                                                     <  1:2263355/139‑1 (MQ=255)
                                                                                            aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa                                                        <  1:108734/150‑1 (MQ=255)
                                                                                            aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa                                                        >  2:1318994/1‑150 (MQ=255)
                                                                                            aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa                                                        <  2:65519/150‑1 (MQ=255)
                                                                                             tttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAt                                                       <  1:1482998/150‑1 (MQ=255)
                                                                                              ttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt                                                      <  2:1796217/150‑1 (MQ=255)
                                                                                                cATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt                                                    <  1:1005715/150‑1 (MQ=255)
                                                                                                cATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt                                                    <  1:1014050/150‑1 (MQ=255)
                                                                                                    tattGATGGGTGGAAACGATCATTGTCATTTTACCAGGCGCAGAGCTTTACGCCTGTTATGTAACTCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGc                                                <  1:2133496/147‑1 (MQ=255)
                                                                                                    aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGc                                                <  1:828585/150‑1 (MQ=255)
                                                                                                     aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCa                                               >  2:1545427/1‑150 (MQ=255)
                                                                                                     aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCa                                               >  1:862795/1‑150 (MQ=255)
                                                                                                       tGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAgt                                             <  2:370605/150‑1 (MQ=255)
                                                                                                        gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAgtt                                            >  2:883859/1‑150 (MQ=255)
                                                                                                        gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAATTTTTGCAgtt                                            >  2:2316027/1‑150 (MQ=255)
                                                                                                           ttctGGAAACTATCATTGTCATTTGACCAGGCTCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGCTTTTGATTCGGCGCGAGCGAGTGTAATAATATCGAAAAAAGAAAGTTCACCTTTACGGACTATCATTTATTGCAgttggt                                         <  2:1257955/147‑3 (MQ=255)
                                                                                                             gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                                       >  1:589518/1‑150 (MQ=255)
                                                                                                                   aCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGa                                 <  1:1658476/150‑1 (MQ=255)
                                                                                                                    cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGaa                                >  2:1843363/1‑150 (MQ=255)
                                                                                                                    cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGaa                                <  1:1843363/150‑1 (MQ=255)
                                                                                                                     gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAAc                               >  1:669259/1‑150 (MQ=255)
                                                                                                                        cATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATc                            <  1:683979/150‑1 (MQ=255)
                                                                                                                         attttaaTTTTACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGAGAGGCTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATcc                           >  2:768770/7‑150 (MQ=255)
                                                                                                                           tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCcgc                         <  2:1831410/150‑1 (MQ=255)
                                                                                                                            gTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCcgcg                        <  1:479687/150‑1 (MQ=255)
                                                                                                                                 cataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcg                   >  2:160281/4‑150 (MQ=255)
                                                                                                                                   taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcgcg                 >  1:1484575/2‑150 (MQ=255)
                                                                                                                                   taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcgcg                 >  2:826927/2‑150 (MQ=255)
                                                                                                                                    aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCgg                >  1:2303583/1‑150 (MQ=255)
                                                                                                                                    aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCgg                <  2:1318205/150‑1 (MQ=255)
                                                                                                                                    aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCgg                <  2:1862874/150‑1 (MQ=255)
                                                                                                                                     ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGga               >  1:979584/1‑150 (MQ=255)
                                                                                                                                     ccACGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTTAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTATACGGACGATCATTTTTTGCAGAAGTTCCTTATCAACATCCGCGACACGCGga               >  2:2423765/1‑150 (MQ=255)
                                                                                                                                      cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGgag              >  1:94269/1‑150 (MQ=255)
                                                                                                                                           gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTc         >  1:2173976/1‑150 (MQ=255)
                                                                                                                                           gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTc         >  2:841355/1‑150 (MQ=255)
                                                                                                                                             agagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCaa       <  1:1545427/150‑1 (MQ=255)
                                                                                                                                              gagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAg      >  1:2015546/1‑150 (MQ=255)
                                                                                                                                              gagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAg      >  1:64066/1‑150 (MQ=255)
                                                                                                                                               agCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAGa     <  2:464433/150‑1 (MQ=255)
                                                                                                                                                gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAGAt    <  2:613128/150‑1 (MQ=255)
                                                                                                                                                gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAGAt    >  1:1437844/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAGATGG  >  NC_000913/478610‑478905

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