Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007779_cfb 4,371,274 (A)5→3 intergenic (‑9/+296) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779_cfb4,371,2710A.83.3% 86.8 / 22.2 36intergenic (‑6/+300)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/6);  new base . (27/3);  total (27/9)
Fisher's exact test for biased strand distribution p-value = 4.31e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_007779_cfb4,371,2720A.83.3% 90.0 / 19.0 36intergenic (‑7/+299)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/6);  new base . (27/3);  total (27/9)
Fisher's exact test for biased strand distribution p-value = 4.31e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTA  >  NC_007779_cfb/4371022‑4371507
                                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                           
aTGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTTTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGa                                                                                                                                                                                                                                              >  1:285876/1‑250 (MQ=255)
   cGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCGGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGATCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                                                                                                                                                                           <  2:213423/250‑2 (MQ=255)
   cGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                                                                                                                                                                           >  2:361693/1‑249 (MQ=255)
       gATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaatagcg                                                                                                                                                                                                                                       <  1:361693/250‑6 (MQ=255)
            cAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                                                                                                                                                                >  1:409471/1‑250 (MQ=255)
             aaTCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTATCGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                                                                                                                                                                                               >  1:145134/1‑250 (MQ=255)
                cTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTATCGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAcc                                                                                                                                                                                                                            <  2:145134/250‑1 (MQ=255)
                 tGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCg                                                                                                                                                                                                                           <  2:90044/250‑1 (MQ=255)
                 tGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCg                                                                                                                                                                                                                           >  1:403400/1‑250 (MQ=255)
                             aCTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATgg                                                                                                                                                                                                               >  1:1309702/1‑250 (MQ=255)
                              cTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGga                                                                                                                                                                                                              >  1:520533/1‑250 (MQ=255)
                                    tagagcgCGGAGGGCAGGTTTTAACTCGTTGTTCCTCCCAAACTTCAGCGAAAATTGGCGGGGCCGCAGAAGAGAGGTTGCCATTACAAGAAAAGATAGTCAGAAAATAGGTGCCGATGGGTGCGGGGTTGGTGATGTATTGGGGGTTACGGCGAACCAGCTTTTCTGTCTGAGGAACAAGATAGTCCAGCCCGCCAGAAACCTGCAGGGC‑‑AAAAATAGCGGCGATAACCGCCATGATAATGGAGATGAc                                                                                                                                                                                                        <  2:1309702/244‑1 (MQ=255)
                                         aGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGa                                                                                                                                                                                                   >  1:798202/1‑250 (MQ=255)
                                             gACGGCAAGGATTAACGCCTTGTTCCTTCGCAACTTCAGCGATAATTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg                                                                                                                                                                                               >  1:1027617/1‑250 (MQ=255)
                                                 tcAAGTTTTAACGCTTTTTCCCTTCGCATCTACAGAGAAAACGGGGAGGGCCCCGAGAGATATGTTGCCACTCAGAAAAAAGATAGTGAGAAAAAAGGTCAGGTGGGGTGAGGGGATCGGGATTTATTGCGGTTAGCGGCGGAAAAGTTTTCCTGCTTGATGAGCCTGATCGTCGAGACCGCCAAAACCCTGCAGGGAAAAAATGGCGCCGATAGCCGCCATGATAATGGAGATGACATCGAACGGGATg                                                                                                                                                                                           <  2:79485/249‑1 (MQ=255)
                                                           aCGCTTTGTTCTTTCGAAACTTCGAGGATAACTGGCAGTTGCGAAGGGGAGATGTCGCCAGTCCCACAAAAGATATCCAGAAATAGGGCCACGGTCGGGGGGAGGATGGTATTGTTTTCGGGTTTCCGGCCCAGCGGTTTTTCTTTCTAATAACCAAATTACTCCCGCCCGCCGGAACCCTGCAGGGAAAAAATGGGGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGtt                                                                                                                                                                                 <  2:613179/250‑1 (MQ=255)
                                                                   ttCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                                                                                                                                                                                               <  2:640272/196‑1 (MQ=255)
                                                                   ttCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                                                                                                                                                                                               >  1:640272/1‑196 (MQ=255)
                                                                       nccacNACTACAGCGATAGCTGGACGTGCCGCCAGGAAAATGTTGCAAGTACCAGTAAAGATAGTAGGAAAATGGGGCGGGATGGGGGCGAGGGTGGTAATGTTTTTGGGTTTAGGGCGCGGCGTTTTTTCTGTCTGATGACCCAGATAGTCAAGCCCGCCACCAACCGGCTTGGAAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAg                                                                                                                                                                     <  2:798202/246‑1 (MQ=255)
                                                                                        aaCTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCAcc                                                                                                                                                    >  1:898729/1‑250 (MQ=255)
                                                                                        aaCTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCAcc                                                                                                                                                    >  1:577257/1‑250 (MQ=255)
                                                                                           tGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCAccccc                                                                                                                                                 >  1:1175630/1‑250 (MQ=255)
                                                                                           tGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCAccccc                                                                                                                                                 >  1:3224/1‑250 (MQ=255)
                                                                                            ggCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCa                                                                                                                                                >  1:981252/1‑250 (MQ=255)
                                                                                                    cggggaggagaTGTGCCCAGCACAAGAAAAGATAGAAAGAAATGGGTCACGGGGGGGGGGGGGCTTGGTGTTGTTTTTGGGTTTCCGGCACAGCAGCTTTTTTTTCGGAGAACCAAGAAAGCCCAGCCCGCCAGCACCCTGCATGGC‑‑AAAAATAGCGGGAATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCt                                                                                                                                        <  2:159791/243‑1 (MQ=255)
                                                                                                                    ccAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTa                                                                                                                        >  2:442842/1‑250 (MQ=255)
                                                                                                                        tACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCATACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTcc                                                                                                                    >  1:475268/1‑250 (MQ=255)
                                                                                                                                        cAGGAAATAGGTCACGATCAGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGccc                                                                                                    >  1:1308912/1‑250 (MQ=255)
                                                                                                                                          ggAAATAGGTCGCGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCACTCTGGCGCCCaa                                                                                                  >  1:798728/1‑250 (MQ=255)
                                                                                                                                              aTAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCaagaag                                                                                              >  1:554007/1‑250 (MQ=255)
                                                                                                                                                         aTCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCa                                                                                    >  1:380952/1‑250 (MQ=255)
                                                                                                                                                          tCGGTTCTAGGACCGTGATGTATTTCGGGTTACGGCGCGGCCGCTTTTCTGTCTGATGAACCAGATCGTCCAGACCGCCAGCAACCTGCCTGGCAG‑‑AAATAGAGGCGATAACCGCCATGATAATGGAGTTTACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGAGGCAAGAACACGCAACCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGAACGNCANGaa                                                                                   >  2:1095835/1‑250 (MQ=255)
                                                                                                                                                              tGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACt                                                                               >  1:469321/1‑250 (MQ=255)
                                                                                                                                                                cGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTatg                                                                            >  1:540425/1‑250 (MQ=255)
                                                                                                                                                                                         ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCTAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATAATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCatat                                                      >  1:203655/1‑250 (MQ=255)
                                                                                                                                                                                             ggCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGc                                                  >  1:1306540/1‑250 (MQ=255)
                                                                                                                                                                                               cgcAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGcct                                                >  1:916220/1‑250 (MQ=255)
                                                                                                                                                                                                          ttCTGTCTGATGAACCAGATAGCCCAGCCCGCCAGCACCCTGCATGGC‑‑AAAAATAGGGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGttttt                                     <  1:510120/250‑1 (MQ=255)
                                                                                                                                                                                                              gtctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATAATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAAc                                 >  2:99919/1‑250 (MQ=255)
                                                                                                                                                                                                                                               cAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTa  >  1:404155/1‑250 (MQ=255)
                                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                           
ATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTA  >  NC_007779_cfb/4371022‑4371507

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: