Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,144,735 G→C 5.6% 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 5.6% 555.6 / 15.4 157A247A (GCG→GCCpflDputative glycine radical domain‑containing pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (66/82);  new base C (4/5);  total (70/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 158 total aligned reads displayed.

CTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGC  >  NC_000913/4144588‑4144883
                                                                                                                                                   |                                                                                                                                                    
cTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGt                                                                                                                                                     >  1:831048/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)
       aaaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaa                                                                                                                                              <  1:1711729/149‑3 (MQ=255)
        aaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaac                                                                                                                                             <  1:764403/149‑4 (MQ=255)
         acacATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGTCAGTTAt                                                                                                                                            >  1:1179554/1‑149 (MQ=255)
          cacaTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGCCGCGAAGGGCCGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTACTTGCACGCCTGCCAgaacgt                                                                                                                                           <  4:138914/149‑6 (MQ=255)
             aTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTg                                                                                                                                        <  1:1688465/149‑1 (MQ=255)
                   cGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACAt                                                                                                                                  <  4:326682/149‑1 (MQ=255)
                     ttACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGa                                                                                                                                >  4:1520658/1‑149 (MQ=255)
                       aCGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAc                                                                                                                              <  3:1485516/149‑1 (MQ=255)
                         gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcat                                                                                                                            >  2:374627/1‑149 (MQ=255)
                         gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcat                                                                                                                            >  3:623762/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)
                                 ggCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGc                                                                                                                    <  2:191027/149‑1 (MQ=255)
                                   cGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCaa                                                                                                                  <  1:620934/149‑1 (MQ=255)
                                         cGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGaa                                                                                                            <  2:1473140/149‑1 (MQ=255)
                                           aTGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATc                                                                                                          >  4:178769/1‑149 (MQ=255)
                                           aTGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATc                                                                                                          >  1:796265/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)
                                                    aaCTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAg                                                                                                 >  4:351835/1‑149 (MQ=255)
                                                     aCTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGt                                                                                                <  2:657963/149‑1 (MQ=255)
                                                      cTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGtt                                                                                               <  2:831048/149‑1 (MQ=255)
                                                      cTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGtt                                                                                               <  4:956690/149‑1 (MQ=255)
                                                      cTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGtt                                                                                               >  3:420731/1‑149 (MQ=255)
                                                       tGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTa                                                                                                                                     <  1:103622/110‑1 (MQ=255)
                                                          acaGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAc                                                                                                                              <  1:435544/114‑1 (MQ=255)
                                                          acaGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCt                                                                                           <  3:1645000/149‑1 (MQ=255)
                                                              aTGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCg                                                                                       <  3:61857/149‑1 (MQ=255)
                                                                  ccAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGATCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                   <  1:479441/149‑1 (MQ=255)
                                                                          gcgAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTc                                                                           >  4:1617539/1‑149 (MQ=255)
                                                                            gAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGa                                                                         <  2:1054756/149‑1 (MQ=255)
                                                                             aaGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGAc                                                                        <  2:1616350/149‑1 (MQ=255)
                                                                               gagCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCa                                                                      >  3:232492/1‑149 (MQ=255)
                                                                               gagCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCa                                                                      <  1:1052482/149‑1 (MQ=255)
                                                                                agCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAg                                                                     <  1:1241500/149‑1 (MQ=255)
                                                                                agCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAg                                                                     <  2:311121/149‑1 (MQ=255)
                                                                                agCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAg                                                                     >  1:1304891/1‑149 (MQ=255)
                                                                                 gctgctGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGt                                                                    <  1:1018884/149‑1 (MQ=255)
                                                                                  ctgctgACTATTGCAGAGATCTCACGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGta                                                                   >  2:1047186/1‑149 (MQ=255)
                                                                                       taCTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGTAATACGAATCCAACGCCAGTTAGCTATCGTTAGGGCGCTTCGACCAGTATATGt                                                              >  3:328443/2‑149 (MQ=255)
                                                                                       gACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGt                                                              <  4:816908/149‑1 (MQ=255)
                                                                                        aCTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                   <  1:332471/127‑1 (MQ=255)
                                                                                        aCTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGCTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGtt                                                             <  3:1417760/149‑1 (MQ=255)
                                                                                          tATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGc                                                           >  2:583474/1‑149 (MQ=255)
                                                                                              gCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGtt                                                       <  1:214959/149‑1 (MQ=255)
                                                                                              gCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGtt                                                       >  4:554588/1‑149 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATa                                                                                                                >  1:172991/1‑90 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                     >  2:456522/1‑149 (MQ=255)
                                                                                                agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGCTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                     >  4:1553491/1‑149 (MQ=255)
                                                                                                 gagaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGaa                                                                                                                               <  4:1320030/74‑1 (MQ=255)
                                                                                                 gagaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGaa                                                                                                                               >  3:1320030/1‑74 (MQ=255)
                                                                                                 gagaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTa                                                    >  1:548155/1‑149 (MQ=255)
                                                                                                  agaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAAt                                                                                                           <  3:70560/93‑1 (MQ=255)
                                                                                                  agaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAAt                                                                                                           >  4:70560/1‑93 (MQ=255)
                                                                                                  agaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTAt                                                   >  3:361164/1‑149 (MQ=255)
                                                                                                   gaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATc                                                  >  4:898441/1‑149 (MQ=255)
                                                                                                   gaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATc                                                  <  1:374627/149‑1 (MQ=255)
                                                                                                    aTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCa                                                 <  3:298787/149‑1 (MQ=255)
                                                                                                      ctcCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGa                                               <  3:958966/149‑1 (MQ=255)
                                                                                                       tcCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAc                                              >  2:983311/1‑149 (MQ=255)
                                                                                                          cGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATTCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcatc                                           >  4:274587/1‑149 (MQ=255)
                                                                                                               tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                      <  3:785656/149‑1 (MQ=255)
                                                                                                               tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                      <  4:371898/149‑1 (MQ=255)
                                                                                                               tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                      >  3:1068812/1‑149 (MQ=255)
                                                                                                               tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCCGCAATACGAATCCAACGCCAGTTCGCTATCGTTg                                                                                    >  3:277754/1‑103 (MQ=255)
                                                                                                                   gcgcAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCa                                  >  4:416605/1‑149 (MQ=255)
                                                                                                                     gcAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAgg                                >  4:668375/1‑149 (MQ=255)
                                                                                                                      cAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTaa                                      <  2:392374/142‑1 (MQ=255)
                                                                                                                      cAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAggg                               <  2:984822/149‑1 (MQ=255)
                                                                                                                       aaCATAAGCCGCAGACGTTCTGGTAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGc                              >  2:1378779/1‑149 (MQ=255)
                                                                                                                        aCATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCg                             >  4:455407/1‑149 (MQ=255)
                                                                                                                          aTAAGCCGCAGACGTTATGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGaa                           <  4:1193502/149‑1 (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)
                                                                                                                                cGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcg                     <  1:1743210/149‑1 (MQ=255)
                                                                                                                                    gaataacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGt                 >  3:1520576/7‑149 (MQ=255)
                                                                                                                                        aacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCt             >  3:1692995/3‑149 (MQ=255)
                                                                                                                                         tCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTg            <  3:1137027/149‑1 (MQ=255)
                                                                                                                                            ggCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGACGCGGCGTTCCTGaaa         <  3:1548881/149‑1 (MQ=255)
                                                                                                                                            ggCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACACCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa         <  1:222836/149‑1 (MQ=255)
                                                                                                                                            ggCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa         >  1:417579/1‑149 (MQ=255)
                                                                                                                                             gCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAg        <  2:1184608/149‑1 (MQ=255)
                                                                                                                                               aGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGaa      <  4:967158/149‑1 (MQ=255)
                                                                                                                                               aGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGaa      >  3:486849/1‑149 (MQ=255)
                                                                                                                                                gtcGTGCCAGTTATTCTCGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGATTCGACCAGTATATGTTGCAGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCAGCGTTCTTGAAAGAAc     <  2:269787/147‑1 (MQ=255)
                                                                                                                                                  cGTGCCAGTTATTCTCGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTg   <  2:547792/149‑1 (MQ=255)
                                                                                                                                                   gTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGc  >  2:1661407/1‑149 (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 < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: