Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,937,802 A→G D112G (GAC→GGC)  fucK → L‑fuculokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,937,8020AG98.0% 336.8 / ‑3.1 98D112G (GAC→GGC) fucKL‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (46/50);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

GTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCT  >  NC_000913/2937667‑2937936
                                                                                                                                       |                                                                                                                                      
gTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCaa                                                                                                                                     <  1:156275/139‑1 (MQ=255)
gTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCaa                                                                                                                                     <  1:624242/139‑1 (MQ=255)
    actgactgAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATAtt                                                                                                                                 >  1:236886/1‑139 (MQ=255)
        actgaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAAc                                                                                                                             >  2:1455021/1‑139 (MQ=255)
                 cacaTCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAAtct                                                                                                                    >  1:837450/1‑139 (MQ=255)
                         cgGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAGATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAg                                                                                                            >  1:901702/1‑139 (MQ=255)
                         cgGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAg                                                                                                            >  1:546731/1‑139 (MQ=255)
                              tCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGtt                                                                                                       >  2:603438/1‑139 (MQ=255)
                               cGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      <  2:837450/139‑1 (MQ=255)
                               cGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      <  1:727925/139‑1 (MQ=255)
                               cGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      <  2:392070/139‑1 (MQ=255)
                               cGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGAGAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      <  2:901702/139‑1 (MQ=255)
                                   gTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAgg                                                                                                  <  2:979782/139‑1 (MQ=255)
                                      accaccTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCa                                                                                                                                      <  2:286895/100‑1 (MQ=255)
                                      accaccTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCa                                                                                                                                      >  1:286895/1‑100 (MQ=255)
                                                    ggATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGa                                                                                 <  1:304290/139‑1 (MQ=255)
                                                     gATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAg                                                                                >  1:729829/1‑139 (MQ=255)
                                                          cgcTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCg                                                                                   >  1:817526/1‑131 (MQ=255)
                                                          cgcTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCg                                                                                   <  2:817526/131‑1 (MQ=255)
                                                          cgcTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCttt                                                                           <  1:807581/139‑1 (MQ=255)
                                                             tctGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGt                                                                        <  2:236886/139‑1 (MQ=255)
                                                               tGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGtt                                                                                                                         <  2:1148524/88‑1 (MQ=255)
                                                               tGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCa                                                                                                    <  2:249896/109‑1 (MQ=255)
                                                               tGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCa                                                                                                    >  1:249896/1‑109 (MQ=255)
                                                               tGGGAGATAAGCAAGGCAATGTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGtt                                                                                                                         >  1:1148524/1‑88 (MQ=255)
                                                                 gTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                                    <  1:603438/139‑1 (MQ=255)
                                                                     aTAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGtt                                                                                                                         >  1:175081/1‑82 (MQ=255)
                                                                     aTAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGtt                                                                                                                         <  2:175081/82‑1 (MQ=255)
                                                                         gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAg                                                                                                            >  2:1509770/1‑91 (MQ=255)
                                                                         gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAg                                                                                                            <  1:1509770/91‑1 (MQ=255)
                                                                         gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                               <  2:974026/136‑1 (MQ=255)
                                                                         gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                               >  1:974026/1‑136 (MQ=255)
                                                                           aaGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTtata                                                          >  2:685602/1‑139 (MQ=255)
                                                                              gCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATaa                                                       >  2:1420440/1‑139 (MQ=255)
                                                                                  tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCa                                                                                                    <  1:879627/90‑1 (MQ=255)
                                                                                  tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCa                                                                                                    >  2:879627/1‑90 (MQ=255)
                                                                                    tGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      >  2:815280/1‑86 (MQ=255)
                                                                                    tGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      <  1:815280/86‑1 (MQ=255)
                                                                                           tCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTg                                          >  1:709140/1‑139 (MQ=255)
                                                                                            ccGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATAtttac                                                                                                                              >  1:1167177/1‑51 (MQ=255)
                                                                                            ccGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGaa                                                                                                                              <  2:1167177/54‑1 (MQ=255)
                                                                                             cGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTgaa                                        <  1:685602/139‑1 (MQ=255)
                                                                                                ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGgtg                                                >  1:464483/1‑128 (MQ=255)
                                                                                                ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGgtg                                                <  2:464483/128‑1 (MQ=255)
                                                                                                ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTgaaaga                                     >  2:1104741/1‑139 (MQ=255)
                                                                                                ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTgaaaga                                     >  1:1510293/1‑139 (MQ=255)
                                                                                                   ttaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGaaaa                                  <  2:729829/139‑1 (MQ=255)
                                                                                                     aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGtt                                                                                                                         >  1:1687845/1‑50 (MQ=255)
                                                                                                     aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGtt                                                                                                                         <  2:1687845/50‑1 (MQ=255)
                                                                                                     aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                                    <  2:90873/103‑1 (MQ=255)
                                                                                                     aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                                    >  1:90873/1‑103 (MQ=255)
                                                                                                     aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                                 <  1:1148619/106‑1 (MQ=255)
                                                                                                     aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                                 >  2:1148619/1‑106 (MQ=255)
                                                                                                     aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGaaaa                                  >  2:260798/1‑137 (MQ=255)
                                                                                                     aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGaaaa                                  <  1:260798/137‑1 (MQ=255)
                                                                                                     aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAatc                                >  1:462535/1‑139 (MQ=255)
                                                                                                        tGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAatcatc                             <  2:709140/139‑1 (MQ=255)
                                                                                                         ggAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGc                                                                                                                <  1:23849/55‑1 (MQ=255)
                                                                                                         ggAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGc                                                                                                                >  2:23849/1‑55 (MQ=255)
                                                                                                         ggAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTtat                                                           <  1:211795/108‑1 (MQ=255)
                                                                                                         ggAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTtat                                                           >  2:211795/1‑108 (MQ=255)
                                                                                                          gAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAg                                                                                                            >  2:244064/1‑58 (MQ=255)
                                                                                                          gAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAg                                                                                                            <  1:244064/58‑1 (MQ=255)
                                                                                                          gAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      >  1:960154/1‑64 (MQ=255)
                                                                                                          gAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      <  2:960154/64‑1 (MQ=255)
                                                                                                           aaaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTaa                                                                                                                       >  1:296919/1‑46 (MQ=255)
                                                                                                           aaaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTaa                                                                                                                       <  2:296919/46‑1 (MQ=255)
                                                                                                           aaaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCcac                          <  2:546731/139‑1 (MQ=255)
                                                                                                           aaaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCcac                          <  2:462535/139‑1 (MQ=255)
                                                                                                           aaaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCcac                          <  1:1420440/139‑1 (MQ=255)
                                                                                                            aaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCcaca                         >  1:1174321/1‑139 (MQ=255)
                                                                                                             aTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCa                                                                                                    >  2:417509/1‑63 (MQ=255)
                                                                                                             aTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCa                                                                                                    <  1:417509/63‑1 (MQ=255)
                                                                                                             aTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACaa                        <  2:1174321/139‑1 (MQ=255)
                                                                                                             aTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACaa                        <  1:1455021/139‑1 (MQ=255)
                                                                                                                tCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                               <  2:1296904/97‑1 (MQ=255)
                                                                                                                tCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                               >  1:1296904/1‑97 (MQ=255)
                                                                                                                tCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAActg                     >  1:842686/1‑139 (MQ=255)
                                                                                                                tCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAg                                                                                                   >  2:1164382/1‑61 (MQ=255)
                                                                                                                tCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAg                                                                                                   <  1:1164382/61‑1 (MQ=255)
                                                                                                                 ccGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATAc                                                                <  1:1079001/95‑1 (MQ=255)
                                                                                                                 ccGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATAc                                                                >  2:1079001/1‑95 (MQ=255)
                                                                                                                 ccGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGt                                            >  2:764517/1‑115 (MQ=255)
                                                                                                                 ccGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGt                                            <  1:764517/115‑1 (MQ=255)
                                                                                                                 ccGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAActgc                    >  2:1476222/1‑139 (MQ=255)
                                                                                                                  cgcgAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      <  2:943246/56‑1 (MQ=255)
                                                                                                                  cgcgAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                      >  1:943246/1‑56 (MQ=255)
                                                                                                                    cgAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                                    >  2:898885/1‑88 (MQ=255)
                                                                                                                    cgAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                                    <  1:898885/88‑1 (MQ=255)
                                                                                                                    cgAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTgaaaga                                     <  2:611582/119‑1 (MQ=255)
                                                                                                                    cgAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTgaaaga                                     >  1:611582/1‑119 (MQ=255)
                                                                                                                      aaCAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATACAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGaa               <  1:1104741/139‑1 (MQ=255)
                                                                                                                         agcagcGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                                    >  2:105806/1‑83 (MQ=255)
                                                                                                                         agcagcGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                                    <  1:105806/83‑1 (MQ=255)
                                                                                                                            agcGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAAcgcgcg         <  1:1476222/139‑1 (MQ=255)
                                                                                                                                   aTGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTtat                                                           <  2:1395884/82‑1 (MQ=255)
                                                                                                                                   aTGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTtat                                                           >  1:1395884/1‑82 (MQ=255)
                                                                                                                                   aTGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCt  >  1:194650/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCT  >  NC_000913/2937667‑2937936

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