Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,044,654 A→G 5.3% E173G (GAG→GGG)  PP_3560 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,044,6540AG 5.3% 1298.5 / 35.0 398E173G (GAG→GGG) PP_3560LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (190/187);  new base G (12/9);  total (202/196)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 402 total aligned reads displayed.

GGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGCCAGGAAG  >  NC_002947/4044505‑4044803
                                                                                                                                                     |                                                                                                                                                     
ggCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGa                                                                                                                                                       <  1:1975804/150‑1 (MQ=255)
 gCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAg                                                                                                                                                      <  2:515216/150‑1 (MQ=255)
 gCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAg                                                                                                                                                      <  2:2343598/150‑1 (MQ=255)
   cTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCt                                                                                                                                                    >  2:50007/1‑150 (MQ=255)
   cTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCt                                                                                                                                                    >  2:586874/1‑150 (MQ=255)
    tCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTg                                                                                                                                                   >  1:7416898/1‑150 (MQ=255)
     cGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCTCAGGGCTgg                                                                                                                                                  <  2:5976174/150‑1 (MQ=255)
      gggCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCTCAGGGCTGGc                                                                                                                                                 <  2:357275/150‑1 (MQ=255)
      gggCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGc                                                                                                                                                 >  1:6587595/1‑150 (MQ=255)
        gCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCg                                                                                                                                               <  1:7050473/150‑1 (MQ=255)
         cATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCgg                                                                                                                                              <  2:2970607/150‑1 (MQ=255)
         cATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCgg                                                                                                                                              <  2:4315403/150‑1 (MQ=255)
          aTGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCggg                                                                                                                                             >  1:7705102/1‑150 (MQ=255)
          aTGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCggg                                                                                                                                             <  2:520336/150‑1 (MQ=255)
              cggcggAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGtata                                                                                                                                         >  1:4014251/1‑150 (MQ=255)
               ggcggAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAg                                                                                                                                        >  2:1258297/1‑150 (MQ=255)
                gcggAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGc                                                                                                                                       <  1:1935718/150‑1 (MQ=255)
                gcggAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGc                                                                                                                                       <  2:5483561/150‑1 (MQ=255)
                  ggAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGt                                                                                                                                     <  2:8132352/150‑1 (MQ=255)
                        ggcgcggcgcGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCa                                                                                                                               <  1:1368689/150‑1 (MQ=255)
                        ggcgcggcgcGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCa                                                                                                                               <  1:8279626/150‑1 (MQ=255)
                           gcggcgcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGtt                                                                                                                            >  1:7034500/1‑150 (MQ=255)
                               cgcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGACCCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTg                                                                                                                        <  1:750959/150‑1 (MQ=255)
                                gcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTgg                                                                                                                       >  1:4349150/1‑150 (MQ=255)
                                gcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTgg                                                                                                                       <  1:1130867/150‑1 (MQ=255)
                                   tGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGgcca                                                                                                                    <  2:5396753/150‑1 (MQ=255)
                                    ggAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGgccag                                                                                                                   >  1:6754269/1‑150 (MQ=255)
                                              gAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTc                                                                                                         <  2:7866969/150‑1 (MQ=255)
                                                 cTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTg                                                                                                      <  1:6082570/150‑1 (MQ=255)
                                                 cTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTg                                                                                                      <  2:6346982/150‑1 (MQ=255)
                                                 cTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTg                                                                                                      >  2:291157/1‑150 (MQ=255)
                                                  tCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTGTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTgc                                                                                                     >  1:6341820/1‑150 (MQ=255)
                                                         aGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGa                                                                                              <  1:3298696/150‑1 (MQ=255)
                                                         aGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGa                                                                                              >  1:5786852/1‑150 (MQ=255)
                                                         aGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGa                                                                                              <  2:8296504/150‑1 (MQ=255)
                                                                accaccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGt                                                                                       >  2:7562235/1‑150 (MQ=255)
                                                                accaccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGt                                                                                       >  1:7283292/1‑150 (MQ=255)
                                                                 ccaccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTg                                                                                      >  2:1890385/1‑150 (MQ=255)
                                                                  caccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACACGCGTGc                                                                                     >  2:4457711/1‑150 (MQ=255)
                                                                  caccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGc                                                                                     <  1:2120321/150‑1 (MQ=255)
                                                                  caccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGc                                                                                     >  1:8800443/1‑150 (MQ=255)
                                                                   accGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGcc                                                                                    >  2:6704638/1‑150 (MQ=255)
                                                                   accGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGcc                                                                                    >  1:5069907/1‑150 (MQ=255)
                                                                    ccGAACCGGTCAGCCCTTTCCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGGCAGCCATTCCTGCGCTTCGACAAGCGTGCCt                                                                                   <  2:8051987/150‑1 (MQ=255)
                                                                    ccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCt                                                                                   >  2:167964/1‑150 (MQ=255)
                                                                      gAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTgg                                                                                 >  2:8615921/1‑150 (MQ=255)
                                                                      gAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTgg                                                                                 <  1:357615/150‑1 (MQ=255)
                                                                      gAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTgg                                                                                 >  1:384718/1‑150 (MQ=255)
                                                                       aaCCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTggg                                                                                >  1:1239853/1‑150 (MQ=255)
                                                                            gTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGggccggc                                                                           <  1:4902341/150‑1 (MQ=255)
                                                                            gTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGggccggc                                                                           >  2:2222094/1‑150 (MQ=255)
                                                                             tCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGggccggcc                                                                          >  1:3489831/1‑150 (MQ=255)
                                                                               aGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGccggccgg                                                                        <  2:526904/150‑1 (MQ=255)
                                                                                   cTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGggccgg                                                                            <  2:1620659/142‑1 (MQ=255)
                                                                                   cTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGggccgg                                                                            >  1:1620541/1‑142 (MQ=255)
                                                                                    ttACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGAt                                                                   >  1:4680431/1‑150 (MQ=255)
                                                                                      aCCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCg                                                                 <  2:2416586/150‑1 (MQ=255)
                                                                                        ccGCGGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCATGCCTGGGCCGGCCGGCAGATCGAg                                                               >  1:8006375/1‑150 (MQ=255)
                                                                                        ccGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAg                                                               >  1:7143251/1‑150 (MQ=255)
                                                                                         cGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGc                                                              >  1:8586174/1‑150 (MQ=255)
                                                                                          gCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCa                                                             <  2:2709262/150‑1 (MQ=255)
                                                                                             gAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGa                                                          >  2:5611239/1‑150 (MQ=255)
                                                                                                 cTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGgccagccagc                                                                                                              >  1:3414819/1‑94 (MQ=255)
                                                                                                    gACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTAcgc                                                   <  1:2524657/150‑1 (MQ=255)
                                                                                                       ttCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTAcgccgc                                                <  1:1548601/150‑1 (MQ=255)
                                                                                                       ttCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTAcgccgc                                                >  1:3638448/1‑150 (MQ=255)
                                                                                                          acacCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCcag                                             >  1:2313758/1‑150 (MQ=255)
                                                                                                            acCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCcagca                                           >  2:4009677/1‑150 (MQ=255)
                                                                                                              cGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAgcagg                                         >  1:2887830/1‑150 (MQ=255)
                                                                                                               gCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAgcaggt                                        >  1:5707246/1‑150 (MQ=255)
                                                                                                                 tGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAgcaggtgc                                      <  1:2620137/150‑1 (MQ=255)
                                                                                                                   ttCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAgcaggtgcag                                    <  1:175549/150‑1 (MQ=255)
                                                                                                                      gAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGcaggtgcaggtg                                 >  1:4598042/1‑150 (MQ=255)
                                                                                                                       aGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGa                                >  2:6204106/1‑150 (MQ=255)
                                                                                                                        ggACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAg                               >  2:6870338/1‑150 (MQ=255)
                                                                                                                         gACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGc                              >  1:8138843/1‑150 (MQ=255)
                                                                                                                              tttCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACt                         <  1:5528408/150‑1 (MQ=255)
                                                                                                                                tCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGa                       <  1:7889332/150‑1 (MQ=255)
                                                                                                                                 cTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGAt                      <  2:1987732/150‑1 (MQ=255)
                                                                                                                                 cTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGAt                      <  2:3466280/150‑1 (MQ=255)
                                                                                                                                 cTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGAt                      <  1:3799075/150‑1 (MQ=255)
                                                                                                                                  tGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATg                     >  2:4239826/1‑150 (MQ=255)
                                                                                                                                   gggCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATgg                    >  2:8773713/1‑150 (MQ=255)
                                                                                                                                    ggCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGTAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGa                   <  1:50007/150‑1 (MQ=255)
                                                                                                                                       cccggCCAGCTCAGGGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACt                >  1:3885002/5‑150 (MQ=255)
                                                                                                                                       cATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACt                <  2:3771702/150‑1 (MQ=255)
                                                                                                                                        ccggCCAGCTCAGGGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTACGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGAGCAGGTGAGCGAACTGATGGAACTg               >  1:3192735/4‑150 (MQ=255)
                                                                                                                                          ggCCAGCTCAGGGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCCGGTGAGCGAACTGATGGAACTGGa             >  2:5561815/2‑150 (MQ=255)
                                                                                                                                          ggCCAGCTCAGGGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGa             >  1:6478062/2‑150 (MQ=255)
                                                                                                                                           gCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGAc            <  2:735506/150‑1 (MQ=255)
                                                                                                                                           aatggCCCAGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGAc            >  2:3335637/5‑150 (MQ=255)
                                                                                                                                            ccAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACa           >  2:341395/1‑150 (MQ=255)
                                                                                                                                            ccAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACa           >  2:2582884/1‑150 (MQ=255)
                                                                                                                                            ccAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACa           >  2:7114346/1‑150 (MQ=255)
                                                                                                                                              aGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGACCTGGACAGc         >  2:6617813/1‑150 (MQ=255)
                                                                                                                                                cccTGAGCGGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGACCTGATGGAACTGGACAGCCa       <  2:4879579/150‑1 (MQ=255)
                                                                                                                                                 ccTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTATGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGCCAg      <  2:7412140/150‑1 (MQ=255)
                                                                                                                                                 ccTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGCCAg      <  2:6095427/150‑1 (MQ=255)
                                                                                                                                                 ccTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGCCAg      <  2:1191383/150‑1 (MQ=255)
                                                                                                                                                  cTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGCCAgg     >  1:4195534/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGCCAGGAAG  >  NC_002947/4044505‑4044803

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: