Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 962,272 Δ1 bp 5.4% coding (146/1035 nt) selD ← selenide, water dikinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947962,2720T. 5.4% 338.6 / 13.3 111coding (146/1035 nt)selDselenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base T (48/56);  new base . (2/4);  total (51/60)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

Only 100 of 115 total aligned reads displayed.

GCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCATACGC  >  NC_002947/962123‑962421
                                                                                                                                                     |                                                                                                                                                        
gCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCAtcgt                                                                                                                                                          >  1:90795/1‑150 (MQ=255)
         cgccCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGGCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATg                                                                                                                                                >  2:260127/1‑150 (MQ=255)
           cccATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCg                                                                                                                                              >  2:746914/1‑150 (MQ=255)
           cccATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGCGCGACGATGcc                                                                                                                                              <  1:3005109/150‑2 (MQ=255)
            ccATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGCGCGACGATGcag                                                                                                                                             <  1:1283951/150‑3 (MQ=255)
                  cATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGt                                                                                                                                             >  1:2407909/1‑144 (MQ=255)
                  cATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGt                                                                                                                                             <  2:2407909/144‑1 (MQ=255)
                   aTAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGAc                                                                                                                                      <  1:1974573/150‑1 (MQ=255)
                   aTAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGAc                                                                                                                                      >  1:1824566/1‑150 (MQ=255)
                        tGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTTCCGACCcaca                                                                                                                                 >  2:1246750/1‑150 (MQ=255)
                              tGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCg                                                                                                                           <  1:3000042/150‑1 (MQ=255)
                               gATGGCGTTGGTGGCGGCCATGCGGCCCCGGCCGTCCGGGTGGATGCCGTCATGCATGAAGAAGGCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                          <  2:674851/150‑1 (MQ=255)
                               gATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                          <  1:2737924/150‑1 (MQ=255)
                                aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCggg                                                                                                                         <  2:2847540/150‑1 (MQ=255)
                                aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCggg                                                                                                                         <  2:2912584/150‑1 (MQ=255)
                                     gTTGGTGGCGGCAATGCGGCCGAAGTCATAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGCGCGACGATGca                                                                                                                                              >  2:1070528/1‑123 (MQ=255)
                                     gTTGGTGGCGGCAATGCGGCCGAAGTCATAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGCGCGACGATGca                                                                                                                                              <  1:1070528/124‑2 (MQ=255)
                                      ttGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAgg                                                                                                                   <  1:214701/150‑1 (MQ=255)
                                      ttGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAgg                                                                                                                   >  2:1685021/1‑150 (MQ=255)
                                      ttGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAgg                                                                                                                   <  1:260127/150‑1 (MQ=255)
                                         gtggCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGcc                                                                                                                >  1:1667536/1‑150 (MQ=255)
                                           ggcggcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTAGGGTCCAGGGCCTg                                                                                                              <  2:559559/150‑1 (MQ=255)
                                              ggcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGt                                                                                                           <  2:1824566/150‑1 (MQ=255)
                                               gcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTg                                                                                                          >  1:588231/1‑150 (MQ=255)
                                                 aaTGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGcc                                                                                                        <  2:1080925/150‑1 (MQ=255)
                                                  aTGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCg                                                                                                       <  1:542384/150‑1 (MQ=255)
                                                         cGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCgg                                                                                                <  1:624221/150‑1 (MQ=255)
                                                          gAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGACCTgg                                                                                                             <  1:2327571/136‑1 (MQ=255)
                                                          gAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCACGATGCGTTGCCGACCCACCGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGc                                                                                               >  2:2082584/1‑150 (MQ=255)
                                                           aaGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCa                                                                                              >  1:1292380/1‑150 (MQ=255)
                                                             gTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAg                                                                                            <  1:2166954/150‑1 (MQ=255)
                                                             gTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAg                                                                                            <  1:1685021/150‑1 (MQ=255)
                                                              tCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAgg                                                                                           >  1:2576620/1‑150 (MQ=255)
                                                               cGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGa                                                                                          <  2:1303183/150‑1 (MQ=255)
                                                                gTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTTCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGAt                                                                                         >  2:914671/1‑150 (MQ=255)
                                                                gTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGAt                                                                                         >  1:2897251/1‑150 (MQ=255)
                                                                   aGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCAc                                                                                      <  2:1667536/150‑1 (MQ=255)
                                                                   aGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGCGCGACGATGcagca                                                                                                                                           >  2:2253515/1‑93 (MQ=255)
                                                                    gggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACg                                                                                     >  2:1133757/1‑150 (MQ=255)
                                                                     ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                          <  2:435807/112‑1 (MQ=255)
                                                                     ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                          >  1:435807/1‑112 (MQ=255)
                                                                     ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGt                                                                                    <  2:1206500/150‑1 (MQ=255)
                                                                     ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGt                                                                                    >  1:236766/1‑150 (MQ=255)
                                                                       tCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTcc                                                                                  >  2:1162145/1‑150 (MQ=255)
                                                                            gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTgg                                                                                                             >  2:247897/1‑118 (MQ=255)
                                                                            gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCAc                                                                                      <  2:2759920/141‑1 (MQ=255)
                                                                            gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCAc                                                                                      >  1:2759920/1‑141 (MQ=255)
                                                                            gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGt                                                                                    <  2:2240294/143‑1 (MQ=255)
                                                                            gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGt                                                                                    >  1:2240294/1‑143 (MQ=255)
                                                                            gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCAc                                                                             >  1:776677/1‑150 (MQ=255)
                                                                               gaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCt                                                                                                               <  1:2377645/113‑1 (MQ=255)
                                                                               gaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCt                                                                                                               >  2:2377645/1‑113 (MQ=255)
                                                                               gaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCg                                                                                                 >  2:2619989/1‑127 (MQ=255)
                                                                               gaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCg                                                                                                 <  1:2619989/127‑1 (MQ=255)
                                                                               gaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCtt                                                                          >  2:868404/1‑150 (MQ=255)
                                                                                aTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg                                                                         <  2:90795/150‑1 (MQ=255)
                                                                                 tCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgg                                                                       <  2:222841/150‑1 (MQ=255)
                                                                                    gCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTggggg                                                                    >  2:1564489/1‑150 (MQ=255)
                                                                                      aTGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGGCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGACAAGGATCACGTCCAGCACATTCGGGgag                                                                 >  2:1253955/1‑150 (MQ=255)
                                                                                       tGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgg                                                                       >  1:11608/1‑144 (MQ=255)
                                                                                        gaagaaGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGAt                                                               <  2:1285602/150‑1 (MQ=255)
                                                                                          agaagTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGCAGATCt                                                             >  2:291808/1‑150 (MQ=255)
                                                                                           gaagTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCtt                                                            >  2:1665966/1‑150 (MQ=255)
                                                                                           gaagTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCtt                                                            <  1:1564489/150‑1 (MQ=255)
                                                                                              gTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                          >  2:2822891/1‑87 (MQ=255)
                                                                                              gTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCCTGGGGGAGATCTTGCa                                                         >  1:1668912/1‑150 (MQ=255)
                                                                                               tCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGAATCTTGGGGGAGATCTTGCaa                                                        <  2:2576620/150‑1 (MQ=255)
                                                                                                cGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCAACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAAc                                                       <  2:513204/150‑1 (MQ=255)
                                                                                                      tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTggg                                                                                                            >  2:1043658/1‑93 (MQ=255)
                                                                                                      tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTggg                                                                                                            <  1:1043658/93‑1 (MQ=255)
                                                                                                      tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGc                                                 >  1:2639262/1‑150 (MQ=255)
                                                                                                              aCGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCa                                         <  2:236766/150‑1 (MQ=255)
                                                                                                                 ccgcGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATgg                                      <  1:1526437/150‑1 (MQ=255)
                                                                                                                  cgcgGTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCTCGATGCTTTGCCGATCCACTGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATGACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGc                                     <  1:2082584/150‑1 (MQ=255)
                                                                                                                     gcTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGc                                                 <  2:2344271/135‑1 (MQ=255)
                                                                                                                     gcTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGc                                                 >  1:2344271/1‑135 (MQ=255)
                                                                                                                         gtcgtcCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                                                     >  1:71987/1‑112 (MQ=255)
                                                                                                                         gtcgtcCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                                                     <  2:71987/112‑1 (MQ=255)
                                                                                                                          tcgtcCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTgg                             >  1:983745/1‑150 (MQ=255)
                                                                                                                          tcgtcCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTgg                             <  2:1635518/150‑1 (MQ=255)
                                                                                                                               aaGCGCCTACACTGCTTCATAGTCGCGCGATGCGTTGCCGACCCACGTTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAg                        <  1:1427655/149‑1 (MQ=255)
                                                                                                                                 gcgcGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGc                      <  2:3068981/150‑1 (MQ=255)
                                                                                                                                       acacTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCg                <  1:1162145/150‑1 (MQ=255)
                                                                                                                                       acacTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCg                <  2:3038234/150‑1 (MQ=255)
                                                                                                                                       acacTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCg                >  1:2997869/1‑150 (MQ=255)
                                                                                                                                         acTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgc              >  1:2322343/1‑150 (MQ=255)
                                                                                                                                           tgctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCaa                                                        <  2:2753932/106‑1 (MQ=255)
                                                                                                                                           tgctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCaa                                                        >  1:2753932/1‑106 (MQ=255)
                                                                                                                                           tgctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGc                      <  1:1518958/140‑1 (MQ=255)
                                                                                                                                           tgctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgctc            >  2:1768594/1‑150 (MQ=255)
                                                                                                                                           tgctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgctc            >  1:511018/1‑150 (MQ=255)
                                                                                                                                           tgctgcATCGTCGCGAGCTGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGc                      >  2:1518958/1‑140 (MQ=255)
                                                                                                                                             ctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCCTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgctcgc          <  2:1668912/150‑1 (MQ=255)
                                                                                                                                              tgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTg                                                                                                              >  2:2262999/1‑51 (MQ=255)
                                                                                                                                              tgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTg                                                                                                              <  1:2262999/51‑1 (MQ=255)
                                                                                                                                              tgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgctcgct         >  2:2329382/1‑150 (MQ=255)
                                                                                                                                                cATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCa       <  1:914671/150‑1 (MQ=255)
                                                                                                                                                 aTCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGT‑‑AACGCTTTGCAGGCATTAACTGGCAACTACAGCCGGCGCCCCGGCTGTCTTAGGTCAGGCGAATCGGCTCGCTCAt      >  1:1595585/1‑150 (MQ=255)
                                                                                                                                                    gtcgCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCATACg   >  1:425632/1‑150 (MQ=255)
                                                                                                                                                     tcgCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCATACGc  >  1:2823314/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
GCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCATACGC  >  NC_002947/962123‑962421

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