Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,113,851 G→A H208Y (CAT→TAT)  fpr ← flavodoxin/ferredoxin‑NADP(+) reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,113,8510GA94.1% 206.4 / ‑0.0 69H208Y (CAT→TAT) fprflavodoxin/ferredoxin‑NADP(+) reductase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (32/32);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTCCTTGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTAG  >  NC_000913/4113606‑4114094
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
aCCTCCTTGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTgg                                                                                                                                                                                                                                                 >  1:710376/1‑251 (MQ=255)
       tGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCtt                                                                                                                                                                                                                                           <  1:694681/250‑1 (MQ=255)
               gTCATCATCACTCATTTAACAACCGTTAGATGCCGGAATATCGCGGGATGAACAATAAGATGATTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGGCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTGCTTTATTCATCg                                                                                                                                                                                                                                  <  2:479180/251‑1 (MQ=255)
               gTCATCATCACTCATATATCAACCGTTTGATGCCGGCATATCGCGAGATGATCAATAAGATGAATGCGGCTTTGCTCATCGTGCCGTTTAACGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTAGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCCGCAACTTTTGTGTATCGAGCACCATCTATGGATTGCCGCACAGCATCACCTCGCTGGTTTCTTTATTCATCg                                                                                                                                                                                                                                  >  2:396054/1‑251 (MQ=255)
                 catcatcaCTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGc                                                                                                                                                                                                                                <  1:228980/251‑1 (MQ=255)
                    caccaCTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACAGCTTACCAGTCATGCTCCGCTGTCATATGGCCAGGTCCGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCGGCAACTCTTGTGTATCGCGCCCCATCTGTGGATTTCCGCACAGCGTCACATGGCTGGCTTCTTTCTTCATCTGCAgg                                                                                                                                                                                                                             >  2:259730/4‑251 (MQ=255)
                         cTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGGGGATTGCCGCACAGCATCAACTAGCTGGGTTCTTTATTCATCGGCAGGCCaac                                                                                                                                                                                                                        >  2:545880/1‑250 (MQ=255)
                                   aaCCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTcc                                                                                                                                                                                                              <  1:400237/251‑1 (MQ=255)
                                     ccGTTTGATGCCGGCATATCGCGGGATGAACAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAg                                                                                                                                                                                                            <  2:260745/251‑1 (MQ=255)
                                            aTGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTccccc                                                                                                                                                                                                     >  1:322638/1‑250 (MQ=255)
                                              gCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCAc                                                                                                                                                                                                    <  1:456963/250‑1 (MQ=255)
                                                     tatCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCCTCGGCAGGCCAATAGTGCTTTCCCATTCCCCCCTTTTAAtt                                                                                                                                                                                            >  2:483323/1‑251 (MQ=255)
                                                        cgcgGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAAt                                                                                                                                                                                         <  1:287755/251‑1 (MQ=255)
                                                         gcgGGATGATCAATAAGATGAGTGCGGATTTGCCCATCTTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCCCTTTCAATTAATg                                                                                                                                                                                        <  2:584089/251‑1 (MQ=255)
                                                           gggATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGcc                                                                                                                                                                                      >  1:670772/1‑251 (MQ=255)
                                                            ggATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCg                                                                                                                                                                                     >  1:485249/1‑251 (MQ=255)
                                                              atgatCAAGAAGATGAGTGCGGGTTTGCCTATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGGCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTGCAGCAACTGGTGTGTATCGCGCACCATCTGTGGATTGCCGCTCAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAAGCGTGCTTTCCAGGTCCTCACGTTCAATTAATGCCgtt                                                                                                                                                                                   <  2:51807/251‑3 (MQ=255)
                                                                 atCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACtt                                                                                                                                                                                                  >  1:123464/1‑233 (MQ=255)
                                                                 atCAAAAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACtt                                                                                                                                                                                                  <  2:123464/233‑1 (MQ=255)
                                                                      tAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTAtccgtcc                                                                                                                                                                          >  1:79489/1‑251 (MQ=255)
                                                                               tGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGCCGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCTCTCCATCTGTGGATTGCCGCACAGCATCCCATAGCTGGTATCTTTATTCATCGGCAGGCCCATGGGGCTTTCCAGTTCCCCCCTTTTCATTACTGCCGTGATCCGACCGGTggacgg                                                                                                                                                                 >  2:24219/1‑246 (MQ=255)
                                                                                   gatttcatCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGAAATGCTCCGCTGTCATATGGCCCGGTAGGCGACGTAAATGTTTCGTCATCTGACGGGTCTCTTTCAGCCACTGTTGTGTATAGCGCAACCTCTGTGGATTGCCGCACAGCATCACCTAGCTGGTTTCTTTATTCATCGGGAGGACAATCGTGCTTTACAGTTTCCCACTTTCAATTAATTCAGGTATCCGTTCGCGGAGCGAcccc                                                                                                                                                             >  2:40518/8‑251 (MQ=255)
                                                                                              cGTGCCGTTTATCGATAAGTATCGGCTTCCTATTAATCTGCCGCTTCCCTATGTCCCTGATGGCGACGGAAATGTTTCGTCCTTTGCCGGGTCTCTTTCAGCACCTGTTGTGTATCGCGCACCATCTGTGGAGTGCCGCACAGCATCACATATCTGTTTTCTGTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTCCAATTAATGCCGGTATCCGTCCGGTGAGCGCCCCCGCTGCCGTTTc                                                                                                                                                  <  2:73583/251‑1 (MQ=255)
                                                                                               gggCCGTTTATCGATAAGTAACCGTTTACCAGTATCTCGTCGCTGTCGTATTGCCCGGTCGGCGACGTAAACGGTTCGTCATCTGCCGGGTCTCTTTCAGTAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACATCATCTCATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCGTTCCATTTCCCCACTTTTAATTAATGCCGGTATCCGGTCTATTAGCGTCCCCGCTGCCGTTTcc                                                                                                                                                 <  2:461534/249‑1 (MQ=255)
                                                                                                   cGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGAc                                                                                                                                             >  1:352593/1‑251 (MQ=255)
                                                                                                        aTCGATAAGTAACCGCTTACCAGCAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCCGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTTacc                                                                                                                                        >  1:30027/1‑251 (MQ=255)
                                                                                                            aTAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCATCTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCg                                                                                                                                    <  1:53161/251‑1 (MQ=255)
                                                                                                                    ccGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATg                                                                                                                            >  1:118525/1‑251 (MQ=255)
                                                                                                                        ttACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTAATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATgcgc                                                                                                                         >  1:454335/1‑250 (MQ=255)
                                                                                                                            cAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGTCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGttt                                                                                                                    <  1:529130/251‑1 (MQ=255)
                                                                                                                            cAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGttt                                                                                                                    >  1:67570/1‑251 (MQ=255)
                                                                                                                             cgTAATGCTCCGCTGTCATATGGCTCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGAGTATCGCGCACCATCTGTGGATTGCCGCACAGCATAACATAGCTGGTTTCTTTATTAATCGGCCGGCAAATCAGGTTTTCCAGTGCCCCCCTTTCAATTAATGCCGTACTCCGTCCGATTAGCGACCACGCTTCCGTTTCCCGAATGCCCACCgcattgatggcactattt                                                                                                                   >  2:443177/2‑234 (MQ=255)
                                                                                                                             attAATGCCCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGGCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGTCCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGtttt                                                                                                                   <  2:288412/249‑1 (MQ=255)
                                                                                                                             aGTAATGCTCCGCTGTCATATGGCCCGTTCGGGGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGtttt                                                                                                                   <  2:412224/251‑1 (MQ=255)
                                                                                                                                     tCCGCTGTCATATGGCCCAGTCGGCGACGTCAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGt                                                                                                            >  1:196235/1‑250 (MQ=255)
                                                                                                                                        cGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGc                                                                                                         <  1:71656/251‑1 (MQ=255)
                                                                                                                                                        ggTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTg                                                                                         >  1:605290/1‑251 (MQ=255)
                                                                                                                                                             gCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCa                                                                                    <  1:26806/251‑1 (MQ=255)
                                                                                                                                                                aCGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCTCTTTCAATTAATGCCGGTATCCTTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTg                                                                                 <  1:156941/251‑1 (MQ=255)
                                                                                                                                                                 cGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTgg                                                                                <  1:405111/251‑1 (MQ=255)
                                                                                                                                                                 cGTAAATGTTTCGTCATCTGCCGGGTATCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTCGCGTTTTTCCAGTTCCTGCATCAGTgg                                                                                >  1:468674/1‑251 (MQ=255)
                                                                                                                                                                     aaTGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCaa                                                                             >  1:137033/1‑250 (MQ=255)
                                                                                                                                                                         gttgccccccTGCCTGGTCTCTTTCAGCAACTGTTGGGTATCTCGCACCATCTGTGGATTGCCGAACAGCGCCACATAGCTGGTTTCTTTATTCATCGGCTGGCCAATCGTGCTTTCCAGGTCCCCACTTTCAATTAATGCCGGGATCCGTCCGGTTAGCGACCCCGCTGCTGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGGGGCAAATAGc                                                                        <  2:610034/242‑1 (MQ=255)
                                                                                                                                                                                  cTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCAGACTGACCACCGTCTGAATGCGCAGTTATCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGAcgg                                                               >  1:566398/1‑251 (MQ=255)
                                                                                                                                                                                         gTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTAtccgtcc                                                                                                                                                                          <  2:7694/137‑1 (MQ=255)
                                                                                                                                                                                         gTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTAtccgtcc                                                                                                                                                                          >  1:7694/1‑137 (MQ=255)
                                                                                                                                                                                         gTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATa                                                        <  1:483323/251‑1 (MQ=255)
                                                                                                                                                                                            tctTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACg                                                     >  1:458465/1‑251 (MQ=255)
                                                                                                                                                                                             ctTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGa                                                                                                                               <  2:137900/176‑1 (MQ=255)
                                                                                                                                                                                             ctTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGa                                                                                                                               >  1:137900/1‑176 (MQ=255)
                                                                                                                                                                                               ttCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCCCAGCATCACATAGCTGGTTTCTTTATTCAACGGCAGGCCAATCGTGCTTTCCAGTTCCCCTCTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTg                                                   <  2:629324/250‑1 (MQ=255)
                                                                                                                                                                                               ttCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTgc                                                  >  1:262441/1‑251 (MQ=255)
                                                                                                                                                                                                   gCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTTAGCGCCCCCGCTGCCGGTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATTGCTTAAGTCGGCGGCATAACGGgcggcg                                              >  2:346487/1‑251 (MQ=255)
                                                                                                                                                                                                     aaCTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTg                                            <  1:396054/251‑1 (MQ=255)
                                                                                                                                                                                                      aCTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCCATTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGCCCACCGTCTGAATGCGCCGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCTTAACGTGCGGCGTGc                                           >  2:32362/1‑251 (MQ=255)
                                                                                                                                                                                                       cTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCCTCCGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGcc                                          >  2:185956/1‑250 (MQ=255)
                                                                                                                                                                                                         gttgtGTCTCCCGCACCATCTGTGCATTGCACCACAGCATCACATAGCTTGTTTCTTTATTCATCTGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCACTTAATGCCGGTATCCTTCCGGTGAGCGACCCCGCTGCCTTTTCCCGACTGACCACCGTCTGAATGCGCACTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATTAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCTGCGTGCAcc                                        <  2:171449/251‑1 (MQ=255)
                                                                                                                                                                                                         gttgtGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCCGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCAcc                                        >  2:275020/1‑251 (MQ=255)
                                                                                                                                                                                                          ttgtGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCaa                                                                                                                                                                                              >  2:203448/1‑101 (MQ=255)
                                                                                                                                                                                                          ttgtGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCaa                                                                                                                                                                                              <  1:203448/101‑1 (MQ=255)
                                                                                                                                                                                                              gtATCGCGCACCATCTGTGGATGGCCGCACAGCATCACATAGCTGGTTGCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGAc                                   <  2:567598/251‑1 (MQ=255)
                                                                                                                                                                                                                aTCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGAcc                                  >  1:106765/1‑250 (MQ=255)
                                                                                                                                                                                                                     gcACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGAttt                            <  1:24219/251‑1 (MQ=255)
                                                                                                                                                                                                                       aCCATCTGTGGATTGCCGCACCGCATCACATAGCTGTTTTCTTTATTCATCGGCAGGCCAATCGTTCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCTCCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGAttttt                          <  1:519067/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gATTTCCGCACAGCATCTCATAGCTGGATTCTTTATTCATCGGCAGGCCAATCGTTCTTTCCAGTTCCCCTCTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCTCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGTTCt                <  2:485249/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ttGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGTCCAATCGTGCTTTCCAGTTCCCCACTTTCACTTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGCCTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGAGGCCAATAGCTTAAGTCGGCTGCATAACGTGCGGCTTGCACCAGGCCCAGATCTTTGAAGCGATCTa               >  2:68194/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ccGCACAGCATCCCATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTGTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGGTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGAGCTAAATc           <  2:67570/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCCGATTTTTGAAGCGATCTAAATCTTTACCt    >  2:594856/1‑249 (MQ=255)
                                                                                                                                                                                                                                               aTCACATAGCTGGTTTCTTTATTCAACGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGGTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTAg  <  2:106765/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ACCTCCTTGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTAG  >  NC_000913/4113606‑4114094

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