Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,579,154 A→C 5.0% A163A (GCA→GCC)  talA → transaldolase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,579,1540AC 5.0% 304.8 / 3.5 98A163A (GCA→GCC) talAtransaldolase A
Reads supporting (aligned to +/- strand):  ref base A (43/50);  new base C (3/2);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTACTA  >  NC_000913/2579008‑2579300
                                                                                                                                                  |                                                                                                                                                  
gtggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACg                                                                                                                                                  <  4:348156/149‑1 (MQ=255)
 tggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACgg                                                                                                                                                 <  8:223233/149‑1 (MQ=255)
 tggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACgg                                                                                                                                                 <  4:34360/149‑1 (MQ=255)
   gACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGc                                                                                                                                               >  4:67174/1‑149 (MQ=255)
     cTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCt                                                                                                                                             >  1:90464/1‑149 (MQ=255)
           tCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCg                                                                                                                                        >  7:142154/1‑148 (MQ=255)
             aGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCgg                                                                                                                                       >  5:110107/1‑147 (MQ=255)
             aGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGaa                                                                                                                                     <  6:227304/149‑1 (MQ=255)
             aGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGCATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGaa                                                                                                                                     <  3:289132/149‑1 (MQ=255)
              gCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAg                                                                                                                                    <  2:296610/149‑1 (MQ=255)
              gCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAg                                                                                                                                    <  2:290419/149‑1 (MQ=255)
              gCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTCGCACAGGCCCGTGCCTGTGcaaaa                                                                                                                                     <  5:259384/148‑5 (MQ=255)
              gCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTCGCACAGGCCCGTGCCTGTGcaa                                                                                                                                       >  6:259384/1‑144 (MQ=255)
                 acaaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAg                                                                                                                                 <  6:397401/149‑1 (MQ=255)
                  caaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAgg                                                                                                                                >  2:239077/1‑149 (MQ=255)
                   aaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGc                                                                                                                               >  8:283139/1‑149 (MQ=255)
                    aGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCg                                                                                                                              >  2:275779/1‑149 (MQ=255)
                       cGTTGAGAAATCTCGCATTCTGATCAAGCTGGCTTCGAGATGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAATGTATTTACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGTCACGGGCCTGTGCGGAAGCAGGCGttt                                                                                                                           <  8:231617/149‑1 (MQ=255)
                       cGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGttt                                                                                                                           <  4:56511/149‑1 (MQ=255)
                         ttGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGttttt                                                                                                                         <  2:287938/149‑1 (MQ=255)
                            agaAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCt                                                                                                                       >  3:401888/1‑148 (MQ=255)
                            agaAATATCGCATTCTGATCAAGCTGGCTTCGACCTGGGATGGATTTCGCTCGGCAGAAGAGCTGGAAAAAGAAGGTCTTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTg                                                                                                                      <  4:309314/149‑1 (MQ=255)
                             gaAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGa                                                                                                                     >  8:115535/1‑149 (MQ=255)
                                aTCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGAttt                                                                                                                  >  2:277087/1‑149 (MQ=255)
                                    cGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGcc                                                                                                              <  7:174213/149‑1 (MQ=255)
                                     gCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCg                                                                                                             >  3:95193/1‑149 (MQ=255)
                                     gCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCg                                                                                                             <  6:107276/149‑1 (MQ=255)
                                       aTTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGtt                                                                                                           >  3:227670/1‑149 (MQ=255)
                                        ttCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGttt                                                                                                          <  2:402280/149‑1 (MQ=255)
                                         tCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTg                                                                                                         >  1:338819/1‑149 (MQ=255)
                                           tGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaaa                                                                                                                                     <  5:286772/119‑5 (MQ=255)
                                           tGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaaa                                                                                                                                     >  6:286772/1‑115 (MQ=255)
                                             aTCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCgg                                                                                                     <  1:275779/149‑1 (MQ=255)
                                                aaGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGc                                                                                                   <  4:401888/148‑1 (MQ=255)
                                                    tGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTAtt                                                                                              >  5:136745/1‑149 (MQ=255)
                                                    tGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTAtt                                                                                              >  4:72700/1‑149 (MQ=255)
                                                      gCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTa                                                                                            <  3:67174/149‑1 (MQ=255)
                                                           gACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGAc                                                                                        <  1:190890/148‑1 (MQ=255)
                                                            aCCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTg                                                                                      <  3:214230/149‑1 (MQ=255)
                                                            aCCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTg                                                                                      <  6:57131/149‑1 (MQ=255)
                                                              cTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTgg                                                                                     >  6:234094/1‑148 (MQ=255)
                                                                gggAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTAt                                                                                  >  7:156993/1‑149 (MQ=255)
                                                                  gaaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAGGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTAt                                                                                  >  8:221179/1‑147 (MQ=255)
                                                                  gaaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCATTTGCACATGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTAGGGCGTATTTATGACTTGAATCa                                                                                >  2:131397/1‑149 (MQ=255)
                                                                   aaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTAt                                                                                           <  7:283767/137‑1 (MQ=255)
                                                                   aaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTAt                                                                                           >  8:283767/1‑137 (MQ=255)
                                                                   aaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAg                                                                               <  6:110107/149‑1 (MQ=255)
                                                                    aggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGACGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAGGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAgg                                                                              <  7:221179/149‑1 (MQ=255)
                                                                        aTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTg                                                                                      >  6:73135/1‑137 (MQ=255)
                                                                        aTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTg                                                                                      <  5:73135/137‑1 (MQ=255)
                                                                        aTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACg                                                                          >  1:43642/1‑149 (MQ=255)
                                                                         ttCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGAACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACg                                                                          >  5:129104/1‑148 (MQ=255)
                                                                            gcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCaa                                                                       <  1:277087/148‑1 (MQ=255)
                                                                            gcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAg                                                                      >  7:416629/1‑149 (MQ=255)
                                                                                 gCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGa                                                                  >  8:424912/1‑148 (MQ=255)
                                                                                 gCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGa                                                                  <  7:424912/148‑1 (MQ=255)
                                                                                 gCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGAt                                                                 <  4:95193/149‑1 (MQ=255)
                                                                                  cAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATg                                                                <  8:156993/149‑1 (MQ=255)
                                                                                      agagCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAcc                                                            >  6:56968/1‑149 (MQ=255)
                                                                                       gagCTGGAAAAAGAAGGTATAAACTGCAACCTGACGCGGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAccc                                                           <  8:405263/149‑1 (MQ=255)
                                                                                             gAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAtgt                                                     >  5:369557/1‑149 (MQ=255)
                                                                                             gAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAtgt                                                     >  7:75913/1‑149 (MQ=255)
                                                                                               aaaaGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtgg                                                   <  3:72700/149‑1 (MQ=255)
                                                                                                 aagaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGc                                                                             <  1:311243/121‑1 (MQ=255)
                                                                                                 aagaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGc                                                                             >  2:311243/1‑121 (MQ=255)
                                                                                                  agaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtggtg                                                 >  3:311142/1‑148 (MQ=255)
                                                                                                     agGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGtggtgg                                                >  6:54548/1‑146 (MQ=255)
                                                                                                      ggTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaaga                                            <  5:56968/149‑1 (MQ=255)
                                                                                                          ttAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaag                                             >  5:426038/1‑144 (MQ=255)
                                                                                                          ttAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAagaaga                                         <  2:338819/148‑1 (MQ=255)
                                                                                                           tAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGaagaag                                          <  2:156774/146‑1 (MQ=255)
                                                                                                           tAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGaagaag                                          >  1:156774/1‑146 (MQ=255)
                                                                                                             aCTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCg                                     <  7:115535/149‑1 (MQ=255)
                                                                                                              cTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCgg                                    >  5:369005/1‑149 (MQ=255)
                                                                                                              cTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTAGAAGAAGATCCgg                                    <  4:227670/149‑1 (MQ=255)
                                                                                                               tGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCggg                                   >  5:260995/1‑149 (MQ=255)
                                                                                                                 cAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGc                                  <  6:369005/148‑1 (MQ=255)
                                                                                                                  aaCCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGt                                <  2:90464/149‑1 (MQ=255)
                                                                                                                        tcGCTGCTGTTTTCTTTTGCACCGGCACGGGCCTGTGCGGAAGCAGGCGGTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAAt                           <  1:131397/147‑1 (MQ=255)
                                                                                                                         cGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCg                         <  7:283139/149‑1 (MQ=255)
                                                                                                                         cGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCg                         <  5:54548/149‑1 (MQ=255)
                                                                                                                          gctgctGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaaga                                            >  6:84488/1‑129 (MQ=255)
                                                                                                                          gctgctGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaaga                                            <  5:84488/129‑1 (MQ=255)
                                                                                                                           ctgctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGt                       <  6:369557/149‑1 (MQ=255)
                                                                                                                              ctgTTTTCTTTTGCCCAGGCCCGGGCCTGTGCGGAAGCAGGCGTTTTTCTGAGTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAGCCGTATGGGGTGGAAGAGGATCCGGGCGTTAAATCGGTgcg                    >  5:299403/1‑149 (MQ=255)
                                                                                                                                 ttttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCgg                        >  8:175257/1‑142 (MQ=255)
                                                                                                                                 ttttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCgg                        <  7:175257/142‑1 (MQ=255)
                                                                                                                                 ttttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCaa                 <  6:426038/149‑1 (MQ=255)
                                                                                                                                  tttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAccc                                                           >  4:102221/1‑106 (MQ=255)
                                                                                                                                  tttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAccc                                                           <  3:102221/106‑1 (MQ=255)
                                                                                                                                      ttttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCaa                 <  1:109165/144‑1 (MQ=255)
                                                                                                                                      ttttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCaa                 >  2:109165/1‑144 (MQ=255)
                                                                                                                                        ttGCACAGGCACGGGCCTGGGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTGGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATTTGGTGGAAGAAGATCCGGGCGTTATATCGGTGCGCAATATc             <  4:194120/146‑1 (MQ=255)
                                                                                                                                         tGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTAc          >  5:183077/1‑148 (MQ=255)
                                                                                                                                            acaGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGac       >  3:303927/1‑148 (MQ=255)
                                                                                                                                             caGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGacta     >  5:320976/1‑149 (MQ=255)
                                                                                                                                             caGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGGGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGacta     <  6:129104/149‑1 (MQ=255)
                                                                                                                                                gCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGActacta  >  7:291870/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTACTA  >  NC_000913/2579008‑2579300

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