Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,084,656 T→G E419A (GAA→GCA)  speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,084,6560TG98.6% 264.7 / ‑4.4 69E419A (GAA→GCA) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (39/29);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GGTTTTGCGGATCCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCGATGATATTAGACACCAGCACG  >  NC_000913/3084510‑3084802
                                                                                                                                                  |                                                                                                                                                  
ggTTTTGCGGATCCAGCCGCTTTTGCACTTCATGGCACATGCTCAACTAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCAc                                                                                                                                                 >  1:812368/1‑150 (MQ=255)
  ttttGCGGATCCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGc                                                                                                                                               >  1:474951/1‑150 (MQ=255)
       cGGATCCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAgaga                                                                                                                                          >  1:783124/1‑150 (MQ=255)
       cGGATCCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAgaga                                                                                                                                          >  2:1187420/1‑150 (MQ=255)
       cGGATCCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGaa                                                                                                                                         <  2:1135755/151‑1 (MQ=255)
         gATCCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACg                                                                                                                                       <  1:220316/151‑1 (MQ=255)
          aTCCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACg                                                                                                                                       <  1:869595/150‑1 (MQ=255)
                 tgctTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTc                                                                                                                               >  2:274071/1‑151 (MQ=255)
                     tttGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTccc                                                                                                                             <  1:1207854/149‑1 (MQ=255)
                     tttGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTccc                                                                                                                             <  2:211644/149‑1 (MQ=255)
                      ttGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGt                                                                                                                          >  2:786476/1‑151 (MQ=255)
                         cACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTc                                                                                                                        >  1:731836/1‑150 (MQ=255)
                         cACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCg                                                                                                                       >  2:656964/1‑151 (MQ=255)
                             tCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTg                                                                                                                     <  1:847220/149‑1 (MQ=255)
                                     caTGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGc                                                                                                           >  2:966729/1‑151 (MQ=255)
                                     caTGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGc                                                                                                           >  1:811459/1‑151 (MQ=255)
                                     caTGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGc                                                                                                           >  2:827216/1‑151 (MQ=255)
                                        gCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAg                                                                                                        <  2:468624/151‑1 (MQ=255)
                                        gCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGcc                                                                                                          >  2:256401/1‑149 (MQ=255)
                                         cTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAgg                                                                                                       >  1:745991/1‑151 (MQ=255)
                                           cAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCCTCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCCCCTGCCAgg                                                                                                       >  1:743007/1‑149 (MQ=255)
                                             aaTAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTc                                                                                                   <  1:1235245/151‑1 (MQ=255)
                                               tAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTccc                                                                                                 >  1:642648/1‑151 (MQ=255)
                                                 aaGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCcac                                                                                               >  2:1244720/1‑151 (MQ=255)
                                                      gctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATg                                                                                            <  1:47133/149‑1 (MQ=255)
                                                      gctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGc                                                                                           <  1:1155885/150‑1 (MQ=255)
                                                        tCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCttt                                                                                        >  1:1097803/1‑151 (MQ=255)
                                                           gCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGc                                                                                      >  2:454731/1‑150 (MQ=255)
                                                             ccATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAgc                                                                                   <  1:1083001/151‑1 (MQ=255)
                                                             ccATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAgc                                                                                   <  1:247656/151‑1 (MQ=255)
                                                              cATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAgcg                                                                                  <  2:2185/151‑1 (MQ=255)
                                                                tGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAgcgc                                                                                 <  2:921145/150‑1 (MQ=255)
                                                                 gCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAgcgcg                                                                                <  2:639031/150‑1 (MQ=255)
                                                                     gTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCgcgg                                                                           <  1:1083162/151‑1 (MQ=255)
                                                                          tGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCAt                                                                      >  1:445800/1‑151 (MQ=255)
                                                                             aGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCAtctt                                                                   >  1:986247/1‑151 (MQ=255)
                                                                              ggCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCAtcttc                                                                  <  1:139129/151‑1 (MQ=255)
                                                                                   aaGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAg                                                             <  2:880567/151‑1 (MQ=255)
                                                                                      aTGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGgcg                                                          >  2:751117/1‑151 (MQ=255)
                                                                                        gCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGcgcg                                                        <  1:228376/151‑1 (MQ=255)
                                                                                           ggAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTc                                                     <  2:55979/151‑1 (MQ=255)
                                                                                            gAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCg                                                    <  1:793783/151‑1 (MQ=255)
                                                                                            gAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCg                                                    <  2:517611/151‑1 (MQ=255)
                                                                                             aaGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCgg                                                   >  2:1014934/1‑151 (MQ=255)
                                                                                              agagTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGc                                                  >  1:1262666/1‑151 (MQ=255)
                                                                                               gagTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCa                                                 <  2:1068715/151‑1 (MQ=255)
                                                                                                    gCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCg                                              >  2:1258699/1‑149 (MQ=255)
                                                                                                     ccGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCgtgt                                           >  1:1154767/1‑151 (MQ=255)
                                                                                                      cGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCgtg                                            >  1:29718/1‑149 (MQ=255)
                                                                                                               tGTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACg                                 >  1:1291313/1‑151 (MQ=255)
                                                                                                                gTCGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTaag                                 <  2:417806/150‑3 (MQ=255)
                                                                                                                  cGTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTc                              >  1:507560/1‑151 (MQ=255)
                                                                                                                   gTGCAGATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTcc                             <  2:411008/151‑1 (MQ=255)
                                                                                                                       aGATCCATCTGACTGTCGTTTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACg                          <  1:623945/150‑1 (MQ=255)
                                                                                                                        gATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGcc                        >  1:515434/1‑151 (MQ=255)
                                                                                                                        gATCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGcc                        <  2:57718/151‑1 (MQ=255)
                                                                                                                         aTCCATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGGGGATTAGTTACGTTCCACGCCg                       >  1:729599/1‑151 (MQ=255)
                                                                                                                            cATCTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCgatg                    >  2:1047269/1‑151 (MQ=255)
                                                                                                                               cTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTGCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGAATTCGTTACGTTCCACGCCgatga                   >  1:423561/1‑149 (MQ=255)
                                                                                                                               cTGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCgatgat                  >  1:1080581/1‑150 (MQ=255)
                                                                                                                                tGACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCGATGatat                >  2:1103560/1‑151 (MQ=255)
                                                                                                                                 gACTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCgatga                   <  1:1222455/147‑1 (MQ=255)
                                                                                                                                  aCTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCGATGATATTa              <  1:514671/151‑1 (MQ=255)
                                                                                                                                   cTGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCGATGATATTa              <  2:278290/150‑1 (MQ=255)
                                                                                                                                    tGTCGTGTAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCGATGATAtt               >  1:822978/1‑148 (MQ=255)
                                                                                                                                        gtgtAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTACTCGTTACGTTCCACGCCGATGATATTAGACAcc        >  1:845049/1‑151 (MQ=255)
                                                                                                                                         tgtAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCGATGATATTAGACAcc        >  1:962032/1‑150 (MQ=255)
                                                                                                                                           tAACCATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCGATGATATTAGACACCAGc     >  2:27208/1‑151 (MQ=255)
                                                                                                                                              ccATGCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCGATGATATTAGACACCAGCACg  >  2:727542/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGTTTTGCGGATCCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGGTTTCCCACATGCTTTGCAGCGCGCGCGGCGCATCTTCTGCAGGCGCGGTCGGCACCGTGTATTCGTTACGTTCCACGCCGATGATATTAGACACCAGCACG  >  NC_000913/3084510‑3084802

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