Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,144,732 G→C 5.6% Q246H (CAG→CAC 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,7320GC 5.6% 550.7 / 16.0 154Q246H (CAG→CACpflDputative glycine radical domain‑containing pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (66/79);  new base C (4/5);  total (70/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

AAGCCTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAAC  >  NC_000913/4144584‑4144880
                                                                                                                                                    |                                                                                                                                                    
aaGCCTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAg                                                                                                                                                      <  1:121445/149‑1 (MQ=255)
 aGCCTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAgg                                                                                                                                                     >  3:894560/1‑149 (MQ=255)
     tCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTg                                                                                                                                                 >  3:912098/1‑149 (MQ=255)
         aGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAg                                                                                                                                             <  3:1713744/149‑1 (MQ=255)
           aaaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGtt                                                                                                                                           >  1:699348/1‑149 (MQ=255)
             acacATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGTCAGTTAt                                                                                                                                         >  1:1179554/1‑149 (MQ=255)
              cacaTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTACTGGCACGCCTGCCAgaacgt                                                                                                                                        <  4:1429766/149‑6 (MQ=255)
                 aTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTg                                                                                                                                     <  1:1688465/149‑1 (MQ=255)
                  ttttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTgg                                                                                                                                    <  1:850024/149‑1 (MQ=255)
                    ttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTa                                                                                                                                  <  3:212621/149‑1 (MQ=255)
                       cGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACAt                                                                                                                               <  4:326682/149‑1 (MQ=255)
                        gTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATg                                                                                                                              >  2:273364/1‑149 (MQ=255)
                         ttACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGa                                                                                                                             >  4:1520658/1‑149 (MQ=255)
                           aCGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAc                                                                                                                           <  3:1485516/149‑1 (MQ=255)
                             gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcat                                                                                                                         >  3:623762/1‑149 (MQ=255)
                             gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcat                                                                                                                         >  2:374627/1‑149 (MQ=255)
                               cGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcatca                                                                                                                       <  4:894560/149‑1 (MQ=255)
                                 aaCTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGTCAGTTATTCTGGTACATGAACATCAtt                                                                                                                     <  2:1179554/149‑1 (MQ=255)
                                  aCTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTc                                                                                                                    <  2:682920/149‑1 (MQ=255)
                                   cTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCGGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCt                                                                                                                   <  4:928527/149‑1 (MQ=255)
                                       cagAAACGATGGCGGCAAACTGCACAGCTCCCCAGCGTCGCGAAGAGCTGCTGCTTATTGCAGAGATCCCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCaa                                                                                                               <  4:367567/147‑1 (MQ=255)
                                       cGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCaa                                                                                                               <  1:620934/149‑1 (MQ=255)
                                             cGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGaa                                                                                                         <  2:1473140/149‑1 (MQ=255)
                                               aTGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATc                                                                                                       >  1:796265/1‑149 (MQ=255)
                                               aTGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATc                                                                                                       >  4:178769/1‑149 (MQ=255)
                                                 ggcggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCa                                                                                                     >  2:1097246/1‑149 (MQ=255)
                                                  gcggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCaa                                                                                                    <  4:568266/149‑1 (MQ=255)
                                                   cggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAAc                                                                                                   <  2:1689685/149‑1 (MQ=255)
                                                    ggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcatca                                                                                                                       <  1:280343/128‑1 (MQ=255)
                                                    ggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcatca                                                                                                                       >  2:280343/1‑128 (MQ=255)
                                                     gcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGc                                                                                                 >  1:859837/1‑149 (MQ=255)
                                                        aaCTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAg                                                                                              <  2:1003108/149‑1 (MQ=255)
                                                         aCTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGt                                                                                             <  2:657963/149‑1 (MQ=255)
                                                          cTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGtt                                                                                            <  2:831048/149‑1 (MQ=255)
                                                          cTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGtt                                                                                            <  4:956690/149‑1 (MQ=255)
                                                           tGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTa                                                                                                                                  <  1:103622/110‑1 (MQ=255)
                                                              acaGATGCCCAGCGTCGCGAAGAGCTTCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAc                                                                                                                           >  2:435544/1‑114 (MQ=255)
                                                              acaGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCt                                                                                        <  3:1645000/149‑1 (MQ=255)
                                                                      ccAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGATCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                <  1:479441/149‑1 (MQ=255)
                                                                              gcgAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTc                                                                        >  4:1617539/1‑149 (MQ=255)
                                                                                 aaGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGAc                                                                     <  2:1616350/149‑1 (MQ=255)
                                                                                   gagCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCa                                                                   <  1:1052482/149‑1 (MQ=255)
                                                                                    agCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAg                                                                  <  1:1241500/149‑1 (MQ=255)
                                                                                     gctgctGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGt                                                                 <  1:1018884/149‑1 (MQ=255)
                                                                                      ctgctgACTATTGCAGAGATCTCACGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGta                                                                >  2:1047186/1‑149 (MQ=255)
                                                                                       tgctgACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGtat                                                               >  1:171638/1‑149 (MQ=255)
                                                                                           taCTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGTAATACGAATCCAACGCCAGTTAGCTATCGTTAGGGCGCTTCGACCAGTATATGt                                                           >  3:328443/2‑149 (MQ=255)
                                                                                           gACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGt                                                           <  4:816908/149‑1 (MQ=255)
                                                                                            aCTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                >  2:332471/1‑127 (MQ=255)
                                                                                            aCTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                <  1:332471/127‑1 (MQ=255)
                                                                                            aCTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGtt                                                          >  1:472177/1‑149 (MQ=255)
                                                                                              tATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGc                                                        >  2:583474/1‑149 (MQ=255)
                                                                                                  gCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGtt                                                    >  4:554588/1‑149 (MQ=255)
                                                                                                  gCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGtt                                                    <  1:214959/149‑1 (MQ=255)
                                                                                                    agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATa                                                                                                             <  2:172991/90‑1 (MQ=255)
                                                                                                    agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATa                                                                                                             >  1:172991/1‑90 (MQ=255)
                                                                                                    agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                  <  1:207861/149‑1 (MQ=255)
                                                                                                     gagaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGaa                                                                                                                            >  3:1320030/1‑74 (MQ=255)
                                                                                                     gagaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGaa                                                                                                                            <  4:1320030/74‑1 (MQ=255)
                                                                                                     gagaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTa                                                 >  1:548155/1‑149 (MQ=255)
                                                                                                      agaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAAt                                                                                                        <  3:70560/93‑1 (MQ=255)
                                                                                                      agaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTAt                                                >  3:361164/1‑149 (MQ=255)
                                                                                                       gaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATc                                               >  4:898441/1‑149 (MQ=255)
                                                                                                        aTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCa                                              <  3:298787/149‑1 (MQ=255)
                                                                                                          ctcCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGa                                            <  3:958966/149‑1 (MQ=255)
                                                                                                            cccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAca                                          <  2:472177/149‑1 (MQ=255)
                                                                                                              cGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATTCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcatc                                        >  4:274587/1‑149 (MQ=255)
                                                                                                                   tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                   <  4:371898/149‑1 (MQ=255)
                                                                                                                   tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                   >  3:1068812/1‑149 (MQ=255)
                                                                                                                   tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                   <  3:785656/149‑1 (MQ=255)
                                                                                                                   tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCCGCAATACGAATCCAACGCCAGTTCGCTATCGTTg                                                                                 >  3:277754/1‑103 (MQ=255)
                                                                                                                       gcgcAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCa                               >  4:416605/1‑149 (MQ=255)
                                                                                                                        cgcAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAg                              >  1:25044/1‑149 (MQ=255)
                                                                                                                         gcAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAgg                             >  4:668375/1‑149 (MQ=255)
                                                                                                                          cAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                   >  1:392374/1‑142 (MQ=255)
                                                                                                                          cAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAggg                            <  2:984822/149‑1 (MQ=255)
                                                                                                                           aaCATAAGCCGCAGACGTTCTGGTAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGc                           >  2:1378779/1‑149 (MQ=255)
                                                                                                                            aCATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCg                          >  4:455407/1‑149 (MQ=255)
                                                                                                                              aTAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAGCGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGaa                        <  4:78792/149‑1 (MQ=255)
                                                                                                                              aTAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGaa                        <  1:1314529/149‑1 (MQ=255)
                                                                                                                              aTAAGCCGCAGACGTTATGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGaa                        <  4:1193502/149‑1 (MQ=255)
                                                                                                                                aaGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATAc                                                                                                            <  3:385184/63‑1 (MQ=255)
                                                                                                                                aaGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGa                      >  1:79420/1‑149 (MQ=255)
                                                                                                                                aaGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGa                      >  1:996731/1‑149 (MQ=255)
                                                                                                                                 aGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGAt                     >  1:1599951/1‑149 (MQ=255)
                                                                                                                                   ccGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACGCAGGGCGAAGATgc                   <  3:130878/149‑1 (MQ=255)
                                                                                                                                   ccGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgc                   <  3:1367799/149‑1 (MQ=255)
                                                                                                                                   ccGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgc                   <  2:1279384/149‑1 (MQ=255)
                                                                                                                                     gCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcgg                 <  4:1587576/149‑1 (MQ=255)
                                                                                                                                      cAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcggc                >  2:1243020/1‑149 (MQ=255)
                                                                                                                                        gaataacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGt              >  3:1520576/7‑149 (MQ=255)
                                                                                                                                           gTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTcc           >  3:613002/1‑149 (MQ=255)
                                                                                                                                            ttCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCt          <  3:1368988/149‑1 (MQ=255)
                                                                                                                                            aacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCt          >  3:1692995/3‑149 (MQ=255)
                                                                                                                                             tCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTg         <  3:1137027/149‑1 (MQ=255)
                                                                                                                                                ggCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa      <  2:87354/149‑1 (MQ=255)
                                                                                                                                                ggCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGACGCGGCGTTCCTGaaa      <  3:1548881/149‑1 (MQ=255)
                                                                                                                                                ggCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa      >  1:417579/1‑149 (MQ=255)
                                                                                                                                                   aGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGaa   >  3:486849/1‑149 (MQ=255)
                                                                                                                                                   aGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGaa   <  4:967158/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGCCTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAAC  >  NC_000913/4144584‑4144880

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: