Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 759,007 T→G 13.8% F101C (TTC→TGC)  sucA → 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913759,0070TG13.8% 212.3 / 21.8 86F101C (TTC→TGC) sucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (34/40);  new base G (7/5);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAG  >  NC_000913/758859‑759150
                                                                                                                                                    |                                                                                                                                               
ggAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATCTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCtt                                                                                                                                                 <  1:345732/149‑1 (MQ=255)
   gTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCAGTAACGCATACCGCTGCCg                                                                                                                                              >  2:133500/1‑149 (MQ=255)
    tCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACTCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAGGTCCTGCCGCGCATTAACGCATACCGCTTCCGt                                                                                                                                             >  6:258372/1‑149 (MQ=255)
    tCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGt                                                                                                                                             >  1:148826/1‑149 (MQ=255)
    tCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGt                                                                                                                                             >  1:403961/1‑149 (MQ=255)
    tCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGt                                                                                                                                             >  2:359951/1‑149 (MQ=255)
    tCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGt                                                                                                                                             >  7:306175/1‑149 (MQ=255)
    tCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGt                                                                                                                                             >  5:291644/1‑149 (MQ=255)
    tCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTGCCGt                                                                                                                                             <  5:173753/149‑1 (MQ=255)
      aaaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTg                                                                                                                                            >  8:365493/1‑148 (MQ=255)
      aaaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTg                                                                                                                                            <  7:365493/148‑1 (MQ=255)
       aaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGt                                                                                                                                          <  7:340011/149‑1 (MQ=255)
       aaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGt                                                                                                                                          <  3:98789/149‑1 (MQ=255)
        aCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCg                                                                                                                                              <  1:121602/144‑1 (MQ=255)
        aCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCg                                                                                                                                              >  2:121602/1‑144 (MQ=255)
         acGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCGGCAGCTCATTACCGCCTACCGCTTCCGTGGTCa                                                                                                                                        <  4:267657/148‑1 (MQ=255)
          cGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCAc                                                                                                                                       <  8:249866/149‑1 (MQ=255)
            gATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTGCCGTGGTCAc                                                                                                                                       >  2:384289/1‑147 (MQ=255)
                aaTTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCa                                                                                                                                  <  6:68496/148‑1 (MQ=255)
                aaTTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCa                                                                                                                                  <  1:359951/148‑1 (MQ=255)
                aaTTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTGCCGTGGTCACCAGCa                                                                                                                                  <  7:329570/148‑1 (MQ=255)
                     cACTCTCATACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCGGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTGCCGTGGTCACCAGCATGCGaa                                                                                                                            <  6:177576/149‑1 (MQ=255)
                     cACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGaa                                                                                                                            >  8:325351/1‑149 (MQ=255)
                        tctcAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtct                                                                                                                         >  7:145261/1‑149 (MQ=255)
                          tcAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtctc                                                                                                                        >  1:26431/1‑148 (MQ=255)
                          tcAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtct                                                                                                                         <  6:2902/147‑1 (MQ=255)
                            aaaCGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGCCCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGAt                                                                                                                     >  2:212065/1‑149 (MQ=255)
                                gcgTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGc                                                                                                                 <  4:147924/149‑1 (MQ=255)
                                  gTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTg                                                                                                               <  8:25437/149‑1 (MQ=255)
                                  gTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTg                                                                                                               >  2:399301/1‑149 (MQ=255)
                                    gAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTAATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTgg                                                                                                              >  4:46012/1‑148 (MQ=255)
                                     aaTATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGa                                                                                                            <  2:209591/149‑1 (MQ=255)
                                     aaTATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTGCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGa                                                                                                            <  1:384289/149‑1 (MQ=255)
                                         tttCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACtgt                                                                                                        >  8:271625/1‑149 (MQ=255)
                                          ttCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACtgtg                                                                                                       >  8:151618/1‑149 (MQ=255)
                                            ccgccgCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTgg                                                                                                      <  2:196543/148‑1 (MQ=255)
                                              gccgccTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTGCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGgcag                                                                                                   >  1:157460/1‑149 (MQ=255)
                                                      gCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCAGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATAAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGAt                                                                                             <  5:250414/147‑1 (MQ=255)
                                                      gCGAAAGACGCTTCACGTTACTCTTCAACGATCACCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGAt                                                                                             <  4:218677/147‑1 (MQ=255)
                                                        gAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGTAAGATAAAg                                                                                         <  1:352606/149‑1 (MQ=255)
                                                        gAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAg                                                                                         <  1:361321/149‑1 (MQ=255)
                                                        gAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTGCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAg                                                                                         >  6:63688/1‑149 (MQ=255)
                                                        gAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACACATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAg                                                                                         <  4:67446/149‑1 (MQ=255)
                                                                cTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGAt                                                                                 >  4:299575/1‑149 (MQ=255)
                                                                   cACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATc                                                                                <  1:130246/147‑1 (MQ=255)
                                                                         aCTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGTAAGATAAAGTGGCCGATCTGGATCCg                                                                        >  5:129851/1‑149 (MQ=255)
                                                                           tctTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTc                                                                      >  3:163679/1‑149 (MQ=255)
                                                                             ttCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCCGCTCATTATCGCATACCGCTTCCGTGGTCCCCAGCATGCGAATCTCGATGCGCTGGGACTGTGGCAGCAAGCTAAAGTGGCCGATCTGGAGACGTCtt                                                                    >  6:386111/1‑149 (MQ=255)
                                                                                aaCGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGAt                                                                                             >  1:115541/1‑121 (MQ=255)
                                                                                aaCGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGAt                                                                                             <  2:115541/121‑1 (MQ=255)
                                                                                        ccGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGCACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCt                                                          >  7:189506/1‑148 (MQ=255)
                                                                                         cgccccTGACACCACTGTGCAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTCGATCCGTCTTTCCACGATc                                                           <  6:291644/143‑1 (MQ=255)
                                                                                         cGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCTTGCGAATCTCGATCCGCTGTGACTGTTGCATCAAGTTAAAGTGGCCGATCTGGATCCGTCTTTCCACGATc                                                           <  2:148826/146‑1 (MQ=255)
                                                                                         cGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGa                                                             <  1:399301/144‑1 (MQ=255)
                                                                                         cGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGa                                                        <  4:163679/149‑1 (MQ=255)
                                                                                         cGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGa                                                        <  7:342914/149‑1 (MQ=255)
                                                                                             ccTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACCCATACCGCTGCCGTGCTCACCAGCATGCGAAACTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGAACTGGATCCGTCTTTCCACGATCTCACCg                                                     >  1:159428/1‑148 (MQ=255)
                                                                                              cTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGgc                                                                                                     <  8:343716/99‑1 (MQ=255)
                                                                                              cTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGgc                                                                                                     >  7:343716/1‑99 (MQ=255)
                                                                                               tGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAg                                                  >  6:375872/1‑149 (MQ=255)
                                                                                               tGACACCAATGTGAAGCAGGTTAAACTCCTGCAGCTCATAAACGCATACCGCTTCCGTGGTCACCAGCAAGCGAATCTCGATCCGCTGGCACTGTGGCAGCAAGATAAACTGGCCCATCTCGATCCGTCTATCCACTATCTGACCGaac                                                  >  5:256057/1‑148 (MQ=255)
                                                                                                gACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTGCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGActgccac                                                 >  7:345505/1‑143 (MQ=255)
                                                                                                 acacCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCa                                                <  7:325351/149‑1 (MQ=255)
                                                                                                 acacCAATGTGAAGCAGGTTAAAGGCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGCGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCa                                                <  6:94769/149‑1 (MQ=255)
                                                                                                   acCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGAt                                                                           >  2:64375/1‑120 (MQ=255)
                                                                                                   acCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGAt                                                                           <  1:64375/120‑1 (MQ=255)
                                                                                                   acCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGa                                              <  3:315404/149‑1 (MQ=255)
                                                                                                    ccAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGAc                                             <  7:271625/149‑1 (MQ=255)
                                                                                                        tgtgAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGCTCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTcc                                         <  4:237961/149‑1 (MQ=255)
                                                                                                          tgAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCa                                        <  4:295095/148‑1 (MQ=255)
                                                                                                          tgAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAg                                       <  3:38007/149‑1 (MQ=255)
                                                                                                          ggacgCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGGGGTCACCCGCATGCGAAGCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCa                                        <  5:86844/144‑1 (MQ=255)
                                                                                                           gAAGCAGGTTTAAGTCCTGCAGCTCATTAACGCGTACCGCTTCCGTGGTCACCGGCATGCGTATGTCGATCCGCGGGGCGTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGGCTTTCCACGATCTGACCGAAGCAGACTTCCAgg                                      >  8:171791/1‑149 (MQ=255)
                                                                                                                 ggTTAAAGTCCTGCCGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTCGCCGATCTGGATACGTCTTTCCACGATATGAACGAAGCAGACTTCCAGgagcc                                  >  5:324295/1‑145 (MQ=255)
                                                                                                                    tAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTc                              >  1:16741/1‑148 (MQ=255)
                                                                                                                     aaaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCa                             <  1:253977/148‑1 (MQ=255)
                                                                                                                     aaaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCAAGGAGACCTTCaa                            >  4:114002/1‑149 (MQ=255)
                                                                                                                      aaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAAc                           <  2:26431/149‑1 (MQ=255)
                                                                                                                             gCAGCTCATTAACGCATACCGCTGCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGtt                    >  6:389674/1‑149 (MQ=255)
                                                                                                                                 cTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCAtt                <  2:403961/149‑1 (MQ=255)
                                                                                                                                  tCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTGTTTCCCCGATCTGACCGAAGCAGACTTCCAGGAGCCCTTCAACGTCGGTTCAttt               >  8:90923/1‑149 (MQ=255)
                                                                                                                                     ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGCTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGc             >  1:265828/1‑148 (MQ=255)
                                                                                                                                     ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGc             >  3:324120/1‑148 (MQ=255)
                                                                                                                                     ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGc             >  8:70861/1‑148 (MQ=255)
                                                                                                                                     ttAACGCATACCGCTGCCTTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGcc            <  7:321119/149‑1 (MQ=255)
                                                                                                                                               ccGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCaaa   <  8:166207/148‑1 (MQ=255)
                                                                                                                                               ccGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGGCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAg  <  7:151618/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAG  >  NC_000913/758859‑759150

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