Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 736,382 T→C intergenic (+163/‑63) ybfC → / → ybfQ putative secreted protein/pseudogene, H repeat‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913736,3820TC100.0% 166.4 / NA 44intergenic (+163/‑63)ybfC/ybfQputative secreted protein/pseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/23);  total (21/23)

CAAATTAATCATTGTGAAAAATTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGTCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGAT  >  NC_000913/736238‑736528
                                                                                                                                                |                                                                                                                                                  
cAAATTAATCATTGTGAAAAATTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTaa                                                                                                                                                >  1:72962/1‑149 (MQ=255)
    ttAATCATTGTGAAAAATTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGt                                                                                                                                          >  1:472111/1‑151 (MQ=255)
      aaTCATTGTGAAAAATTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGttt                                                                                                                                        <  1:618016/151‑1 (MQ=255)
        tCATTGTGAAAAATTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATATTACTTATTGCCTAAAGCTGttttt                                                                                                                                      >  2:407133/1‑151 (MQ=255)
          aTTGTGAAAAATTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTc                                                                                                                                     >  2:440849/1‑150 (MQ=255)
              tgAAAAATTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTACCTAAAGCTGttttt                                                                                                                                      >  2:635878/1‑145 (MQ=255)
                 aaaaTTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTta                                                                                                                             <  1:186192/151‑1 (MQ=255)
                   aaTTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTtata                                                                                                                           <  2:544857/151‑1 (MQ=255)
                       atataACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATg                                                                                                                       >  1:292629/1‑151 (MQ=255)
                            aCTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATc                                                                                                                    >  2:256390/1‑149 (MQ=255)
                               cATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAAt                                                                                                                <  1:36555/150‑1 (MQ=255)
                               cATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATa                                                                                                               <  1:531011/151‑1 (MQ=255)
                                     aTAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCtt                                                                                                         <  2:72959/151‑1 (MQ=255)
                                       aGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTAc                                                                                                       >  1:416743/1‑151 (MQ=255)
                                        ggTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACa                                                                                                      <  2:356304/151‑1 (MQ=255)
                                        ggTGTTTATTTTACAGTCCAATTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCt                                                                                                          <  2:472099/147‑1 (MQ=14)
                                           gtTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACata                                                                                                    <  1:260684/150‑1 (MQ=255)
                                                 tttGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAgg                                                                                             >  1:234102/1‑151 (MQ=255)
                                                    gCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTg                                                                                           <  2:266447/150‑1 (MQ=255)
                                                           aTTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATAATCAAATACTCCTTACATAAATAAGGGGAACAAAt                                                                                    >  2:245836/1‑150 (MQ=255)
                                                             ttGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATgg                                                                                  <  1:111665/150‑1 (MQ=255)
                                                               gCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGaa                                                                                <  2:381375/150‑1 (MQ=255)
                                                                cAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAgagatcggaatag                                                                                  >  2:196610/1‑135 (MQ=255)
                                                                       aGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTaaaaaa                                                                       <  2:240428/151‑1 (MQ=255)
                                                                           cACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTg                                                                    <  1:308973/150‑1 (MQ=255)
                                                                               tAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATgg                                                                >  1:95092/1‑150 (MQ=255)
                                                                                aaCGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGa                                                               >  2:558888/1‑150 (MQ=255)
                                                                                 aCGCTACTTTTTTGTACTAAATAATTCGAATCTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACGCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGaa                                                              <  2:84517/150‑1 (MQ=255)
                                                                                     tACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACatat                                                         >  1:126139/1‑151 (MQ=255)
                                                                                             tGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTAtt                                                  <  1:488605/150‑1 (MQ=255)
                                                                                                   aaaTAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCg                                           >  2:310107/1‑151 (MQ=255)
                                                                                                       aatTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTa                                       <  2:555854/151‑1 (MQ=255)
                                                                                                         ttCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTAc                                      >  1:186662/1‑150 (MQ=255)
                                                                                                                tttAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGc                              >  2:366636/1‑151 (MQ=255)
                                                                                                                tttAAGTTCAACAATTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTAGTACACCCAATgacagaca                                 >  2:386337/1‑148 (MQ=255)
                                                                                                                           caaTTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTAGACCCGATTACAGACAAGCCTGGAAAgt                     >  2:48847/1‑149 (MQ=255)
                                                                                                                           caaTTGAGATACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAgtgg                   >  2:451266/1‑151 (MQ=255)
                                                                                                                                gagaTACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCa               <  1:584566/150‑1 (MQ=255)
                                                                                                                                gagaTACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCAt              <  1:366646/151‑1 (MQ=255)
                                                                                                                                  gaTACTACTTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATaa            <  1:229754/151‑1 (MQ=255)
                                                                                                                                        actTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTg        <  1:30091/149‑1 (MQ=255)
                                                                                                                                        actTATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGAGCATAAATTGt       >  1:543657/1‑150 (MQ=255)
                                                                                                                                           tATTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAg    <  2:242931/150‑1 (MQ=255)
                                                                                                                                            aTTGCCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGat  <  1:607757/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CAAATTAATCATTGTGAAAAATTATATAACTCATGGAATAGGTGTTTATTTTGCGTTCCATTTGCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGAATTTTAAGTTCAACAATTGAGATACTACTTATTGTCTAAAGCTGTTTTTCATAGCTTATACATGATCAAATACTCCTTACATAAATAAGGTGAACAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGAT  >  NC_000913/736238‑736528

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