Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,003,920 G→T 7.7% L46L (CTG→CTT allE → S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,003,9200GT 7.7% 595.3 / 28.6 182L46L (CTG→CTTallES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (97/71);  new base T (9/5);  total (106/76)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 185 total aligned reads displayed.

TCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAT  >  NC_002947/4003771‑4004069
                                                                                                                                                     |                                                                                                                                                     
    gaTCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCt                                                                                                                                                   >  2:2984819/1‑150 (MQ=255)
      tCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTc                                                                                                                                                 >  2:4921788/1‑150 (MQ=255)
        gTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTg                                                                                                                                               <  2:2210649/150‑1 (MQ=255)
         tCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTgg                                                                                                                                              >  2:2553808/1‑150 (MQ=255)
              gTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACaa                                                                                                                                         >  2:1865318/1‑150 (MQ=255)
                 gAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGAt                                                                                                                                      >  1:3570654/1‑150 (MQ=255)
                  aaCCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATg                                                                                                                                     >  2:1720952/1‑150 (MQ=255)
                  aaCCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATg                                                                                                                                     >  2:1721167/1‑150 (MQ=255)
                   aCCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATgc                                                                                                                                    >  1:1508700/1‑150 (MQ=255)
                    ccATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATgcg                                                                                                                                   >  2:42928/1‑150 (MQ=255)
                     cATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGTCTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATgcgc                                                                                                                                  <  2:44973/150‑1 (MQ=255)
                       ttCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAt                                                                                                                                >  1:3388725/1‑150 (MQ=255)
                        tCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtg                                                                                                                               >  1:4642582/1‑150 (MQ=255)
                         cATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgt                                                                                                                              >  2:3560100/1‑150 (MQ=255)
                          aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgtg                                                                                                                             >  2:3467228/1‑150 (MQ=255)
                           tACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTgg                                                                                                                            <  2:1151549/150‑1 (MQ=255)
                           aaCTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTgg                                                                                                                            <  2:2829910/149‑1 (MQ=255)
                              ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                         >  2:1232953/1‑150 (MQ=255)
                              ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                         >  2:307162/1‑150 (MQ=255)
                              ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                         <  2:460045/150‑1 (MQ=255)
                                 gCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATc                                                                                                                      <  2:3801571/150‑1 (MQ=255)
                                 gCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATc                                                                                                                      >  2:3406472/1‑150 (MQ=255)
                                  ccccccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcg                                                                                                                     >  2:1933113/1‑150 (MQ=255)
                                     cccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgccc                                                                                                                  >  2:3140630/1‑150 (MQ=255)
                                     cccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgccc                                                                                                                  >  2:4709194/1‑150 (MQ=255)
                                      ccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgcccg                                                                                                                 >  2:3829381/1‑150 (MQ=255)
                                       cACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgcccgc                                                                                                                <  1:4131207/150‑1 (MQ=255)
                                           gTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCccgcccgc                                                                                                            <  1:1158952/150‑1 (MQ=255)
                                            tGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCt                                                                                                           >  2:1315838/1‑150 (MQ=255)
                                            tGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCt                                                                                                           <  2:4387421/150‑1 (MQ=255)
                                               gCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGAc                                                                                                        <  1:1158131/150‑1 (MQ=255)
                                                   ccGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGt                                                                                                    <  1:2031496/150‑1 (MQ=255)
                                                   ccGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGt                                                                                                    >  2:445107/1‑150 (MQ=255)
                                                      gCTCAAACCGAGCTGCTGACTAACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGtttt                                                                                                 <  2:5029107/150‑1 (MQ=255)
                                                       cTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTg                                                                                                <  2:2083002/150‑1 (MQ=255)
                                                        tCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGc                                                                                               >  2:536406/1‑150 (MQ=255)
                                                        tCAAACCGAGCTGCTGACTGACCGTGCCATGTTAACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTCACATCGTCACCAGCCACCTGCACTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGc                                                                                               >  2:2641609/1‑150 (MQ=255)
                                                           aaCCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGa                                                                                            >  2:3007008/1‑150 (MQ=255)
                                                            aCCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGaa                                                                                           >  1:4034274/1‑150 (MQ=255)
                                                                 actcctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACttt                                                                                      <  1:3411106/146‑1 (MQ=255)
                                                                  ctgctgACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTtc                                                                                     <  2:2489910/150‑1 (MQ=255)
                                                                   tgctgACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTtct                                                                                    >  1:4619820/1‑150 (MQ=255)
                                                                   tgctgACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTtct                                                                                    <  1:4921788/150‑1 (MQ=255)
                                                                   tgctgACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTtct                                                                                    >  1:1729427/1‑150 (MQ=255)
                                                                    gctgACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTtctc                                                                                   >  2:629676/1‑150 (MQ=255)
                                                                     ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                  <  2:748429/150‑1 (MQ=255)
                                                                     ctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcg                                                                                  >  2:3235826/1‑149 (MQ=255)
                                                                      tgaccgatCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcgc                                                                                 >  1:1393565/1‑148 (MQ=255)
                                                                        actgacCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGACGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAg                                                                               >  1:1188148/1‑150 (MQ=255)
                                                                          tgacCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTa                                                                             >  2:4825433/1‑150 (MQ=255)
                                                                          tgacCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTa                                                                             >  1:2710222/1‑150 (MQ=255)
                                                                            acCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACa                                                                           <  1:4509426/150‑1 (MQ=255)
                                                                               gcgCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatatta                                                                        <  2:2495698/148‑3 (MQ=255)
                                                                               gTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcg                                                                        <  1:1755983/150‑1 (MQ=255)
                                                                                 gCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatg                                                                      <  2:1041160/150‑5 (MQ=255)
                                                                                 gCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatg                                                                      >  2:1629832/1‑146 (MQ=255)
                                                                                 gCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatg                                                                      >  2:1389616/1‑146 (MQ=255)
                                                                                 gCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcgtc                                                                      <  2:1468286/150‑1 (MQ=255)
                                                                                 gCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAGGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGGACAtcgtc                                                                      >  2:4249419/1‑150 (MQ=255)
                                                                                      gTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACt                                                                 >  2:3669666/1‑150 (MQ=255)
                                                                                       ttCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTg                                                                >  2:3990228/1‑150 (MQ=255)
                                                                                         cACCGAAGCCTATGCCGTCATCCCCAACGGCGTGATGCGTGACCTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGc                                                              <  1:1267036/150‑1 (MQ=255)
                                                                                         cACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGc                                                              <  2:1513868/150‑1 (MQ=255)
                                                                                          aCCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgc                                                                                                                    >  1:49787/1‑95 (MQ=255)
                                                                                           ccGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGccc                                                            <  2:3755584/150‑1 (MQ=255)
                                                                                            cGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGcccc                                                           >  1:1377804/1‑150 (MQ=255)
                                                                                            cGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGcccc                                                           >  1:1006774/1‑150 (MQ=255)
                                                                                             gAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCg                                                          <  2:1913504/150‑1 (MQ=255)
                                                                                             gAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCg                                                          >  1:3399080/1‑150 (MQ=255)
                                                                                             gAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCg                                                          >  1:827064/1‑150 (MQ=255)
                                                                                              aaGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCgg                                                         >  1:1065803/1‑150 (MQ=255)
                                                                                                gCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGaa                                                       <  1:994288/150‑1 (MQ=255)
                                                                                                 ccTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAg                                                      >  2:1768455/1‑150 (MQ=255)
                                                                                                   tATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACCAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggc                                                    >  1:3604700/1‑150 (MQ=255)
                                                                                                   tATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGCTGCGCCTGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGACACGTTCTCCCAGTACCTCGTCGCACTGGCCCCGGAAggc                                                    >  2:4767332/1‑150 (MQ=255)
                                                                                                     tGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggcgg                                                  <  1:5042728/150‑1 (MQ=255)
                                                                                                      gCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggcggc                                                 <  2:1932025/150‑1 (MQ=255)
                                                                                                       ccGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCa                                                <  1:2862685/150‑1 (MQ=255)
                                                                                                       ccGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCa                                                <  1:2756525/150‑1 (MQ=255)
                                                                                                               cccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAgcgc                                        >  1:129890/1‑150 (MQ=255)
                                                                                                                ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                       >  1:3161221/1‑150 (MQ=255)
                                                                                                                    aGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCaga                                   <  2:3656677/150‑1 (MQ=255)
                                                                                                                        gTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTg                               >  1:2983776/1‑150 (MQ=255)
                                                                                                                           aTGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGAc                            <  1:4297842/150‑1 (MQ=255)
                                                                                                                                  aCATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATg                     <  1:4911491/150‑1 (MQ=255)
                                                                                                                                   cATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGc                    <  2:2878853/150‑1 (MQ=255)
                                                                                                                                        tCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAg               >  1:4648118/1‑150 (MQ=255)
                                                                                                                                            cAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgt           >  1:3996444/1‑150 (MQ=255)
                                                                                                                                            cAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgt           <  1:4843962/150‑1 (MQ=255)
                                                                                                                                            cAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgt           <  2:3388725/150‑1 (MQ=255)
                                                                                                                                            cAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgt           >  1:4414544/1‑150 (MQ=255)
                                                                                                                                             aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGATACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg          <  1:643758/150‑1 (MQ=255)
                                                                                                                                             aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg          >  2:178796/1‑150 (MQ=255)
                                                                                                                                             aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg          <  1:2501669/150‑1 (MQ=255)
                                                                                                                                             aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg          <  1:3560100/150‑1 (MQ=255)
                                                                                                                                              gCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtgg         <  1:2677489/150‑1 (MQ=255)
                                                                                                                                                 aCCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGgtgt      >  1:1869513/1‑150 (MQ=255)
                                                                                                                                                 aCCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGgtgt      <  1:1777536/150‑1 (MQ=255)
                                                                                                                                                 aCCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGgtgt      >  1:1869542/1‑150 (MQ=255)
                                                                                                                                                     gCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAt  >  1:4796740/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAT  >  NC_002947/4003771‑4004069

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: