Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,144,735 G→C 7.3% 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 7.3% 400.2 / 15.6 117A247A (GCG→GCCpflDputative glycine radical domain‑containing pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (51/57);  new base C (3/6);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 119 total aligned reads displayed.

CTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGC  >  NC_000913/4144588‑4144883
                                                                                                                                                   |                                                                                                                                                    
cTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGt                                                                                                                                                     <  2:105148/149‑1 (MQ=255)
 tCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTg                                                                                                                                                    >  4:44494/1‑149 (MQ=255)
   gCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGcc                                                                                                                                                  >  4:365652/1‑149 (MQ=255)
   gCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATAGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTgac                                                                                                                                                  >  4:565265/1‑147 (MQ=255)
    cAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCa                                                                                                                                                 <  3:96724/149‑1 (MQ=255)
    cAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCa                                                                                                                                                 <  3:281096/149‑1 (MQ=255)
       aaaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaa                                                                                                                                              <  3:838012/149‑3 (MQ=255)
        aaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaac                                                                                                                                             <  2:551779/149‑4 (MQ=255)
        aaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGAACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaac                                                                                                                                             <  1:249867/149‑4 (MQ=255)
         acacATTTTGCGTTACGCCGAACTGGCGGAAACGATGTCGGCAAACTGCACATATGCACAGAGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTAt                                                                                                                                            >  2:587189/1‑149 (MQ=255)
         acacATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTAt                                                                                                                                            >  2:99554/1‑149 (MQ=255)
           acaTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGt                                                                                                                                               <  3:553682/144‑1 (MQ=255)
           acaTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGt                                                                                                                                               >  4:553682/1‑144 (MQ=255)
           acaTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTc                                                                                                                                          >  4:251075/1‑149 (MQ=255)
             aTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTg                                                                                                                                        <  1:217670/149‑1 (MQ=255)
              ttttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTgg                                                                                                                                       >  4:53711/1‑149 (MQ=255)
                ttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTa                                                                                                                                     <  3:519585/149‑1 (MQ=255)
                ttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTa                                                                                                                                     <  1:99554/149‑1 (MQ=255)
                    gTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATg                                                                                                                                 >  1:491642/1‑149 (MQ=255)
                    gTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATg                                                                                                                                 >  1:165650/1‑149 (MQ=255)
                         gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcat                                                                                                                            >  4:469448/1‑149 (MQ=255)
                         gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcat                                                                                                                            >  2:482950/1‑149 (MQ=255)
                           cGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcatca                                                                                                                          <  3:749299/149‑1 (MQ=255)
                             aaCTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCAtt                                                                                                                        <  1:718338/149‑1 (MQ=255)
                             aaCTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCAtt                                                                                                                        <  3:656963/149‑1 (MQ=255)
                                   cGGAAACGATGGCGGAAAACTGCACAGATGCCAAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCaa                                                                                                                  <  2:895397/149‑1 (MQ=255)
                                     gAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATa                                                                                                                >  4:808406/1‑149 (MQ=255)
                                     gAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGACATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcatca                                                                                                                          >  1:106642/1‑139 (MQ=255)
                                     gAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTAGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcatca                                                                                                                          <  2:106642/139‑1 (MQ=255)
                                         cGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGaa                                                                                                            <  1:30083/149‑1 (MQ=255)
                                            tGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATcc                                                                                                         <  1:583091/149‑1 (MQ=255)
                                               cggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAAc                                                                                                      <  3:316842/149‑1 (MQ=255)
                                                ggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACg                                                                                                     >  1:694952/1‑149 (MQ=255)
                                                 gcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGc                                                                                                    <  1:482950/149‑1 (MQ=255)
                                                   aaaCTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAc                                                                                                                              >  2:308426/1‑121 (MQ=255)
                                                   aaaCTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAc                                                                                                                              <  1:308426/121‑1 (MQ=255)
                                                           caGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTa                                                                                          <  3:251075/149‑1 (MQ=255)
                                                           caGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTa                                                                                          <  4:652125/149‑1 (MQ=255)
                                                           caGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTACCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTa                                                                                          >  1:517784/1‑149 (MQ=255)
                                                            aGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTAt                                                                                         >  4:192701/1‑149 (MQ=255)
                                                            aGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTAt                                                                                         >  4:730361/1‑149 (MQ=255)
                                                             gATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATc                                                                                        <  3:469448/149‑1 (MQ=255)
                                                              aTGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTTCATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCg                                                                                       >  2:74701/1‑149 (MQ=255)
                                                                  ccAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                   >  4:844734/1‑149 (MQ=255)
                                                                  ccAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                   <  1:771827/149‑1 (MQ=255)
                                                                   cAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTggg                                                                                  <  1:666026/149‑1 (MQ=255)
                                                                         cgcgAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCtt                                                                            <  4:242279/149‑1 (MQ=255)
                                                                          gcgAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATGAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATGCGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTc                                                                           >  4:824123/1‑149 (MQ=255)
                                                                            gAAGTGCTGCTGACTATTGCAGCGATCTCCCGCCATGACGCGCTACATAAGCCGCAGACGTTCTGGAAGGCGTGCCCGGTATCCTGGTACATGATCATCATTCTTCAATACGAATCCAACGCCAGTTCGCTCTTGCTGGAGCGCTTCGa                                                                         >  1:48735/1‑149 (MQ=255)
                                                                             aaGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGAc                                                                        >  4:830701/1‑149 (MQ=255)
                                                                             aaGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGAc                                                                        <  3:853631/149‑1 (MQ=255)
                                                                               gagCTGCTGACTATTGCAGAGATCCCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCa                                                                      <  2:695844/149‑1 (MQ=255)
                                                                                agCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAg                                                                     >  3:108340/1‑149 (MQ=255)
                                                                                   tgctgACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCa                                                                                                  <  1:721120/117‑1 (MQ=255)
                                                                                   tgctgACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGtat                                                                  >  1:165029/1‑149 (MQ=255)
                                                                                   tgctgACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGtat                                                                  >  4:715457/1‑149 (MQ=255)
                                                                                    gctgACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCTCTATCGTTGGGGCGCTTCGACCAGtata                                                                 <  3:565265/149‑1 (MQ=255)
                                                                                     ctgACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTatat                                                                <  2:694952/149‑1 (MQ=255)
                                                                                       gACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGt                                                              >  4:534464/1‑149 (MQ=255)
                                                                                        aCTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGtt                                                             >  4:785947/1‑149 (MQ=255)
                                                                                         cTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTg                                                            <  3:684887/149‑1 (MQ=255)
                                                                                            ttGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCg                                                         >  2:102195/1‑149 (MQ=255)
                                                                                              gCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGAc                                                                        >  2:849371/1‑132 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTGTATGTTGCCGTTCt                                                     >  2:411903/1‑149 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                     >  3:515522/1‑149 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                     >  2:171955/1‑149 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                     <  1:347151/149‑1 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                     >  3:293115/1‑149 (MQ=255)
                                                                                                  agaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTAGGGCGCTTCGACCAGTATATGTTGCCGTTCTAt                                                   >  1:778381/1‑149 (MQ=255)
                                                                                                     tctcCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAg                                                <  3:785549/149‑1 (MQ=255)
                                                                                                        cccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTTCATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAca                                             <  1:74701/149‑1 (MQ=255)
                                                                                                         ccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGGATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcat                                            >  4:497848/1‑149 (MQ=255)
                                                                                                         ccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTAt                                                   <  1:340042/142‑1 (MQ=255)
                                                                                                         ccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTAt                                                   >  2:340042/1‑142 (MQ=255)
                                                                                                          cGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcatc                                           <  3:53711/149‑1 (MQ=255)
                                                                                                           gCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcatca                                          <  1:171955/149‑1 (MQ=255)
                                                                                                           gCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcatca                                          <  1:566150/149‑1 (MQ=255)
                                                                                                            ccATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcatcat                                         >  4:772606/1‑149 (MQ=255)
                                                                                                             cATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCAtt                                        >  2:620616/1‑149 (MQ=255)
                                                                                                               tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAGCGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                      <  2:533403/149‑1 (MQ=255)
                                                                                                               tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                      <  1:491026/149‑1 (MQ=255)
                                                                                                                aaCGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAAc                                     <  4:746910/149‑1 (MQ=255)
                                                                                                                 aCGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAAcc                                    >  1:621615/1‑149 (MQ=255)
                                                                                                                           tAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATAc                                                                                                               <  2:228540/64‑1 (MQ=255)
                                                                                                                             aGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGAt                        >  4:659740/1‑149 (MQ=255)
                                                                                                                               ccGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgc                      <  3:844734/149‑1 (MQ=255)
                                                                                                                                cGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcg                     >  1:166214/1‑149 (MQ=255)
                                                                                                                                  cAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCGGTTCGCGATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcggc                   >  1:640464/1‑149 (MQ=255)
                                                                                                                                  cAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcggc                   <  2:620080/149‑1 (MQ=255)
                                                                                                                                   aGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCACTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcggcg                  >  4:44695/1‑149 (MQ=255)
                                                                                                                                    gaataacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcg                     >  3:395754/7‑145 (MQ=255)
                                                                                                                                    gACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGt                 >  4:695428/1‑149 (MQ=255)
                                                                                                                                      cGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTc               <  1:859051/149‑1 (MQ=255)
                                                                                                                                       taacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTcc              >  2:21607/4‑149 (MQ=255)
                                                                                                                                          cTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGa           >  4:227822/1‑149 (MQ=255)
                                                                                                                                            ggCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa         <  2:469973/149‑1 (MQ=255)
                                                                                                                                             gCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAAc                                                                                                      >  3:841523/1‑55 (MQ=255)
                                                                                                                                             gCAGGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAAc                                                                                                      <  4:841523/55‑1 (MQ=255)
                                                                                                                                                   gTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGc  >  3:802709/1‑149 (MQ=255)
                                                                                                                                                   gTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGc  <  2:621615/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGC  >  NC_000913/4144588‑4144883

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: