Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,125,147 T→G K105Q (AAA→CAA)  glcF ← glycolate oxidase 4Fe‑4S iron‑sulfur cluster subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,125,1470TG98.7% 300.4 / ‑3.5 79K105Q (AAA→CAA) glcFglycolate oxidase 4Fe‑4S iron‑sulfur cluster subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (36/42);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

TCAGCAGGCAGTTTTGCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTTCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTG  >  NC_000913/3124999‑3125291
                                                                                                                                                    |                                                                                                                                                
tCAGCAGGCAGTTTTGCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCt                                                                                                                                                >  2:1158620/1‑151 (MQ=255)
 cagcagGCAGTTTTGCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTg                                                                                                                                               <  1:1188237/151‑1 (MQ=255)
   gcaggcagCTTTGCTTTGATCTGTTCCGGGAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCCCTACCTCGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTgctgct                                                                                                                                             >  2:551656/1‑151 (MQ=255)
     aggcagTTTTGCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCg                                                                                                                                           >  1:1198523/1‑151 (MQ=255)
        cagTTTTGCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGAc                                                                                                                                         >  1:416974/1‑150 (MQ=255)
        cactgttGCTCTGACCTGATCCGGTAAAAACGGTCGCAGCACCAGCCCTGCTTGCGTCGGCGCACGGAAGCCCGCCGGACGCGGCTCTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTc                                                                                                                                            <  2:362274/142‑1 (MQ=255)
         agTTTTGCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACa                                                                                                                                        >  2:471439/1‑150 (MQ=255)
             ttGCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAatat                                                                                                                                    <  1:405235/150‑1 (MQ=255)
             ttGCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAatat                                                                                                                                    <  2:596679/150‑1 (MQ=255)
               gCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCa                                                                                                                                  <  2:967303/150‑1 (MQ=255)
               gCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCAc                                                                                                                                 <  1:888725/151‑1 (MQ=255)
                    gACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCg                                                                                                                            >  1:202279/1‑151 (MQ=255)
                    gACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCg                                                                                                                            <  2:207176/151‑1 (MQ=255)
                      ccTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGat                                                                                                                          <  1:186957/151‑1 (MQ=255)
                      ccTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCACCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGat                                                                                                                          >  1:145124/1‑151 (MQ=255)
                         gTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGata                                                                                                                         <  2:199812/149‑1 (MQ=255)
                          ttCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATcc                                                                                                                      >  2:879116/1‑151 (MQ=255)
                             cGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCa                                                                                                                     <  1:557908/149‑1 (MQ=255)
                             cGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGc                                                                                                                   <  1:989464/151‑1 (MQ=255)
                                 aaaaaCGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAAt                                                                                                               >  1:995538/1‑151 (MQ=255)
                                  aaaaCGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAAtt                                                                                                              <  2:1198201/151‑1 (MQ=255)
                                  aaaaCGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAAtt                                                                                                              <  2:236970/151‑1 (MQ=255)
                                   aaaCGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATtg                                                                                                             <  2:809716/151‑1 (MQ=255)
                                    aaCGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATtgt                                                                                                            <  1:495131/151‑1 (MQ=255)
                                      cGGTCGCAGCACCAGCCCTGCCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTGCCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGa                                                                                                          <  2:1202121/151‑1 (MQ=255)
                                      cGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATtgtg                                                                                                           >  2:479921/1‑150 (MQ=255)
                                      cGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATtgt                                                                                                            <  2:618409/149‑1 (MQ=255)
                                               caCCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGcaca                                                                                                 <  1:997811/151‑1 (MQ=255)
                                               caCCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGcaca                                                                                                 <  1:678030/151‑1 (MQ=255)
                                                      ccTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGa                                                                                            >  1:975385/1‑149 (MQ=255)
                                                      ccTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAACCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGa                                                                                            >  2:734276/1‑149 (MQ=255)
                                                            tGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAgg                                                                                    >  1:917879/1‑151 (MQ=255)
                                                             gCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAggt                                                                                   <  1:134944/151‑1 (MQ=255)
                                                               gTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAggtgg                                                                                 >  2:29722/1‑151 (MQ=255)
                                                                 cAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAggtggt                                                                                <  2:528066/150‑1 (MQ=255)
                                                                 cAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGtt                                                                               <  1:159141/151‑1 (MQ=255)
                                                                   gcgcACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTc                                                                             <  1:61845/151‑1 (MQ=255)
                                                                     gcACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTcac                                                                           >  2:1184561/1‑151 (MQ=255)
                                                                      cACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTcaca                                                                          <  1:557458/151‑1 (MQ=255)
                                                                           aaGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTAc                                                                     >  1:363342/1‑151 (MQ=255)
                                                                               ccgccgGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCa                                                                 <  1:176863/151‑1 (MQ=255)
                                                                                 gccgGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCa                                                                 >  2:507069/1‑149 (MQ=255)
                                                                                 gccgGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGAGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGAGGTTTCACAATTACGGc                                                                  >  2:665333/1‑148 (MQ=255)
                                                                                    ggACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCaa                                                                >  2:704684/1‑147 (MQ=255)
                                                                                    ggACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGa                                                            <  1:610536/151‑1 (MQ=255)
                                                                                    ggACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGa                                                            <  2:194176/151‑1 (MQ=255)
                                                                                     gACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAg                                                           <  2:47799/151‑1 (MQ=255)
                                                                                        gcgGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGc                                                         >  1:628766/1‑150 (MQ=255)
                                                                                        gcgGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCa                                                        <  1:447632/151‑1 (MQ=255)
                                                                                            cACTACCTGGCGCAATCCTTCGCGCAGTAGTCGCTCCGGCAGTGGGCGTTTCTCTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCTCACCATAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGa                                                    <  1:778152/151‑1 (MQ=255)
                                                                                             aCTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGAt                                                   <  1:929311/151‑1 (MQ=255)
                                                                                             aCTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGAt                                                   >  1:1195121/1‑151 (MQ=255)
                                                                                             aCTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGAt                                                   >  2:71973/1‑151 (MQ=255)
                                                                                             aCTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGAt                                                   >  1:277332/1‑151 (MQ=255)
                                                                                               tACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCg                                                 >  1:897466/1‑151 (MQ=255)
                                                                                               tACCTGGCGCAATCCTTCGCGCAGTATTCGCCCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCg                                                 >  1:627185/1‑151 (MQ=255)
                                                                                                aCCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCga                                                <  1:74327/151‑1 (MQ=255)
                                                                                                aCCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCg                                                 >  1:208032/1‑150 (MQ=255)
                                                                                                 ccTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCgag                                               <  1:725384/151‑1 (MQ=255)
                                                                                                         aaTCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCt                                       >  1:50647/1‑151 (MQ=255)
                                                                                                              ttCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGCTATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGtt                                  >  2:826177/1‑151 (MQ=255)
                                                                                                              ttCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGtt                                  >  1:534924/1‑151 (MQ=255)
                                                                                                               tCGCGCAGTACTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGtt                                  <  2:303516/150‑1 (MQ=255)
                                                                                                                cgcgCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGtttt                                <  2:726576/151‑1 (MQ=255)
                                                                                                                 gcgcAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGtttt                                >  1:434198/1‑150 (MQ=255)
                                                                                                                   gcAGTATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAg                             <  1:606175/151‑1 (MQ=255)
                                                                                                                       tATTCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTg                         >  1:913001/1‑151 (MQ=255)
                                                                                                                        agtCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTg                         <  2:642976/148‑1 (MQ=255)
                                                                                                                          tCGCTCCGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACt                      >  2:704266/1‑151 (MQ=255)
                                                                                                                             cgccGGCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCg                   <  2:100585/149‑1 (MQ=255)
                                                                                                                                  gCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGcc              >  2:331727/1‑151 (MQ=255)
                                                                                                                                  gCAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGCGACTTCGTTGcc              >  2:955262/1‑151 (MQ=255)
                                                                                                                                   cAGTGGGCGTTTCACTTTCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCt             >  2:255761/1‑151 (MQ=255)
                                                                                                                                   cAGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCt             <  2:284821/151‑1 (MQ=255)
                                                                                                                                    aGTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCtt            >  2:1078564/1‑151 (MQ=255)
                                                                                                                                     gTGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTc           >  1:391680/1‑151 (MQ=255)
                                                                                                                                      tGGGCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTcc          >  1:905088/1‑151 (MQ=255)
                                                                                                                                       gggCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTcc          <  2:145117/150‑1 (MQ=255)
                                                                                                                                         gCGTTTCACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGc       <  2:955150/151‑1 (MQ=255)
                                                                                                                                               cACTTGCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTg  <  1:286603/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCAGCAGGCAGTTTTGCTCTGACCTGTTCCGGTAAAAACGGTCGCAGCACCAGCCCTACCTGCGTCAGCGCACGGAAGACCGCCGGACGCGGCACTACCTGGCGCAATCCTTCGCGCAGTATTCGCTCCGGCAGTGGGCGTTTCACTTTCTGCTCGACAATATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTG  >  NC_000913/3124999‑3125291

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