Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 336,866 G→C 7.6% G30A (GGC→GCC)  yahF → putative NAD(P)‑binding succinyl‑CoA synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913336,8660GC 7.6% 268.6 / 12.3 79G30A (GGC→GCC) yahFputative NAD(P)‑binding succinyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base G (36/37);  new base C (3/3);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

TGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGG  >  NC_000913/336724‑336998
                                                                                                                                              |                                                                                                                                    
tGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCg                                                                                                                                    >  1:54828/1‑145 (MQ=255)
tGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCg                                                                                                                                    <  2:54828/145‑1 (MQ=255)
    gTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCa                                                                                                                             <  3:949937/148‑1 (MQ=255)
          catggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAtt                                                                                                                       >  1:130600/1‑148 (MQ=255)
          catggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAtt                                                                                                                       >  4:1100058/1‑148 (MQ=255)
                        gtgCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatggcga                                                                                                         >  2:1091730/1‑148 (MQ=255)
                                 gCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGa                                                                                                <  3:1100058/148‑1 (MQ=255)
                                     aaCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATaa                                                                                            <  1:897268/148‑1 (MQ=255)
                                           cGGAGGGAATTATGTCAGTTAAAATAGTCATTAAAACGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTc                                                                                                                                  >  4:963596/1‑104 (MQ=255)
                                           cGGAGCGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTc                                                                                                                                  <  3:963596/104‑1 (MQ=255)
                                             gAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAg                                                                                                                               >  4:911366/1‑105 (MQ=255)
                                             gAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAg                                                                                                                               <  3:911366/105‑1 (MQ=255)
                                                ggAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTga                                                                                 <  1:1091730/148‑1 (MQ=255)
                                                 gAATTATGTCAGTTAAAATAGTCATTAAACCGGATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                            <  3:539102/104‑1 (MQ=255)
                                                 gAATTATGTCAGTTAAAATAGTCATTAAACCGGATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                            >  4:539102/1‑104 (MQ=255)
                                                 gAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAtt                                                                                                                       <  2:1001297/109‑1 (MQ=255)
                                                 gAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAtt                                                                                                                       >  1:1001297/1‑109 (MQ=255)
                                                        gTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                         >  2:717075/1‑148 (MQ=255)
                                                          cAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATaaa                                                                                           <  3:863905/128‑1 (MQ=255)
                                                          cAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATaaa                                                                                           >  4:863905/1‑128 (MQ=255)
                                                          cAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                       <  2:130600/148‑1 (MQ=255)
                                                                    aGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCa                                                                                                                         >  3:918919/1‑88 (MQ=255)
                                                                    aGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCa                                                                                                                         <  4:918919/88‑1 (MQ=255)
                                                                      tCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGcc                                                           <  1:717075/148‑1 (MQ=255)
                                                                         ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaa                                                                                                    >  1:793700/1‑104 (MQ=255)
                                                                         ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaa                                                                                                    <  2:793700/104‑1 (MQ=255)
                                                                         ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaaga                                                                              <  4:362449/126‑1 (MQ=255)
                                                                         ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaaga                                                                              >  3:362449/1‑126 (MQ=255)
                                                                         ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaaga                                                                              >  1:19676/1‑126 (MQ=255)
                                                                         ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaaga                                                                              <  2:19676/126‑1 (MQ=255)
                                                                          tAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTACCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATg                                                                                                 <  3:23746/106‑1 (MQ=255)
                                                                          tAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTACCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATg                                                                                                 >  4:23746/1‑106 (MQ=255)
                                                                                 aaTACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatg                                                                                                             >  1:255156/1‑87 (MQ=255)
                                                                                 aaTACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatg                                                                                                             <  2:255156/87‑1 (MQ=255)
                                                                                  aTACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaaa                                                                                                   >  1:626292/1‑96 (MQ=255)
                                                                                  aTACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaaa                                                                                                   <  2:626292/96‑1 (MQ=255)
                                                                                       tATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGa                                          >  1:870163/1‑148 (MQ=255)
                                                                                          tttGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGa                                                                                                >  3:1094510/1‑91 (MQ=255)
                                                                                          tttGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGa                                                                                                <  4:1094510/91‑1 (MQ=255)
                                                                                               ttCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatggcg                                                                                                          <  4:80639/76‑1 (MQ=255)
                                                                                               ttCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatggcg                                                                                                          >  3:80639/1‑76 (MQ=255)
                                                                                                tCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGt                                                                                                                                   <  2:857792/50‑1 (MQ=255)
                                                                                                tCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGt                                                                                                                                   >  1:857792/1‑50 (MQ=255)
                                                                                                tCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                         <  4:198014/108‑1 (MQ=255)
                                                                                                tCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                         >  3:198014/1‑108 (MQ=255)
                                                                                                     ctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                         >  3:753972/1‑87 (MQ=255)
                                                                                                     ctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAttt                                                                          >  4:943230/1‑102 (MQ=255)
                                                                                                     ctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAttt                                                                          <  3:943230/102‑1 (MQ=255)
                                                                                                     ctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCACTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                         <  4:753972/87‑1 (MQ=255)
                                                                                                           gATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                       <  1:235407/99‑1 (MQ=255)
                                                                                                           gATGTCTATCTCCACGCGGGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                       >  2:235407/1‑99 (MQ=255)
                                                                                                             tGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGt                    >  4:989630/1‑148 (MQ=255)
                                                                                                                cTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGAcc                              >  3:559878/1‑135 (MQ=255)
                                                                                                                cTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGAcc                              <  4:559878/135‑1 (MQ=255)
                                                                                                                cTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGtcat                 <  3:989630/148‑1 (MQ=255)
                                                                                                                 tATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCa                                               >  4:447079/1‑117 (MQ=255)
                                                                                                                 tATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCa                                               <  3:447079/117‑1 (MQ=255)
                                                                                                                  aTCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaaa                                                                                                   <  1:680150/64‑1 (MQ=255)
                                                                                                                  aTCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaaa                                                                                                   >  2:680150/1‑64 (MQ=255)
                                                                                                                      ccACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaa                                                                                                    >  3:589936/1‑59 (MQ=255)
                                                                                                                      ccACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaa                                                                                                    <  4:589936/59‑1 (MQ=255)
                                                                                                                       cACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAGGGCg                                                                                       >  2:493181/1‑71 (MQ=255)
                                                                                                                       cACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCg                                                                                       <  1:493181/71‑1 (MQ=255)
                                                                                                                        aCGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTg                                                                                  >  1:77311/1‑75 (MQ=255)
                                                                                                                        aCGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTg                                                                                  <  2:77311/75‑1 (MQ=255)
                                                                                                                         cgcgTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                         <  1:565666/67‑1 (MQ=255)
                                                                                                                         cgcgTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                         >  2:565666/1‑67 (MQ=255)
                                                                                                                                 aaTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                       <  2:884625/77‑1 (MQ=255)
                                                                                                                                 aaTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAAGAAAGGCGTGCAGAAGAATTTAgg                                                                       >  1:884625/1‑77 (MQ=255)
                                                                                                                                  aTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTgatgat                       >  4:524963/1‑124 (MQ=255)
                                                                                                                                  aTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTgatgat                       <  3:524963/124‑1 (MQ=255)
                                                                                                                                      gcTCGACGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                       <  4:311699/71‑1 (MQ=255)
                                                                                                                                      gcTCGACGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                       >  3:311699/2‑72 (MQ=255)
                                                                                                                                       cTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCa                                               >  4:428470/1‑95 (MQ=255)
                                                                                                                                       cTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCa                                               <  3:428470/95‑1 (MQ=255)
                                                                                                                                        tCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAt                                                                            >  1:583611/1‑65 (MQ=255)
                                                                                                                                        tCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAt                                                                            <  2:583611/65‑1 (MQ=255)
                                                                                                                                           aCGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCggg  <  2:633108/136‑1 (MQ=255)
                                                                                                                                           aCGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCggg  >  1:633108/1‑136 (MQ=255)
                                                                                                                                              |                                                                                                                                    
TGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGG  >  NC_000913/336724‑336998

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