Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,085,619 G→T 7.7% A68E (GCA→GAA)  fdoG ← formate dehydrogenase‑O, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,085,6190GT 7.7% 229.8 / 2.4 86A68E (GCA→GAA) fdoGformate dehydrogenase‑O, large subunit
Reads supporting (aligned to +/- strand):  ref base G (34/44);  new base T (6/2);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGG  >  NC_000913/4085471‑4085766
                                                                                                                                                    |                                                                                                                                                    
aCCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTg                                                                                                                                                      <  1:935004/149‑1 (MQ=255)
aCCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTg                                                                                                                                                      <  1:25497/149‑1 (MQ=255)
   tGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTc                                                                                                                                                   <  2:990762/149‑1 (MQ=255)
    ggcggccGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGTCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTcc                                                                                                                                                  <  2:1054664/143‑1 (MQ=255)
     gcgcACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGAGAGATGCTTTGGCGTTTATCGCTCCg                                                                                                                                                 >  1:27988/1‑149 (MQ=255)
      cgcACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGt                                                                                                                                                <  2:135566/149‑1 (MQ=255)
        cACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCa                                                                                                                                              <  2:1379862/149‑1 (MQ=255)
        cACGGTAGTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGTAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCa                                                                                                                                              <  2:269767/149‑1 (MQ=255)
           ggTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCg                                                                                                                                           >  2:452897/1‑149 (MQ=255)
              agaccGGAAGCTGCGGAGGGCTTTCGGAGGGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCGTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAgg                                                                                                                                        <  2:260787/146‑1 (MQ=255)
                  cGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGt                                                                                                                                    <  1:1261167/149‑1 (MQ=255)
                   ggAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTa                                                                                                                                   <  1:496565/149‑1 (MQ=255)
                    gAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTAc                                                                                                                                  >  1:1406154/1‑149 (MQ=255)
                       aCTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGCTCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATc                                                                                                                               <  1:973011/149‑1 (MQ=255)
                        cTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTGGGCGTTTTTGGCTCCGTCACCGAGGCGGtcggcca                                                                                                                              >  1:441320/1‑143 (MQ=255)
                         ttCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCCTCGATAGTGAAGATAGATGCCTTGGCGTTTTTTGCTCCGTCACCGAGGCGGTACATCaa                                                                                                                             >  2:447575/1‑149 (MQ=255)
                           cAGACGGCTTTCGGGGTGGATGAAATCCACCAGGCCAGCGCCTTGCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTCTGCTCCGTCACCGAGGCTGTACATCAACa                                                                                                                           <  1:1026072/149‑1 (MQ=255)
                             gACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGc                                                                                                                         <  1:542961/149‑1 (MQ=255)
                             gACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGc                                                                                                                         >  2:1279731/1‑149 (MQ=255)
                                ggCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCGGTACATCAACAGCCca                                                                                                                      >  1:1006826/1‑149 (MQ=255)
                                   tttCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAg                                                                                                                   >  2:842051/1‑149 (MQ=255)
                                   tttCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAg                                                                                                                   <  1:859925/149‑1 (MQ=255)
                                   tttCGGAGTGGATGAAATCCACCAGGCCAACGCCTTTCGGACAAAGTGCACCGAGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACAGAGACTGTACATCAACAgacgacgg                                                                                                                   >  1:889280/1‑142 (MQ=255)
                                    ttCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGc                                                                                                                  <  1:426145/149‑1 (MQ=255)
                                    ttCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGc                                                                                                                  <  1:67753/149‑1 (MQ=255)
                                            ggATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCCGATCGCCTTCGATATGGAAGATAGATGCTTTTGCGTTTTTTGCTCCGTCACCGAGGCCGTACATCACGAGCCCCAAGCCCACgcgc                                                                                                          >  1:1132159/1‑146 (MQ=255)
                                             gATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAAc                                                                                                         >  2:99181/1‑149 (MQ=255)
                                               tGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACaa                                                                                                       >  2:259960/1‑149 (MQ=255)
                                                 aaaTCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATa                                                                                                     >  2:240313/1‑149 (MQ=255)
                                                 aaaCCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGCGGCAGTACATCAACAGCCCACAGCCTACGGAACAATa                                                                                                     <  2:676835/149‑1 (MQ=255)
                                                        ccaGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAg                                                                                              <  2:627070/149‑1 (MQ=255)
                                                         caGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTCACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGCGGCTGCCCACCAACAGCCCACAGCCTCCGGAACAATAGGTGCAgg                                                                                             >  2:305125/1‑149 (MQ=255)
                                                            gCCAGCGCCTTTCGGACAAAGTGCACCGCGGTGGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTTTTTTTTTTCCTAGGCTGTACATCAACAGCCCACAGCCTACGGAGCAATAGGGGCAGGtct                                                                                          >  2:1043209/1‑147 (MQ=255)
                                                              cAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAGGATAGATGCTTTGGCGTTTTTTGACCTTGGATCGAGGCTGTGCAACCGCAGCCCAGAGCCTCCGGAACGATAGGGTCAggcgtta                                                                                        >  1:699014/1‑144 (MQ=255)
                                                              cAGCGCCCTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTACGTCACCGAGGCTGTACATCCACAGCCCACAGCCTACGGAACAATAGGTGCAgggagta                                                                                        >  2:1220204/1‑144 (MQ=255)
                                                               aGCGCCTTTCGGACAAAGTGCGCCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATAGGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTAc                                                                                       >  2:353088/1‑149 (MQ=255)
                                                                gcgcCCTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTTCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACg                                                                                      <  1:1431092/149‑1 (MQ=255)
                                                                    cTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGGTTTTTGCTCCGTCACCGAGGCTGTGCATCGACAGCCCACAGCCTACGGGACAAGAGGTGGGGGTATTCCGGGtt                                                                                  >  1:1103958/1‑149 (MQ=255)
                                                                        cGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTcgc                                                                              <  1:356754/149‑1 (MQ=255)
                                                                         ggACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTcgcg                                                                             <  1:1233935/149‑1 (MQ=255)
                                                                           aCAAAGTGCACCGCGGTTGACCGGGTGATACGGATCGCCTTCGATATGGAAGACAGATGCTTTGGCGTTTTTTGCTGCGTCACCGAGGCTGTACATACACAGCCCACAGCCTACGGTACAATAGGGGCAGGTATTACGGGTTTCGCggg                                                                           >  2:803544/1‑149 (MQ=255)
                                                                               aGTGCACCGGGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCCACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTgcg                                                                       <  2:83889/149‑1 (MQ=255)
                                                                                    agcgcgGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCgccgca                                                                  <  2:179992/147‑1 (MQ=255)
                                                                                      cgggtttGACCGGTTGAGCCGGACCGCCTTCGATATGGAAGCTAGATGCTTTGGCGTTTTTTTCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGt                                                                <  2:997821/144‑1 (MQ=255)
                                                                                      cgcgGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTTTTTCTTTTTCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCGGt                                                                >  2:1442210/1‑149 (MQ=255)
                                                                                      cgcgGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTTCGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGt                                                                <  1:938543/149‑1 (MQ=255)
                                                                                        cgGGTGGCCGGGTGACCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGttt                                                              <  2:1235562/149‑1 (MQ=255)
                                                                                         gggttaCCGGGTGAGTCGGATGGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCCCCGAGGCTGTACACCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTta                                                             <  2:554335/144‑1 (MQ=255)
                                                                                         ggTTGACGGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTta                                                             <  1:452897/149‑1 (MQ=255)
                                                                                           ttGACCGGGTGATCCAGATCGCGTTCGATATAGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTtata                                                           <  1:1386844/149‑1 (MQ=255)
                                                                                               ccGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCACCAGTTTACACTGc                                                       >  2:122633/1‑149 (MQ=255)
                                                                                                  ggTGCTCCGGATCGGCTTCGATATGGAAGATAGATGCTTTGGCGTTTTCTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACCGCCTACGGAACAAGAGGTGCAGGGATTACGGGTTTCGCGGGTGGGAAGCAGTTTACACTGACgg                                                    <  2:256585/149‑1 (MQ=255)
                                                                                                  ggTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCgg                                                    <  2:1269964/149‑1 (MQ=255)
                                                                                                  ggTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCgg                                                    <  2:154205/149‑1 (MQ=255)
                                                                                                  ggTGATCCGGATCGCCTTCGAGATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCgg                                                    <  1:482717/149‑1 (MQ=255)
                                                                                                  ggTGATCCGGATCGCCGTCGATATGGAAGATAGACGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGTCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCgg                                                    <  1:349860/149‑1 (MQ=255)
                                                                                                      aTCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTTCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGCGCAAGTATCACGGGTTTCGCGGGCGCGCAGCAGTTTATACTCCCgggcta                                                >  2:571208/1‑146 (MQ=255)
                                                                                                      aTCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCGGTACATCAACAGCCCACAGCCTACGGAACACTAGGTGCAGGTATTGCGGGGTTCGCGGGTGCGCCGCAGTTTATACTGCCGGGgtt                                                >  2:560450/1‑149 (MQ=255)
                                                                                                           gATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGCTTCcgcg                                           >  2:213767/1‑149 (MQ=255)
                                                                                                            aTCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGa                                          <  2:1406154/149‑1 (MQ=255)
                                                                                                                ccTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTTTTTTTTTTTCGAGGCTGTACATCAACAGCCCACAGCATACGGGACAATGGGTGCAGGGATTGCGGGTTTCGCGGGTGCGCCGCCGTTTGTACTGCCGGGGTTCCGCGCGTgt                                     >  2:123985/1‑148 (MQ=255)
                                                                                                                ccTTCGATATGGAAGATAGATGCTTTAGCGCCTCACGAGGCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGc                                     >  1:1440973/1‑149 (MQ=255)
                                                                                                                 cTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCCCAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGGGCGAGGCAGTTTATACTGCCGGGTTTCCGCGGGTGCt                                    >  1:19850/1‑149 (MQ=255)
                                                                                                                 cTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTGCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCt                                    >  1:492290/1‑149 (MQ=255)
                                                                                                                     gATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGGTTATACTGCCCGGATTCCGCGCGTGCTACGc                                >  2:1311091/1‑149 (MQ=255)
                                                                                                                        atGGAAGATGGATGCTTGGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTgg                             <  2:441320/149‑1 (MQ=255)
                                                                                                                         tGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTggg                            <  2:192706/149‑1 (MQ=255)
                                                                                                                         tGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGGTTCCGCGAGTGCTACGCTggg                            >  1:718785/1‑149 (MQ=255)
                                                                                                                           gAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGCGTGCTACGCTGGGTg                          >  2:157169/1‑149 (MQ=255)
                                                                                                                           gAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTg                          >  1:1151191/1‑149 (MQ=255)
                                                                                                                           gAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTg                          >  1:810473/1‑149 (MQ=255)
                                                                                                                               atagatGCTTTGGCGTTTTTTGCTTCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTGTATACTGCCGGGTTTCCCCGAGTGCTCCGCTGGGTGCaaa                      >  2:749974/1‑149 (MQ=255)
                                                                                                                                tagatGCTTTGGCGTTTTTTTTTTTTTTTCCGAGGCTGTACCTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTAGGGGTTTCGCGGGTTCGCAGCAGTTTATACTACCGTGTTTCCGCGTGTGGTACGCTGGGGGCaaaa                     >  2:58142/1‑149 (MQ=255)
                                                                                                                                 agatGCTTTGGCGTTTTTTTTTTTTTTTTTGAGGCGGTACATCAACAGCCCACAGCCTTCGGTGGAATGGGTGCAGGTATTTCGGGTTTCGCGGGTGGGAAGCTTTGTCTTCTTTCTGGTTTTCGTGAGTGCTAGGCTCGGTGCGAAAc                    >  2:516458/1‑149 (MQ=255)
                                                                                                                                    tGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCa                 <  1:259960/149‑1 (MQ=255)
                                                                                                                                     gCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAg                <  2:1279274/149‑1 (MQ=255)
                                                                                                                                      cTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGAAGGTTATACTACCGGGTTTCCGGGGGGGCTACGCTGAGTGCAATACCGAGt               >  2:776813/1‑149 (MQ=255)
                                                                                                                                       tttGGCGTTTTTTGCTCCGTCACAGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTg              >  1:809799/1‑149 (MQ=255)
                                                                                                                                          ggCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCg           <  2:1426360/149‑1 (MQ=255)
                                                                                                                                           gCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGc          <  1:240313/149‑1 (MQ=255)
                                                                                                                                            cGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCt         <  1:1243382/149‑1 (MQ=255)
                                                                                                                                              ttttttGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTCCCCGCGTAAGTACGCTGGGTGCAAAACCCAGTGCCGCTGc       <  1:353088/149‑1 (MQ=255)
                                                                                                                                              ttttttGCCCCTTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGc       <  1:550283/149‑1 (MQ=255)
                                                                                                                                               ttttgGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGTGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAACCCTCGTGTCGCTGcc      <  2:27988/149‑1 (MQ=255)
                                                                                                                                                ttttGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCg     >  2:267610/1‑149 (MQ=255)
                                                                                                                                                 tttGCCCCGTCACCGAGGCTGTACATCAACAGCCCAAAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCgt    <  2:699014/149‑1 (MQ=255)
                                                                                                                                                  ttGCGCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCgtg   <  1:1454344/149‑1 (MQ=255)
                                                                                                                                                  tcgCTCCGTCCCCGAGGCTGTACAGCAACAGCCCACAGCCGACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCgtg   <  2:888884/147‑1 (MQ=255)
                                                                                                                                                   tGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCgtgg  <  1:1437637/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGG  >  NC_000913/4085471‑4085766

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: