Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 968,837 G→A D158N (GAT→AAT)  lpxK → tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913968,8370GA97.3% 237.1 / NA 74D158N (GAT→AAT) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (26/46);  total (27/47)
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

ATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGA  >  NC_000913/968594‑969077
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
aTTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGGTTAAACCGGGTTGATTTTTACAAGCAATTATGGGCCTGGTTCGGTTTCTACCGTTTGTTTTTATGCGGTAAAAACCATTCTGGCGGCACACCCACATGTGCCCATCATCGTAAACGGCGACGGGTTTCAGCATTTCCGCCTGGGGCGGGATGTggg                                                                                                                                                                                                                                              >  1:681972/1‑250 (MQ=255)
   cgcgTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCCTTCTGGCGCCACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAACATTACCCTCTGGCGCGTAATGTGGAAAt                                                                                                                                                                                                                                          >  1:108849/1‑251 (MQ=255)
    gcgTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCCCTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCCGGTGTGCAGATCATCGTAACCGACGACGGGTTTCAGCATTACCGTCTGGGGCGGAATGTGGGAAtt                                                                                                                                                                                                                                         >  1:697763/1‑251 (MQ=255)
            gTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTAt                                                                                                                                                                                                                                 <  1:17287/251‑1 (MQ=255)
               gTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGCTGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGCCGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTgt                                                                                                                                                                                                                              >  1:411761/1‑250 (MQ=255)
                 aTCGCGGGGATATGGTGGTAAGGCTGAAGCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGAt                                                                                                                                                                                                                             <  2:286558/250‑1 (MQ=255)
                      ggggagagggtggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGtgtg                                                                                                                                                                                                                       <  2:287947/243‑1 (MQ=255)
                          atagggtggtAAGGCTGAATCTTAGCCGCTGTTATTGTTGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTc                                                                                                                                                                                                                   <  2:582914/247‑1 (MQ=255)
                                 ggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAAACTGTGTTGATTTATCAACGCACTGATGCGCCTGTTTCGGGTTCTCCCGGTCGTTCTGATGCGGTAAAAACCATTCTTGCGCCACACCCTGATGTGCAGATCATCGTAACCGCCGACGGGTTACAGAATTACCGTCTGGCGCGGAATGTGGAAATTGTCGTTTTTTGTGGTGTGGGTCGCGTTgg                                                                                                                                                                                                            >  1:671448/1‑251 (MQ=255)
                                           aaTCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATTAACGCCCCGATGCGCCTGGTTCGGTTTCTCCCGTTCCTTCTGATGCGGTAAAAGCCATTTTGGCGCAACACCCCTGTTTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTCGTTTTTGATGGTGTGGGTCGCTTTGGCAATGGCtgg                                                                                                                                                                                                  >  1:160538/1‑251 (MQ=255)
                                           aaTCTTATCCGCTGTTATTGTCGGCAGATAACACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCtg                                                                                                                                                                                                   <  1:213886/250‑1 (MQ=255)
                                                tATCCGTTGTTATTGTCTGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGTCGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCtggtg                                                                                                                                                                                                <  1:791082/248‑1 (MQ=255)
                                                tATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCCTCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCtggtg                                                                                                                                                                                                <  1:812065/248‑1 (MQ=255)
                                                  tCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTTTCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACCCCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGGTATTGATGGGGTGCGTCGCTTTGGCAATGGCTGGGGGTTGc                                                                                                                                                                                          >  1:609775/1‑251 (MQ=255)
                                                    cGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAACCTGTGTTGATTTATTAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTTATGCGGTAAAAACCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTTTTTATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCcg                                                                                                                                                                                        >  1:246919/1‑251 (MQ=255)
                                                     gCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCcgg                                                                                                                                                                                       <  2:411761/251‑1 (MQ=255)
                                                        gTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCcggc                                                                                                                                                                                      <  1:106789/249‑1 (MQ=255)
                                                         ttATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGGACCTGTGTTGATTTATCAACGCAATGATGCGCCTGTTGGGGTTTCTCCCGTTCGTTCTGATTCGGGAAAAGCCATTCTGGCGCCACACCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTACCGTCTGGCGGGTAATGTGGAAATTGTCGTTATTTGTGGTGTGCGTCGCTTTGGCAATGGCTGGTGGGTGCCGGCggg                                                                                                                                                                                   >  1:237586/1‑251 (MQ=255)
                                                           aTTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGGTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTAACGTCTGGCGCGTAATGTGGGAATTGTCGTTATTGGTGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGGCGggggggg                                                                                                                                                                                 >  1:623745/1‑250 (MQ=255)
                                                                   cAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGt                                                                                                                                                                        <  1:466458/251‑1 (MQ=255)
                                                                    aGATACCACAACAGCACAGGCGGGTGGTGAACCTGTGTTGATTTTTCAAAGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGCAAAAGCCCTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAACGACGAAGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTTGCGGCGGGGCCacagcgtg                                                                                                                                                                       >  1:57273/1‑244 (MQ=255)
                                                                    aGATACCACAACAGCACAGGCGGGTGATGAAACTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATTTGCAGATCCTCGTAACCGACGACGGTTTAAAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTTATGGGGTGCGGCGCTTTTGCAATGGCTGGGGGGTTGCGGCgggggccctgcggg                                                                                                                                                                        >  1:611064/1‑241 (MQ=255)
                                                                     gATACCACAACAGCACAGGCGGGTGATGAACCTGGGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTga                                                                                                                                                                      <  1:229627/251‑1 (MQ=255)
                                                                        aCCACAACAGCACAGGCGGGTGATGAACCTGTGTTGGTTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTTTGATGCGGGAAAAGCCCTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGGGTGCGTCGCTTTGGCAATGGCTGGGGGTTGCCGGCGGGGGCCATGCGTgaggg                                                                                                                                                                   >  1:85045/1‑249 (MQ=255)
                                                                                aGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCGGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGgc                                                                                                                                                           <  2:770209/251‑1 (MQ=255)
                                                                                 gCACAGGCGGGGGATGAACCTGTGTTGATTTATTAAAGCACTGATGCGCCTGTTGCGGTTTTTCCCGTTCGTTTTGATGCGGGAAAAACCCTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAACGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTCGTTATTGTTGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGGCGGCGGGGGCCAAtcgtgggcgggggggggg                                                                                                                                                           >  1:352474/1‑234 (MQ=255)
                                                                                     aGGCGGGTGATGAACCTGTGTTGATTTAGCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCGGATGCGGTAAAAGCCATTCTGGCGCAACTCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTAGTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTa                                                                                                                                                      <  2:101596/251‑1 (MQ=255)
                                                                                      ggcggggtatgaACCTGTGTTGATTTACCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTTATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGGCGTTAGTGATGGTGTGCGTCGCTGTGGCAATGGCTTGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTaa                                                                                                                                                     <  2:761707/243‑1 (MQ=255)
                                                                                              atgaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGtt                                                                                                                                             <  1:415623/251‑1 (MQ=255)
                                                                                                gaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTGCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGAGGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGGGCGTCGCTGTGGCAATGGCTGGGGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGCGGTTg                                                                                                                                            <  1:96816/250‑1 (MQ=255)
                                                                                                gaACCTGTGTTGATTGATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTg                                                                                                                                            <  1:314774/250‑1 (MQ=255)
                                                                                                  aCCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAAAGGGTGAGCGGGGGGGGGGCGTTAAGTCgggttgtg                                                                                                                                         >  2:37286/1‑245 (MQ=255)
                                                                                                   ccTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACCCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGAGGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGCGGTTGATGc                                                                                                                                        <  2:654222/251‑1 (MQ=255)
                                                                                                   ccTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGTTGc                                                                                                                                        >  1:21569/1‑251 (MQ=255)
                                                                                                   ccTGTGTTGATTTATCAACGCACTGATGAGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATAGTAACCGACGAAGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGAGTGCGACGCTTTGGCAAGGGCTGGAGGATGCCGGCGGGGCCAATACGTGAGCGAGCGGGGGGGTTaaacacggtggaggg                                                                                                                                        >  1:182218/1‑239 (MQ=255)
                                                                                                      gtggtGATTTATCAACGCACTGATGCGCCTGTTGCGGTTGCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGGCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAATTCGGTTGATGCGGt                                                                                                                                     <  2:737288/247‑1 (MQ=255)
                                                                                                               gttgcAACGCATTGGTGCGCCTGTTGCGGTTTCGCCCGGTCGTTCGGATGCGGTAAAAGCCATGCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCAGTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGAGGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTc                                                                                                                              <  2:375021/246‑1 (MQ=255)
                                                                                                                     aCGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcggcg                                                                                                                       >  1:129324/1‑251 (MQ=255)
                                                                                                                        cACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTAGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCCTGAGCGAGCGGGGCGCTTTAATGCGGGTGGTGCGGTAATCGTCacaggcgggg                                                                                                                    >  1:398506/1‑242 (MQ=255)
                                                                                                                               gcgcCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACCCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTc                                                                                                              <  1:614620/250‑1 (MQ=255)
                                                                                                                                cgacTGTTGCGGTTCCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGTGCATCGTAACCGACGACGGTTTACAGCATTACCGGCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGc                                                                                                            <  2:727760/248‑1 (MQ=255)
                                                                                                                                 gcCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGCCGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGgtgttccccccc                                                                                                           >  1:408186/1‑243 (MQ=255)
                                                                                                                                 gcCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTTAAGGCGGTTGATGCGGTAATCGTCAAcggcggggttccctgcc                                                                                                           >  1:812536/1‑240 (MQ=255)
                                                                                                                                     gTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGGGGGAATCGCCAACGGCGGGGTccgccggcgcgg                                                                                                       >  1:507352/1‑241 (MQ=255)
                                                                                                                                          ggTTTCTCCCGTTCGGTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCt                                                                                                               <  1:326765/238‑1 (MQ=255)
                                                                                                                                           gTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGTGAGCGGGCGGGGGGCTTAAAGGCGGTTGATGCGGGAATCGGCAAAGGCGGGGTCCCCCGCAACgggggaa                                                                                                  >  1:181255/1‑245 (MQ=255)
                                                                                                                                           gTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCAGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAAt                                                                                                 <  1:171409/251‑1 (MQ=255)
                                                                                                                                                 cccGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATcccc                                                                                             <  1:642899/249‑1 (MQ=255)
                                                                                                                                                  ccGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGGCCCTCGCAGCgggggaatacccctg                                                                                          >  1:703131/1‑238 (MQ=255)
                                                                                                                                                  ccGTGCGTTCGGATGCGGTAAAAGCCAGTCGGGCGCAACCCCCTGAGGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATg                                                                                          <  2:611064/251‑1 (MQ=255)
                                                                                                                                                          tcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAACGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATctgc                                                                                   <  1:240569/250‑1 (MQ=255)
                                                                                                                                                                       aaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCg                                                                     <  1:634610/251‑1 (MQ=255)
                                                                                                                                                                        aGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTAGTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCgg                                                                    <  2:697763/251‑1 (MQ=255)
                                                                                                                                                                                gggCGCAACACCCTGATGTGCAGAGCATCGTAACCGACGACGGTTGACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTa                                                            <  2:255132/250‑1 (MQ=255)
                                                                                                                                                                                     cAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGGGGGGCGCTTAAAGTCGGTTGATGAGGTAATCGTCAAAGGCGGGGTCCCTCCCAGCGGGGGAAACCCCATGCATCTGCTGCCgggggcggcgggggatttaagtac                                                       >  1:541062/1‑230 (MQ=255)
                                                                                                                                                                                         ccccTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGt                                                   <  2:85045/250‑1 (MQ=255)
                                                                                                                                                                                         ccccTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCAGGTCAGGCGGTGAATTTACGTACCGGt                                                   <  1:655316/250‑1 (MQ=255)
                                                                                                                                                                                         ccccTGAGGTGCATATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTTGCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGt                                                   <  2:681972/250‑1 (MQ=255)
                                                                                                                                                                                            ccttagtggcagactagcGTAACCGAAGACGGTTTACATCATTAACGTCTGGCGCGTAATGTCGTAAGTGGCGTTAGTTATGGTGTGCGTCGTTTTGGCAATGGTTGGTGGTTGCCGGCTGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGCCGGATGTTGCGGTAATCGTCAACGGCGGTGTCCCTCGCGACGGTCAAGTCCCCATTCATCTGCTGCCGGGGCAGGCTGTGAATTTACCTCCCGGttc                                                 <  2:379094/234‑3 (MQ=255)
                                                                                                                                                                                               gATGTGCAGATCAGCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTAcgcg                                              <  2:246919/251‑1 (MQ=255)
                                                                                                                                                                                                atgggCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGt                                             <  2:108849/247‑1 (MQ=255)
                                                                                                                                                                                                   tgCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTTCCCCGCAGCGGCGGAAACCCCACTCCTCTGCCGCCGGGTCCGGCGGGGGATTTACGGAACGGGGCGCGTtgt                                          >  1:769408/1‑251 (MQ=255)
                                                                                                                                                                                                                aaCCCACGAGGGATTACAGCATTACCGTCTGGCTCTGCATTGGAATATTGTCTTTATGGATGGTGTGCTTCTCGTTGGCAATGGCTGTTGGGTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAACTCGTGTGAAGCGGTAATCGTCAACGGCGTTTTCCCTCGCAGCGGTGAAATCCCCATGCAGCTGCTGCCGGGTCATGCGGTGAATTTTCGTACCGGTACGCGTTTTGACCTTCCTCAGc                             <  2:516962/251‑1 (MQ=255)
                                                                                                                                                                                                                 aCCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCt                            <  1:99/251‑1 (MQ=255)
                                                                                                                                                                                                                 aCCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCt                            <  2:160538/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTAGTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGAc                                       <  2:591919/239‑1 (MQ=255)
                                                                                                                                                                                                                         cgGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATg                    <  1:598487/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGGTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGt                   <  1:760433/251‑1 (MQ=255)
                                                                                                                                                                                                                              tACAGCATTACCGTCTGGCGCGGAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCTTTGTGACGTTGCTCAGCTTGAACATGTAGt                <  1:625060/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gCATTACCGTCTGGCGCGTAATGGGGAAATTGTCGTTTTTGATGGTGGGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGc             <  2:541062/249‑1 (MQ=255)
                                                                                                                                                                                                                                   cATTTCCGTCTGGCGCGTAATGTGGAAATTGGCGTTATTGAGGGTGTGCGTCGCTTTGGCAATGGCTTGTGGTTGCCGGCGGGGCCAATGCGTGCGCGAGCGGGGCGCTTAAAGGCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGAt          <  2:507352/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cggttGGCGCGTAATGTGGAAATTGTCGGTATTGATGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCggg    <  2:597127/247‑1 (MQ=255)
                                                                                                                                                                                                                                          gTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGGGGAATCCCCCTGCCTCTGCTGCCGGGGCCAGCGGTGGATTTACGTACCGGTACGCGTTGTGGCGTTGTTCAGCttttacaatttgtggcggttgcgggg   >  1:245846/1‑227 (MQ=255)
                                                                                                                                                                                                                                           gcTGGCTCGTAATGTGGAAATTGTCGTGATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGATAGCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAAGTTACGGACCTGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGa  <  1:475893/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
ATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGA  >  NC_000913/968594‑969077

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