Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,144,735 G→C 9.0% A247A (GCG→GCC)  pflD → putative glycine radical domain‑containing pyruvate formate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,144,7350GC 9.0% 295.4 / 16.2 89A247A (GCG→GCC) pflDputative glycine radical domain‑containing pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (32/48);  new base C (4/4);  total (37/52)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACT  >  NC_000913/4144588‑4144881
                                                                                                                                                   |                                                                                                                                                  
cTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGt                                                                                                                                                   >  1:374778/1‑149 (MQ=255)
cTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGt                                                                                                                                                   >  4:199091/1‑149 (MQ=255)
   gCAGAAACACATTTTGCGTTACGCCGAACTTGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGcc                                                                                                                                                <  1:449711/149‑1 (MQ=255)
       aaaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaa                                                                                                                                            <  3:642979/149‑3 (MQ=255)
       aaaCACATCTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaa                                                                                                                                            <  1:438694/149‑3 (MQ=255)
        aaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaac                                                                                                                                           <  1:250594/149‑4 (MQ=255)
          cacaTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTAtt                                                                                                                                         >  2:428412/1‑149 (MQ=255)
           acaTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTc                                                                                                                                        >  4:262618/1‑149 (MQ=255)
           acaTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTc                                                                                                                                        >  2:425062/1‑149 (MQ=255)
              ttttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTgg                                                                                                                                     <  4:237734/149‑1 (MQ=255)
              ttttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTgg                                                                                                                                     <  4:346449/149‑1 (MQ=255)
                ttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTa                                                                                                                                   <  3:139042/149‑1 (MQ=255)
                  gCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACa                                                                                                                                 <  1:689780/149‑1 (MQ=255)
                   cGTTACGCTGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGGGCAGAGCTGCTGACTAGTGCAGAGATCTCCCGCCATAACGCGCAACATAAGTCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACAt                                                                                                                                <  2:469457/149‑1 (MQ=255)
                    gTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATg                                                                                                                               >  2:603167/1‑149 (MQ=255)
                      tACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGTAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGaa                                                                                                                             <  4:401226/149‑1 (MQ=255)
                       aCGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCTAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAc                                                                                                                            <  4:65839/149‑1 (MQ=255)
                         gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaacgt                                                                                                                                         <  3:83252/134‑6 (MQ=255)
                         gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaacgt                                                                                                                                         >  4:83252/1‑129 (MQ=255)
                                 ggCGGAAACGATGGCGGCAAACTGCACAGCTGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGc                                                                                                                  <  3:127601/149‑1 (MQ=255)
                                   cGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCaa                                                                                                                <  1:253036/149‑1 (MQ=255)
                                   cGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCaa                                                                                                                <  1:425062/149‑1 (MQ=255)
                                   cGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCAGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCaa                                                                                                                <  2:593851/149‑1 (MQ=255)
                                              gcggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCaa                                                                                                     <  4:629901/149‑1 (MQ=255)
                                              gcggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCaa                                                                                                     <  3:304600/149‑1 (MQ=255)
                                                 gcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGc                                                                                                  <  4:171315/149‑1 (MQ=255)
                                                       tGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTg                                                                                                                                      >  2:369423/1‑107 (MQ=255)
                                                       tGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTg                                                                                                                                      <  1:369423/107‑1 (MQ=255)
                                                       tGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTc                                                                                            <  4:62162/149‑1 (MQ=255)
                                                       tGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTc                                                                                            <  3:589220/149‑1 (MQ=255)
                                                       tGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTc                                                                                            >  4:265172/1‑149 (MQ=255)
                                                       tGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTc                                                                                            >  4:363083/1‑149 (MQ=255)
                                                          acaGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATTACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCt                                                                                         >  1:228915/1‑149 (MQ=255)
                                                          acaGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACa                                                                                                                                 <  4:362968/109‑1 (MQ=255)
                                                          acaGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACa                                                                                                                                 >  3:362968/1‑109 (MQ=255)
                                                                gCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGtt                                                                                   <  4:56705/149‑1 (MQ=255)
                                                                  ccAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTAc                                                                                                                                  <  4:204711/100‑1 (MQ=255)
                                                                  ccAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTAc                                                                                                                                  >  3:204711/1‑100 (MQ=255)
                                                                  ccAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                 <  4:629335/149‑1 (MQ=255)
                                                                        tCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCt                                                                           <  4:412000/149‑1 (MQ=255)
                                                                        tCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCt                                                                           <  1:159062/149‑1 (MQ=255)
                                                                            gAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGa                                                                       <  1:603167/149‑1 (MQ=255)
                                                                            gAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGa                                                                       >  1:502842/1‑149 (MQ=255)
                                                                               gagCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCa                                                                    <  4:253380/149‑1 (MQ=255)
                                                                                   tgctgACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGtat                                                                >  2:651850/1‑149 (MQ=255)
                                                                                    gctgACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGtata                                                               <  1:526278/149‑1 (MQ=255)
                                                                                       gACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGt                                                            <  3:262618/149‑1 (MQ=255)
                                                                                        aCTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGtt                                                           <  3:280619/149‑1 (MQ=255)
                                                                                          tATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATcc                                                                                                       <  4:615389/103‑1 (MQ=255)
                                                                                          tATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGc                                                         >  4:505969/1‑149 (MQ=255)
                                                                                          tATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGCTATTCTGGTACATGAACATCATTCTGCAATACGAATcc                                                                                                       >  3:615389/1‑103 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                   >  3:441838/1‑149 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                   >  3:648181/1‑149 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                   >  4:302867/1‑149 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                   <  1:513507/149‑1 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                   <  3:587677/149‑1 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                   >  1:16421/1‑149 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                   >  2:670219/1‑149 (MQ=255)
                                                                                                    aTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                   <  2:448524/145‑1 (MQ=255)
                                                                                                    aTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCa                                               >  2:36344/1‑149 (MQ=255)
                                                                                                    aTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCa                                               <  3:265172/149‑1 (MQ=255)
                                                                                                    aTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGGGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                   >  1:448524/1‑145 (MQ=255)
                                                                                                        cccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAca                                           <  4:153365/149‑1 (MQ=255)
                                                                                                         ccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTAt                                                 <  1:453982/142‑1 (MQ=255)
                                                                                                         ccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTAt                                                 >  2:453982/1‑142 (MQ=255)
                                                                                                         ccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcat                                          >  3:306178/1‑149 (MQ=255)
                                                                                                              aTAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTa                                     >  3:213676/1‑149 (MQ=255)
                                                                                                               tAACGCGCAACATCAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                    <  1:654127/149‑1 (MQ=255)
                                                                                                                      cAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAggg                             <  2:374778/149‑1 (MQ=255)
                                                                                                                       aaCATAAGCCGCAGACGTTCTGGCAGGCATGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGc                            >  1:371074/1‑149 (MQ=255)
                                                                                                                          aTAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcat                                          >  1:53609/1‑132 (MQ=255)
                                                                                                                          aTAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcat                                          <  2:53609/132‑1 (MQ=255)
                                                                                                                          aTAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGaa                         <  3:505969/149‑1 (MQ=255)
                                                                                                                           tAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAg                        >  1:385304/1‑149 (MQ=255)
                                                                                                                            aaGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAggg                             >  3:514901/1‑143 (MQ=255)
                                                                                                                            aaGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAggg                             <  4:514901/143‑1 (MQ=255)
                                                                                                                               ccGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGATCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgc                    <  3:199091/149‑1 (MQ=255)
                                                                                                                                cGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcg                   <  1:384183/149‑1 (MQ=255)
                                                                                                                                 gCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                 <  2:631673/86‑1 (MQ=255)
                                                                                                                                 gCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                 >  1:631673/1‑86 (MQ=255)
                                                                                                                                 gCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcgg                  <  2:451102/149‑1 (MQ=255)
                                                                                                                                  cAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcggc                 >  2:504091/1‑149 (MQ=255)
                                                                                                                                     aCGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGtt              >  2:204643/1‑149 (MQ=255)
                                                                                                                                      ataacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCGGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTc             >  1:255004/5‑149 (MQ=255)
                                                                                                                                      ataacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTc             >  3:627239/5‑149 (MQ=255)
                                                                                                                                          cTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGa         >  3:311954/1‑149 (MQ=255)
                                                                                                                                            ggCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa       <  3:302867/149‑1 (MQ=255)
                                                                                                                                               aGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGaa    <  2:16421/149‑1 (MQ=255)
                                                                                                                                                 gCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATTAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACt  <  2:228915/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACT  >  NC_000913/4144588‑4144881

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑