Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,299,062 G→C 5.3% A102G (GCC→GGC)  fdhF ← formate dehydrogenase‑H, selenopolypeptide subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,299,0620GC 5.3% 392.4 / 8.5 112A102G (GCC→GGC) fdhFformate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base G (54/52);  new base C (3/3);  total (57/55)
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.99e-01

Only 100 of 113 total aligned reads displayed.

ATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTT  >  NC_000913/4298916‑4299210
                                                                                                                                                  |                                                                                                                                                    
aTTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCa                                                                                                                                                    >  3:247278/1‑149 (MQ=255)
 ttGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCAt                                                                                                                                                   >  2:196750/1‑149 (MQ=255)
       cAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGAc                                                                                                                                             <  2:1230657/149‑1 (MQ=255)
       cAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGAc                                                                                                                                             <  4:247278/149‑1 (MQ=255)
           ccTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTa                                                                                                                                         <  1:497738/149‑1 (MQ=255)
            cTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTAc                                                                                                                                        <  4:909930/149‑1 (MQ=255)
            cTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCCGGATGCCATCCAGAc                                                                                                                                             <  1:624214/144‑1 (MQ=255)
            cTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCCGGATGCCATCCAGAc                                                                                                                                             <  1:287106/144‑1 (MQ=255)
            cTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCCGGATGCCATCCAGAc                                                                                                                                             >  2:624214/1‑144 (MQ=255)
            cTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCCGGATGCCATCCAGAc                                                                                                                                             >  2:287106/1‑144 (MQ=255)
             tGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCt                                                                                                                                       >  3:1448231/1‑149 (MQ=255)
             tGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACt                                                                                                                                       >  4:291149/1‑149 (MQ=255)
             tGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACt                                                                                                                                       >  3:879755/1‑149 (MQ=255)
                 aCCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTtctc                                                                                                                                   >  1:511324/1‑149 (MQ=255)
                   cGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCtt                                                                                                                                 <  4:545564/149‑1 (MQ=255)
                    gATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCTAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCttt                                                                                                                                <  4:1219910/149‑1 (MQ=255)
                    gATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCttt                                                                                                                                <  4:1448231/149‑1 (MQ=255)
                         gCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATgg                                                                                                                           <  2:511324/149‑1 (MQ=255)
                          ccGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCAt                                                                                                                                                   >  2:513048/1‑124 (MQ=255)
                          ccGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGgc                                                                                                                          >  3:1285330/1‑149 (MQ=255)
                           cGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGgcg                                                                                                                         <  4:1285330/149‑1 (MQ=255)
                                 aGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTctct                                                                                                                                  >  1:1319686/1‑134 (MQ=255)
                                 aGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTctct                                                                                                                                  <  2:1319686/134‑1 (MQ=255)
                                 aGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAgg                                                                                                                   <  3:291149/149‑1 (MQ=255)
                                  gACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGgc                                                                                                                  <  3:929174/149‑1 (MQ=255)
                                       gAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCg                                                                                                             >  3:1415985/1‑149 (MQ=255)
                                             agcagTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGAt                                                                                                                             >  2:16888/1‑127 (MQ=255)
                                                  tCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAAtt                                                                                                  >  3:928659/1‑149 (MQ=255)
                                                    aaCGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCa                                                                                                <  2:944947/149‑1 (MQ=255)
                                                     aCGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAg                                                                                               >  3:113596/1‑149 (MQ=255)
                                                       gTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATg                                                                                                                            <  2:1365222/118‑1 (MQ=255)
                                                       gTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATg                                                                                                                            >  1:1365222/1‑118 (MQ=255)
                                                           ttGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCCGGATGCCATcca                                                                                                                                                >  2:61897/1‑93 (MQ=255)
                                                           ttGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCCGGATGCCATcca                                                                                                                                                <  1:61897/94‑2 (MQ=255)
                                                              gTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTctct                                                                                                                                  >  1:185914/1‑105 (MQ=255)
                                                              gTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTctct                                                                                                                                  <  2:185914/105‑1 (MQ=255)
                                                              gTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATc                                                                                      <  4:113596/149‑1 (MQ=255)
                                                               tACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATcc                                                                                     <  4:1427609/149‑1 (MQ=255)
                                                                aCCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATccc                                                                                    >  4:619047/1‑149 (MQ=255)
                                                                  cAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAg                                                                                  <  4:1415985/149‑1 (MQ=255)
                                                                   aaTAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGa                                                                                                                                              <  2:906390/88‑1 (MQ=255)
                                                                   aaTAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGa                                                                                                                                              >  1:906390/1‑88 (MQ=255)
                                                                   aaTAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAgg                                                                                 >  1:659668/1‑149 (MQ=255)
                                                                         ggCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAg                                                                           <  2:659668/149‑1 (MQ=255)
                                                                                gcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGa                                                                                                                              >  3:124142/1‑91 (MQ=255)
                                                                                gcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGa                                                                                                                              <  4:124142/91‑1 (MQ=255)
                                                                                      ttttGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGa                                                                                <  1:1148602/131‑1 (MQ=255)
                                                                                      ttttGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGa                                                                                >  2:1148602/1‑131 (MQ=255)
                                                                                           cATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAg                                                                                                                    <  2:964305/90‑1 (MQ=255)
                                                                                           cATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAg                                                                                                                    <  3:831863/90‑1 (MQ=255)
                                                                                           cATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAg                                                                                                                    >  4:831863/1‑90 (MQ=255)
                                                                                           cATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAg                                                                                                                    >  1:964305/1‑90 (MQ=255)
                                                                                             ttACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGtt                                                                        >  3:1144857/1‑132 (MQ=255)
                                                                                             ttACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGtt                                                                        <  4:1144857/132‑1 (MQ=255)
                                                                                             ttACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGc                                                       >  2:519935/1‑149 (MQ=255)
                                                                                               aCATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTg                                                     >  3:1195302/1‑149 (MQ=255)
                                                                                                cATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTgg                                                    <  4:1195302/149‑1 (MQ=255)
                                                                                                 aTAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTtctc                                                                                                                                   >  2:88254/1‑69 (MQ=255)
                                                                                                    gTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCt                                                                                                               <  2:622034/86‑1 (MQ=255)
                                                                                                    gTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCt                                                                                                               >  1:622034/1‑86 (MQ=255)
                                                                                                     ttGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGa                                                                     <  1:1372024/127‑1 (MQ=255)
                                                                                                     ttGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGa                                                                     >  2:1372024/1‑127 (MQ=255)
                                                                                                         tttCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATc                                           >  2:687577/1‑149 (MQ=255)
                                                                                                           tCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCAt                                         >  2:770868/1‑149 (MQ=255)
                                                                                                             gTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGt                                                                                              >  1:1369630/1‑94 (MQ=255)
                                                                                                               tACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATgggg                                     <  1:770868/149‑1 (MQ=255)
                                                                                                               tACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATgggg                                     <  1:196750/149‑1 (MQ=255)
                                                                                                               tACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATgggg                                     <  3:619047/149‑1 (MQ=255)
                                                                                                                  ccGTACCACGCTAGGAGCCGGTCGTCTGGATGGCATCCGCACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGCaa                                                                                                    <  2:166839/83‑1 (MQ=255)
                                                                                                                  ccGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTaa                                                                                                    >  1:166839/1‑83 (MQ=255)
                                                                                                                  ccGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAAtt                                                                                                  <  1:1145216/85‑1 (MQ=255)
                                                                                                                  ccGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAAtt                                                                                                  >  2:1145216/1‑85 (MQ=255)
                                                                                                                  ccGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAAtt                                                                                                  <  4:1137928/85‑1 (MQ=255)
                                                                                                                  ccGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAAtt                                                                                                  >  3:1137928/1‑85 (MQ=255)
                                                                                                                      aCCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCa                              >  1:479145/1‑149 (MQ=255)
                                                                                                                       ccACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATcc                                                                                     <  4:1288179/93‑1 (MQ=255)
                                                                                                                       ccACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATcc                                                                                     >  3:1288179/1‑93 (MQ=255)
                                                                                                                       ccACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATgggg                                     >  4:198142/1‑141 (MQ=255)
                                                                                                                       ccACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATgggg                                     <  3:198142/141‑1 (MQ=255)
                                                                                                                       ccACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAg                             <  1:687577/149‑1 (MQ=255)
                                                                                                                         aCGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGgc                           >  3:336276/1‑149 (MQ=255)
                                                                                                                         aCGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGgc                           >  1:374619/1‑149 (MQ=255)
                                                                                                                                 agCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCa                              <  2:1400568/138‑1 (MQ=255)
                                                                                                                                 agCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCa                              >  1:1400568/1‑138 (MQ=255)
                                                                                                                                 agCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTc                   <  1:519935/149‑1 (MQ=255)
                                                                                                                                   ccGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTc                                                                       >  3:212917/1‑95 (MQ=255)
                                                                                                                                   ccGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTc                                                                       <  4:212917/95‑1 (MQ=255)
                                                                                                                                     ggTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGa               >  4:741432/1‑149 (MQ=255)
                                                                                                                                      gtcgtcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGAt              >  4:1014484/1‑149 (MQ=255)
                                                                                                                                       tcgtcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGcgcg                        >  2:193020/1‑138 (MQ=255)
                                                                                                                                       tcgtcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGcgcg                        <  1:193020/138‑1 (MQ=255)
                                                                                                                                       tcgtcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATc             >  3:178513/1‑149 (MQ=255)
                                                                                                                                          tcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTgg          <  3:741432/149‑1 (MQ=255)
                                                                                                                                          tcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTgg          <  4:879755/149‑1 (MQ=255)
                                                                                                                                            tGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGt        >  4:623860/1‑149 (MQ=255)
                                                                                                                                            tGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGt        <  3:1014484/149‑1 (MQ=255)
                                                                                                                                             ggATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTa       >  1:865842/1‑149 (MQ=255)
                                                                                                                                              gATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAg                                                                           <  4:23648/80‑1 (MQ=255)
                                                                                                                                              gATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAg                                                                           >  3:23648/1‑80 (MQ=255)
                                                                                                                                                  gCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGtt  >  2:784790/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTT  >  NC_000913/4298916‑4299210

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