Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,003,897 C→A 6.6% R39S (CGT→AGT) ‡ allE → S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,003,8970CA 6.6% 655.1 / 25.7 198R39S (CGT→AGT) ‡allES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base C (91/94);  new base A (8/5);  total (99/99)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 199 total aligned reads displayed.

GCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCC  >  NC_002947/4003750‑4004046
                                                                                                                                                   |                                                                                                                                                     
gCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGt                                                                                                                                                     <  2:665439/150‑1 (MQ=255)
gCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGt                                                                                                                                                     <  1:2026604/150‑1 (MQ=255)
   tGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGAc                                                                                                                                                  <  1:3834628/150‑1 (MQ=255)
    ggggAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACa                                                                                                                                                 >  1:4070382/1‑150 (MQ=255)
      ggAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATc                                                                                                                                               >  1:3845615/1‑150 (MQ=255)
      ggAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTCCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATc                                                                                                                                               <  2:3694299/150‑1 (MQ=255)
              caacaaGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAg                                                                                                                                       >  1:3454503/1‑150 (MQ=255)
               aacaaGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGc                                                                                                                                      >  1:3110517/1‑150 (MQ=255)
                 caaGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCa                                                                                                                                    >  1:2050816/1‑150 (MQ=255)
                  aaGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCAc                                                                                                                                   >  2:3609734/1‑150 (MQ=255)
                    gTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCt                                                                                                                                 <  1:3987448/150‑1 (MQ=255)
                    gTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCt                                                                                                                                 <  2:2757169/150‑1 (MQ=255)
                      cGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTcacctgccacctgc                                                                                                                               >  2:4965157/1‑150 (MQ=255)
                       gagaTCTTCCATGTAGAACCAATCATACTTCTCCCCCCACGGTGGGAACCCGGCTCAAACAGAGCTACTGACTTACCGTGCCATGTTCACAGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGAAATAGTCAACAGCAACCAGcc                                                                                                                              >  2:4051780/1‑150 (MQ=255)
                       gagaTCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGcc                                                                                                                              >  1:4451819/1‑150 (MQ=255)
                          aTCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCtt                                                                                                                           <  1:1287992/150‑1 (MQ=255)
                             gTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCGCTTCTg                                                                                                                        <  2:121754/150‑1 (MQ=255)
                              tCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTgg                                                                                                                       <  1:4054427/150‑1 (MQ=255)
                              tCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTgg                                                                                                                       <  2:2575144/150‑1 (MQ=255)
                              tCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTgg                                                                                                                       <  1:2959255/150‑1 (MQ=255)
                              tCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTgg                                                                                                                       <  1:1244575/150‑1 (MQ=255)
                                cATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGa                                                                                                                     <  1:4873112/150‑1 (MQ=255)
                                       aaCCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATg                                                                                                              >  2:1396312/1‑150 (MQ=255)
                                          cATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATgcgc                                                                                                           <  1:4719703/150‑1 (MQ=255)
                                            ttCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCa                                                                                                                                        >  2:4265271/1‑119 (MQ=255)
                                            ttCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAt                                                                                                         <  2:2018834/150‑1 (MQ=255)
                                              cATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgt                                                                                                       <  2:3911073/150‑1 (MQ=255)
                                               aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgtg                                                                                                      >  2:3054397/1‑150 (MQ=255)
                                               aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgtg                                                                                                      >  1:2789145/1‑150 (MQ=255)
                                               aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtggg                                                                                                      >  2:4015401/1‑148 (MQ=255)
                                                tACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTgg                                                                                                     <  2:5197552/150‑1 (MQ=255)
                                                tACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTgg                                                                                                     <  1:1268723/150‑1 (MQ=255)
                                                 aCTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTggg                                                                                                    <  1:3054397/150‑1 (MQ=255)
                                                 aCTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTggg                                                                                                    <  2:1271128/150‑1 (MQ=255)
                                                   ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTgg                                                                                                                       <  1:2394365/129‑1 (MQ=255)
                                                   ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTgg                                                                                                                       >  2:2394365/1‑129 (MQ=255)
                                                   ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                  >  2:3836936/1‑150 (MQ=255)
                                                   ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                  <  1:1555934/150‑1 (MQ=255)
                                                   ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                  >  2:1987956/1‑150 (MQ=255)
                                                   ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                  >  1:4627158/1‑150 (MQ=255)
                                                   ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                  >  1:3153213/1‑150 (MQ=255)
                                                   ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                  >  1:841712/1‑150 (MQ=255)
                                                    tCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTa                                                                                                 <  2:2533162/150‑1 (MQ=255)
                                                      gCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATc                                                                                               <  2:3297606/150‑1 (MQ=255)
                                                      gCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATc                                                                                               <  1:4965157/150‑1 (MQ=255)
                                                        cccccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgc                                                                                             <  1:1582951/150‑1 (MQ=255)
                                                        cccccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgc                                                                                             <  1:1606906/150‑1 (MQ=255)
                                                          cccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgccc                                                                                           <  1:427785/150‑1 (MQ=255)
                                                            cACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgcccgc                                                                                         >  2:3926111/1‑150 (MQ=255)
                                                                   ggCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGa                                                                                  <  2:4503876/150‑1 (MQ=255)
                                                                    gCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACaa                                                                                                                  >  2:3988410/1‑117 (MQ=255)
                                                                      aCCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCg                                                                               >  2:3681589/1‑150 (MQ=255)
                                                                        ccGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGt                                                                             >  2:5229861/1‑150 (MQ=255)
                                                                          ggCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGttt                                                                           <  1:1173233/150‑1 (MQ=255)
                                                                          ggCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGttt                                                                           >  1:3475071/1‑150 (MQ=255)
                                                                             tCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGc                                                                        <  1:4371376/150‑1 (MQ=255)
                                                                             tCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGc                                                                        >  1:4526266/1‑150 (MQ=255)
                                                                              cAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGcc                                                                       <  1:828689/150‑1 (MQ=255)
                                                                                  ccGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGaaa                                                                   <  2:1342641/150‑1 (MQ=255)
                                                                                   cGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAAc                                                                  <  1:291817/150‑1 (MQ=255)
                                                                                   cGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAAc                                                                  >  1:3881331/1‑150 (MQ=255)
                                                                                   cGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAAc                                                                  >  1:1060396/1‑150 (MQ=255)
                                                                                      gtcgatgactgacCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGtt                                                               <  1:5079356/145‑1 (MQ=255)
                                                                                      gctgctGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGtt                                                               <  1:590523/150‑1 (MQ=255)
                                                                                       ctcctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtc                                                              <  1:3749418/147‑1 (MQ=255)
                                                                                        tcctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtct                                                             >  1:861330/3‑150 (MQ=255)
                                                                                         gctgACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTtctc                                                            >  1:2726046/1‑150 (MQ=255)
                                                                                          ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                           <  1:4612988/150‑1 (MQ=255)
                                                                                          ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                           <  1:3341939/150‑1 (MQ=255)
                                                                                          ctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcg                                                           >  2:3070560/1‑149 (MQ=255)
                                                                                           tgaccgatCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcgc                                                          >  2:174779/1‑148 (MQ=255)
                                                                                            gactgacCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCa                                                         <  2:3153213/150‑1 (MQ=255)
                                                                                                  ccGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAt                                                   <  2:3115479/150‑1 (MQ=255)
                                                                                                   cGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtc                                                  >  2:3831530/1‑150 (MQ=255)
                                                                                                     tGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcgt                                                <  1:4559789/150‑1 (MQ=255)
                                                                                                      gCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatg                                               <  1:739164/150‑5 (MQ=255)
                                                                                                      gCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatg                                               <  2:3466720/150‑5 (MQ=255)
                                                                                                      gCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcgtc                                               <  2:1931985/150‑1 (MQ=255)
                                                                                                       ccATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatgg                                              <  1:2254002/150‑6 (MQ=255)
                                                                                                         aTGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGaa                                            >  2:920635/1‑150 (MQ=255)
                                                                                                          tGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTGCATCGTCGAAc                                           >  1:2627360/1‑150 (MQ=255)
                                                                                                            ttCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTg                                         >  2:5369028/1‑150 (MQ=255)
                                                                                                               aCCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGcc                                      >  1:4821180/1‑150 (MQ=255)
                                                                                                                        tATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggc                             >  2:2835337/1‑150 (MQ=255)
                                                                                                                             cGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCag                        <  1:1354566/150‑1 (MQ=255)
                                                                                                                              gTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagc                       >  2:3035391/1‑150 (MQ=255)
                                                                                                                              gTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagc                       <  1:3584146/150‑1 (MQ=255)
                                                                                                                                cATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcga                     <  1:5227830/150‑1 (MQ=255)
                                                                                                                                 aTCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCa                                                         <  1:635896/113‑1 (MQ=255)
                                                                                                                                 aTCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCa                                                         >  2:635896/1‑113 (MQ=255)
                                                                                                                                  tCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcgagc                   >  2:5038068/1‑150 (MQ=255)
                                                                                                                                   ccccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcgagcg                  <  2:4821180/150‑1 (MQ=255)
                                                                                                                                     ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                >  2:1086081/1‑150 (MQ=255)
                                                                                                                                      cAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGccc               >  1:1421481/1‑150 (MQ=255)
                                                                                                                                      cAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGccc               >  2:3522754/1‑150 (MQ=255)
                                                                                                                                      cAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCGGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGccc               >  2:4406907/1‑150 (MQ=255)
                                                                                                                                       aaaGGCGTGATGCGTGACATCGTCAACAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCa              <  2:4160265/150‑1 (MQ=255)
                                                                                                                                             gTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTg        <  1:938432/150‑1 (MQ=255)
                                                                                                                                              tGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTgg       >  1:5322214/1‑150 (MQ=255)
                                                                                                                                               gATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATAGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTATCCCAGTACATCGACGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGa      >  1:4710504/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCC  >  NC_002947/4003750‑4004046

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