Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 759,716 T→C R337R (CGT→CGC 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,7160TC98.9% 349.0 / ‑4.3 93R337R (CGT→CGCsucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (50/42);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

AACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGTCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTAT  >  NC_000913/759568‑759862
                                                                                                                                                    |                                                                                                                                                  
aaCACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCt                                                                                                                                                  <  2:95225/151‑1 (MQ=255)
aaCACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCt                                                                                                                                                  <  4:110287/151‑1 (MQ=255)
  cacCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTgg                                                                                                                                                >  2:75403/1‑151 (MQ=255)
    ccTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGAc                                                                                                                                              >  3:67883/1‑151 (MQ=255)
     cTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGAc                                                                                                                                              >  4:1386593/1‑150 (MQ=255)
         gCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACt                                                                                                                                         >  4:501464/1‑151 (MQ=255)
           aCGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTtg                                                                                                                                       >  2:766596/1‑151 (MQ=255)
             gggTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTtgat                                                                                                                                     >  3:349999/1‑151 (MQ=255)
                tGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGCGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTtgat                                                                                                                                     >  3:176171/1‑148 (MQ=255)
                  aCGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGCTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGcc                                                                                                                                <  1:332580/151‑1 (MQ=255)
                  aCGTGAAATACCACATGGGCTTCTCGTCTGACCTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGcc                                                                                                                                >  3:107651/1‑151 (MQ=255)
                       aaaTACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGagc                                                                                                                            >  4:1620747/1‑150 (MQ=255)
                        aaTACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGagc                                                                                                                            >  2:980283/1‑149 (MQ=255)
                          tACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGagcagc                                                                                                                         >  2:1607096/1‑150 (MQ=255)
                           aCCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAgcagca                                                                                                                        >  2:1277292/1‑150 (MQ=255)
                            ccACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGcaac                                                                                                                      <  2:1129450/151‑1 (MQ=255)
                            ccACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGcaac                                                                                                                      >  3:304000/1‑151 (MQ=255)
                             cacaTACGCTTCTCGTCTGACTTCCATACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGGTTGTAAGCCCGGGCGTTATCAGTGCTGTTCGTGCCCGCATGGACAGTCTTGATGCGCCGAacagca                                                                                                                        >  4:587727/1‑148 (MQ=255)
                               caTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGTCTGGACAGACTTGATGAGCCGAGCAGCAACaaa                                                                                                                   <  2:428414/151‑1 (MQ=255)
                               caTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGcaacaa                                                                                                                    >  2:1663661/1‑150 (MQ=255)
                               caTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACaaa                                                                                                                   <  4:894527/151‑1 (MQ=255)
                                   ggCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGtgc                                                                                                               >  3:257877/1‑151 (MQ=255)
                                   ggCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGtgc                                                                                                               >  4:1652284/1‑151 (MQ=255)
                                       tctcGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGtgctgc                                                                                                            <  2:229083/150‑1 (MQ=255)
                                       tctcGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGtgctgc                                                                                                            <  4:326290/150‑1 (MQ=255)
                                        ctcGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGcc                                                                                                           >  2:507758/1‑150 (MQ=255)
                                          cGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCa                                                                                                          <  4:639993/149‑1 (MQ=255)
                                           gTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATc                                                                                                       <  2:319069/151‑1 (MQ=255)
                                               gACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCAcc                                                                                                    >  2:81589/1‑150 (MQ=255)
                                               gACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCAcc                                                                                                    >  4:545732/1‑150 (MQ=255)
                                                 cTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATc                                                                                                 >  2:995669/1‑151 (MQ=255)
                                                   tCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATcc                                                                                                <  3:305340/150‑1 (MQ=255)
                                                      aGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACgg                                                                                            >  2:892671/1‑151 (MQ=255)
                                                           gATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGAcg                                                                                       >  4:612595/1‑151 (MQ=255)
                                                                cggcCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCa                                                                                  <  3:1284/151‑1 (MQ=255)
                                                                       gTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCg                                                                           <  1:67315/151‑1 (MQ=255)
                                                                         gCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCggg                                                                         >  1:35180/1‑151 (MQ=255)
                                                                         gCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCggg                                                                         <  4:680701/151‑1 (MQ=255)
                                                                            ccTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCTCGGTGACGCCGCAGTGACCGGGCAg                                                                      <  2:475185/151‑1 (MQ=255)
                                                                            ccTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCa                                                                       <  1:349982/150‑1 (MQ=255)
                                                                            ccTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGCAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGc                                                                        <  2:171938/149‑1 (MQ=255)
                                                                             ctggcgctggcgTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGCGACGCCGCAGTGACCGGGCAgg                                                                     >  4:96283/1‑151 (MQ=255)
                                                                              tggcgctggcgTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAggg                                                                    >  1:225178/1‑151 (MQ=255)
                                                                                gcgctggcgTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCg                                                                  >  4:12537/1‑151 (MQ=255)
                                                                                 cgctggcgTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGt                                                                 >  4:994814/1‑151 (MQ=255)
                                                                                 cgctggcgTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGt                                                                 >  4:1267921/1‑151 (MQ=255)
                                                                                     ggcgTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTgg                                                               <  4:299158/149‑1 (MQ=255)
                                                                                     ggcgTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTgg                                                               <  4:764874/149‑1 (MQ=255)
                                                                                      gcgTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTc                                                            <  2:626842/151‑1 (MQ=255)
                                                                                      gcgTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTc                                                            <  4:696412/151‑1 (MQ=255)
                                                                                         tttAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAgg                                                         >  1:375002/1‑151 (MQ=255)
                                                                                          ttAACCCGTCTCTCCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAgg                                                         <  4:10277/150‑1 (MQ=255)
                                                                                          ttAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAgg                                                         <  1:271017/150‑1 (MQ=255)
                                                                                             aCCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAAc                                                     >  2:1575877/1‑151 (MQ=255)
                                                                                               ccGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAAccc                                                   >  4:930663/1‑151 (MQ=255)
                                                                                                 gTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTg                                                 <  1:6954/151‑1 (MQ=255)
                                                                                                 gTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTg                                                 >  4:579884/1‑151 (MQ=255)
                                                                                                  tctcACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTg                                                 >  2:241449/1‑150 (MQ=255)
                                                                                                   ctcACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGa                                                >  2:250888/1‑150 (MQ=255)
                                                                                                     cACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAAc                                              >  3:106713/1‑150 (MQ=255)
                                                                                                       ccTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACAt                                            <  4:1397738/150‑1 (MQ=255)
                                                                                                          tGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTc                                         <  4:1645146/150‑1 (MQ=255)
                                                                                                          tGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTc                                         <  2:1376069/150‑1 (MQ=255)
                                                                                                           gagaTTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGt                                          <  4:936785/148‑1 (MQ=255)
                                                                                                           gagaTTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGa                                       <  4:767616/151‑1 (MQ=255)
                                                                                                           gagaTTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGa                                       <  1:399958/151‑1 (MQ=255)
                                                                                                               ttGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAgc                                   >  4:873172/1‑151 (MQ=255)
                                                                                                               ttGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAgc                                   >  4:305406/1‑151 (MQ=255)
                                                                                                               ttGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAg                                    <  2:611760/150‑1 (MQ=255)
                                                                                                                tGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAgc                                   <  2:1029198/150‑1 (MQ=255)
                                                                                                                  tAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGcgcg                                >  1:137032/1‑151 (MQ=255)
                                                                                                                    aGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTg                              >  4:1382464/1‑151 (MQ=255)
                                                                                                                    aGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTg                              >  2:96805/1‑151 (MQ=255)
                                                                                                                        cGGTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTa                          <  1:420021/151‑1 (MQ=255)
                                                                                                                         ggTAGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTAt                         >  3:428401/1‑151 (MQ=255)
                                                                                                                            aGTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGaa                      <  2:1618774/151‑1 (MQ=255)
                                                                                                                             gTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGaa                      <  4:1124675/150‑1 (MQ=255)
                                                                                                                             gTTATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGa                       >  1:292524/1‑149 (MQ=255)
                                                                                                                               tATCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGt                    <  1:106699/150‑1 (MQ=255)
                                                                                                                                 tCGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTgg                 >  2:1589565/1‑151 (MQ=255)
                                                                                                                                  cGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGc                <  4:450456/151‑1 (MQ=255)
                                                                                                                                  cGGTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGc                <  4:1278903/151‑1 (MQ=255)
                                                                                                                                   ggTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGcg               >  2:698881/1‑151 (MQ=255)
                                                                                                                                   ggTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGcg               >  2:719472/1‑151 (MQ=255)
                                                                                                                                    gTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGcgg              <  3:419559/151‑1 (MQ=255)
                                                                                                                                    gTTCTGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGcgg              >  3:252233/1‑151 (MQ=255)
                                                                                                                                        tGTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGcggtacg          >  4:169985/1‑151 (MQ=255)
                                                                                                                                         gTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGcggtacgg         <  2:544976/151‑1 (MQ=255)
                                                                                                                                         gTTCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGcggtac           <  4:308197/149‑1 (MQ=255)
                                                                                                                                          ttCGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGcggtacggt        >  4:1531228/1‑151 (MQ=255)
                                                                                                                                            cGTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCggtacggtac      <  3:1413/151‑1 (MQ=255)
                                                                                                                                             gTGCCCGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCggtacggtac      <  4:1601829/150‑1 (MQ=255)
                                                                                                                                                cccGCCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTAt  >  2:1105577/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTAGTTATCGGTTCTGTTCGTGCCCGTCTGGACAGACTTGATGAGCCGAGCAGCAACAAAGTGCTGCCAATCACCATCCACGGTGACGCCGCAGTGACCGGGCAGGGCGTGGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTAT  >  NC_000913/759568‑759862

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